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Uber freie Biegungsschwingungen des axial belasteten Stabes 
mit innerer und duferer Dimpfung. 


Von H. Eschler. 


1, Aufgabenstellung. Das Problem der freien Biegungsschwingungen eines Stabes mit 
innerer Dampfung wurde von K. Muto! und K.Sezawa? behandelt. Beide Verfasser legten 
ihren Untersuchungen im Anschlu8 an E. J. Routh? die Annahme zugrunde, da sich die Damp- 


fungsspannung in der Form 
0g = ines (1) 


ausdriicken laBt, wobei FE den Elastizitatsmodul, r; den Koeffizienten der inneren Reibung, der 
eine reine Materialkonstante ist*, und d¢/dt die Deformationsgeschwindigkeit bedeuten. 

Der Routhsche Ansatz wurde von R. Giovannozzi® anlaBlich der Untersuchung von erzwun- 
genen Stabschwingungen etwas modifiziert. Wenn w die Kreisfrequenz der pulsierenden auBeren 
Kraft und damit also auch der erzwungenen Schwingungen ist, wird fiir die Dampfungsspannung 

k de 
angenommen. Der Koeffizient k hangt auBer von der chemischen Zusammensetzung des Werk- 
stoffes auch von der Schwingweite der Deformation, der Temperatur des Stabes und unter 
Umstanden von der Frequenz ab. Bei der Integration der Differentialgleichung hat Giovannozxi 
fiir k einen konstanten mittleren Wert gesetzt. Wenn man den Ansatz (2) auch bei den freien 
Biegungsschwingungen von Staben benutzen will, so wird w offenbar als die noch unbekannte 
Kreisfrequenz der jeweils betrachteten Eigenschwingung aufzufassen sein. 

Wie man zwei experimentellen Arbeiten von Giovannozzi® entnehmen kann, ist der Ansatz (2) 
allgemein anwendbar, wahrend (1) anscheinend nur bei kleinen Schwingweiten angenahert gilt. 
Um im folgenden die Verwendung beider Ansatze zu erméglichen, schreiben wir 


0g ITA Se (3) 


mit K=r; geht (3) in (1) und mit K= k/w in (2) iber. 

Zur Berechnung des Gliedes, um welches die Differentialgleichung der Biegungsschwingungen 
des Stabes bei Beriicksichtigung der inneren Dampfung zu vermehren ist, bedenken wir, da 
fiir die Dehnung ¢ und die Spannung o im Abstand z von der Nullinie die Beziehungen 

Oy 2. 
—— aa und. O—2 an 
gelten, wenn mit y(t,x) die Auslenkung der Stabachse zur Zeit t an der Stelle x, mit M das 
Biegungsmoment und mit J das Tragheitsmoment des Querschnittes bezeichnet wird. Durch 
Einsetzen in (3) ergibt sich das durch die innere Dampfung hervorgerufene Moment 


oe 


a Py 
M,=—EJK>e5 ’ 


aus welchem der die innere Dampfung beriicksichtigende Term in der Differentialgleichung 
durch zweimalige Differentiation nach x folgt. 


1 K. Muto, Z. angew. Math. Mech. 10 (1930) S. 346. 
2 K. Sezawa, Z. angew. Math. Mech. 12 (1932) S. 275. 2 
3 E. J. Routh, Die Dynamik der Systeme starrer Kérper, Bd. II, S. 491. Leipzig 1898. 
4 Werte von € = Eri fiir einige Werkstoffe finden sich in der unter Fufnote 2 zitierten Arbeit auf 
S$. 279 
5 R. Giovannozzi, Aerotecnica 17 (1937) S. 1047. 
6 R. Giovannozzi, Aerotecnica 18 (1938) und 19 (1939) S. 245. 
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Die AuBere Dampfung kann proportional der Geschwindigkeit der Stabachse gesetzt werden}, 
wobei ein Koeffizient 7, der Gestalt des Stabes und dem umgebenden Medium Rechnung tragt. 

Bei Beriicksichtigung der inneren und d4uBeren Dimpfung sowie der Rotationstragheit lautet 
dann die Differentialgleichung der Biegungsschwingungen eines Stabes, der unter dem EinfluB 
einer axialen Druckkraft P steht, 

2 4 5 4 : 
oF 8Y 4 EF PY 4 BIK PY, + 1 2 ol sede 4 eee iy (4) 

Diese Differentialgleichung ist nur fiir den auf beiden Enden frei aufliegenden Stab ohne weiteres 
integrierbar, wobei der den Randbedingungen geniigende Ansatz y= q(t) sin nzx/l, mit 1 als 
Stablange, die Lésung vermittelt. Bei Vorliegen von anderen Randbedingungen wird die Inte- 
gration von (4) schwierig. Man itberzeugt sich leicht, daB® der Bernoullische Ansatz y= q(t) X(x) 
nur in folgenden Sonderfallen zum Ziel fiihrt: 

a) bei Vernachlassigung der inneren und duBeren Dampfung, 

b) bei Vernachlassigung der inneren Dampfung und Rotationstragheit, 

c) bei fehlender Langskraft und Vernachliassigung der aufSeren Dampfung, 

d) bei fehlender Langskraft und Vernachlassigung der Rotationstragheit. 
Dieser Sachverhalt 1aBt aber anderseits vermuten, daB es auch im allgemeinen Fall sehr ange- 
nahert zu stehenden Schwingungen kommt. Eine weitere Stiitze fiir diese Annahme liefert das 
Experiment. In der Wirklichkeit wird man mit innerer und auBerer Dampfung sowie mit dem 
Einflu®B der Rotationstragheit immer rechnen miissen. Obwohl die Differentialgleichung in 
diesem Fall durch den Bernoullischen Ansatz nicht gelést werden kann, ergibt der Versuch doch 
praktisch stehende Schwingungen. Zur Auffindung einer Naherungslésung im allgemeinen 
Fall kann daher gewif ein Produktansatz 


Yn = In (t) Xn (2) (5) 
fiir die n-te Eigenschwingung verwendet werden. 

In erster Linie werden die Eigenschwingungszahlen und Dampfungsexponenten interessieren. 
Solange Reibungskrafte nicht beriicksichtigt werden, erméglicht bekanntlich das Rayleighsche 
Verfahren, welches sich darauf griindet, da in der Energiegleichung die Eigenfunktionen durch 
benachbarte Funktionen ersetzt werden, auf einfache Weise die naherungsweise Berechnung der 
Eigenschwingungszahlen sowohl des Grundtones als auch der Oberténe. Im folgenden soll nun 
gezeigt werden, wie sich auch in unserem Fall eine Naherungsliésung auf analoge Art auffinden 
1aBt. Der dabei beschrittene Weg kann als eine Erweiterung des Rayleighschen Verfahrens auf 
Systeme, die unter der Wirkung von Reibungskraften stehen, bezeichnet werden. 


2, Naherungs!ésung. Den Ausgangspunkt fiir unser Vorhaben bilden die Lagrangeschen 
Gleichungen in der erweiterten Form 
cyOT nO le claw eC 
x\34,) cca RANE Tin emigre 8) (6) 


n 


worin 7’ die kinetische, U die potentielle Energie und Z die Zerstreuungsfunktion, welche gleich 
der halben Arbeit der Reibungskrafte je Zeiteinheit zu setzen ist, bedeuten. Die der Gleichung (6) 
zugrunde liegende Voraussetzung, daf die Reibungskrafte den Geschwindigkeiten direkt pro- 
portional sein miissen, ist fiir die 4uhere Dampfung erfillt und trifft mit dem Ansatz (5) auch 
fiir die innere Dampfung zu, da sich 


IV 
By, Fn ~*n 


Ot 0x4 Ge XX, 


ergibt. : 
Die kinetische Energie mit Beriicksichtigung der Rotationstrigheit ist durch den Ausdruck 


d igo 
Pees is | (By ax+ Fes | (srae) ae 
0 a 


gegeben, welcher fiir die n-te Higenschwingung mit (5) iibergeht in 
l l 


1 a 1 mabe eee 
T=FeP a | Xtde+t oss [ xe ax. (7) 
0 0 


1 Vel. FuBnote ISvons >. el. 
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Vermindert man die Biegungsenergie um die Energie der Zusammendriickung durch die Axial- 
druckkraft P, so ergibt sich die Bezichung fiir die potentielle Energie 
l l 


U=LEJ [ (2) ax ei:p f(y dx 


0 0 
oder fiir die n-te Eigenschwingung 
; 1 i 
1 


2 7 112, 1 2 79 
Ua aE! a. | Xp dx— + P qi | das, (8) 
ft 


6 

Die Zerstreuungsfunktion Z setzt sich aus dep Anteilen Z; fiir die innere und Z, fiir die AauBere 
Dampfung zusammen. Um die Beziehung fiir Z; zu finden, bedenken wir, daB® die von der 
Normalspannung o am Einheitswiirfel bei einer VergréBerung der Dehnung ¢ auf ¢-|- de geleistete 
Arbeit ode betragt. Dann ergibt sich mit dV als Volumelement nach obiger Definition fiir die 


Zerstreuungsfunktion 
APY tT Sha Tal Oa mrss 
=z] | forge a’. 


Setzt man den Ausdruck fiir die Dampfungsspannung (3) ein, so folgt mit Hilfe der schon frither 
benutzten Beziehung ¢—=— z d*y/dx? sowie wegen [ede J weiter 


1 
raha LO p NE 
ee Ks [ (ede dar 


woraus man nach Ejinsetzen von (5) 
l 
Z,= = EIK¢ [ xa dx (9) 
0 


erhalt. 
Fur den Anteil Z, ergibt sich, da vom Langenelement dx des Stabes langs des Weges dy die 
Reibungsarbeit r, dy/dt dx Oy zu leisten ist, definitionsgemaB 
l 


Berl ay\2 
Tig Ta [ (=e) dx, 
0 


und fiir die n-te Eigenschwingung wird daraus 
l 


rn | X2dx. (10) 
0 

Einfiihren von (7), (8), (9) und (10) in die Lagrangesche Gleichung (6) liefert dann mit den 

Abkiirzungen 


1 
l= > 


1 I I 
Sy= | X2dx, S| Xe da, Se es (11) 
0 0 0 


fiir q, die Differentialgleichung 
@ (FS) + JS1) In + (ET KS. + ta Sp) dn + (EIS, — PSs) dn = 0. 


Vermittels des Ansatzes 9, = ef—*n + *n)* ergeben sich daraus fiir die Eigenfrequenzen w, und 


die Diampfungsexponenten 6, die Beziehungen 
es IS PS) /HIKS, tras? 
On O(FSo+ IS) \2e(FS) + IS) 
SPEIRS: 108; 
2 e(FS)+ JS) 
Bei Verwendung des Ansatzes (2) ist zu beachten, da K noch Funktion von w, ist. 


Diese Formeln erméglichen nun tatsachlich die naherungsweise Berechnung der Kigenschwin- 
gungszahlen w, und der Grifen 6, in denjenigen Fallen, in welchen die strenge Liésung der Diffe- 


1* 


(12) 


On (n=1,2,3,...) | 
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rentialgleichung (4) nicht angegeben werden kann. Man braucht nur, analog wie beim Rayleigh- 
schen Verfahren, in den Integralen (11) die Schwingungsformen des Stabes durch benachbarte 
Funktionen zu ersetzen. Das Verfahren, das sich iibrigens auch ohne weiteres bei anderen Rand- 
wertproblemen der Elastokinetik, wie z.B. der schwingenden Platte, anwenden 1a48t, soll an 
Hand eines Beispieles noch niaher erlautert werden. 

Zuvor wollen wir noch kurz auf den beiderseits frei aufliegenden Stab zuriickkommen. Mit 
den Eigenfunktionen X,— sin nzx/lI lassen sich die Integrale (11) leicht berechnen und aus (12) 
folgt 


2 


On —G 2 (BE nn? J) 


— ntat EJK + ral! @ 
~ 202 (PF+ nn? J) 


n2at? (n2? EJ — I? P) nint EJK + ralt \2 
( o2(PF+n'xJ)) ’ 


On 


Dieses Ergebnis, das mit dem Ansatz y,—= ef nt FOr)! cin nzx/l auch unmittelbar aus der 


Differentialgleichung (4) erhalten werden kann, stellt in diesem Sonderfall schon die strenge 
Lésung dar. 
Fir K=r,= 0 vereinfacht sich (12) auf 
Es =p 
GOLe 

" @(FS,+JS;) ’ 
eine Beziehung, die auch aus der Energiegleichung folgt und den Ausgangspunkt fiir das Rayleigh- 
sche Verfahren bildet. 

3. Beispiel. Wir behandeln als Beispiel fiir das angegebene Naherungsverfahren den an 
beiden Enden eingespannten Stab. Als benachbarte Funktionen X, kénnen die Eigenfunktionen 
des Stabes bei vernachlassigter Rotationstragheit und fehlender Langskraft verwendet werden. 
Es sind dies die Integrale der Differentialgleichung 


ax 
dx* 
welche aus (4) mit dem Bernoullischen Ansatz leicht erhalten wird. Bemerkenswerterweise 
besitzt namlich der rotationstragheitsfreie, axial unbelastete Stab mit und ohne Beriicksichtigung 
der inneren sowie auferen Dampfung die gleichen Eigenfunktionen. 

Die Rechnung vereinfacht sich, wenn man die symmetrischen und antisymmetrischen Schwin- 
gungen gesondert betrachtet. Die Funktionen X,, fiir die ersteren lauten im Falle eines bei x—0 
und x= eingespannten Stabes 


EJ —RRoFX=0, 


ey en! l gnl l 
X,= Co} “> COS Pn (* 3) cos Co] Pn (x = z) > (13) 
Se ee l 1 
worin die ¢,/ aus tg ses + Sg ee = 0 zu bestimmen sind. Fiir die antisymmetrischen Schwin- 
gungen gilt 
— ee Qnl és l . gnl ~: l 
N= Olt “9 S10 Qn (* x] ies Sin Pn (-— 3 (14) 


mit der Periodengleichung tg oe — fq a ==A(), 


Setzt man (13) bzw. (14) in die Integrale (11) ein, so ergibt sich fiir die symmetrischen 
Schwingungen 


I l 
Ge cost pit 
h cos? Coj ae 


Fall mae n 
5S, = Vn E (sor? a — cos? ae 2 sin © 


lh nl Qnl 
9 cos Co}? 


2 2 2 


und fiir die antisymmetrischen Schwingungen 


ito Gall parame 
1 q Yn ~ 2 in 
So = [sin ae Sin a 


— gnl ~: nl 5 all a Pl nl 2 ve 
51 = $n 5} (Sine + sin? 5 ) 2sin S cos 3 Sin? | : 


XX. Band 1952. Eschler: Uber freie Biegungsschwingungen des axial belasteten Stabes usw. 5) 


wahrend fiir beide Schwingungsformen 

SG 
ist. Einsetzen dieser Ausdriicke in (12) liefert mit Benutzung der Periodengleichungen die ge- 
suchten Beziehungen 


bie == (gnl)4 EID; iP. ( (gnl)4 EJK + Tal \2 
a oP (PEF + GJ) | 20 P(E + DJ) , | 
a (Qnl)4 EJK+ rau lt 

Se Pon (eh OD, J) 


(15) 


(eed, 3.04.) 4 | 


worin zur Abkiirzung geschrieben wurde 


nl 
pnlg 2 (gal Xg 


nl 3 : 
ae a 2) fir die symmetrischen Schwingungen, 


oo 


nl nl : 
Pnl Ctg e (aio e — 2) fiir die antisymmetrischen Schwingungen. 


In der folgenden Tabelle sind die ersten fiinf Werte von grt und @®, zusammen- 
gestellt. Fiir die symmetrischen Schwingungen ist n ungerade, fiir die antisymmetrischen 
gerade. 


n 1 | 2 | | 4 5 
gnl 4,73004| 7,85320| 10,99561 | ~ 4,5 : 
®, | 12,303 46,050 | 98,906 171,5 85 263,998 


Vernachlassigt man die auBere Dampfung, so ergibt sich aus (15) fiir die Grundschwin- 
gung 


1~ oP (PF + 12,30J) 202 (PF + 12,30 
500.6 EJK | 


ae 500,6 EJ — 12,3012 P =| 500,6 EJ K ) 
(16) 


oi 202 (2F + 12,30J) * 

Dieses Ergebnis ]aBt sich anhand der zitierten Arbeit von Sezawa kontrollieren. Dort wurde die 
Differentialgleichung naherungsweise fiir den Grundton integriert, wobei sich also auf 
einem Weg, der von dem hier eingeschlagenen grundsatzlich abweicht, fiir die Eigenwerte a, 
und 6, in unserer Schreibweise die Beziehung! 


: 500 EJ ri 5 oa 2P 500 EJri 2 
Sy ane ale .7 /500 EJ — 12,74 zal ; ) 
ie 20 1%(PF + 12,74 J) 7 i|/ o(2F + 12,74J) 2 ol(PF + 12,74) 


ergab. Der Vergleich mit Formel (16) zeigt, da®B beide Beziehungen genau denselben Aufbau 


besitzen und nur in den Zahlenwerten etwas voneinander abweichen. 


4, Zusammenfassung. Die strenge Behandlung der Differentialgleichung des frei schwingenden 
Stabes, welcher unter dem EinfluB einer Langskraft steht, sté8t auf Schwierigkeiten, wenn die 
innere und auBere Dampfung beriicksichtigt wird. Bei Fehlen einer Langskraft, aber Beachtung 
der Rotationstragheit und Dampfung ist es ebenfalls schwierig, die strenge Lésung zu finden. 
Wie gezeigt wurde, kann eine Naherungslésung fiir simtliche Eigenschwingungszahlen und 
Dampfungsexponenten auch fiir den allgemeinen Fall, namlich bei Vorhandensein einer Langs- 
kraft und Beriicksichtigung der Dampfung sowie Rotationstragheit, gefunden werden, wenn 
man nicht von der Differentialgleichung sondern von den Lagrangeschen Gleichungen ausgeht, 
wobei durch Hinzunahme eines Gliedes, der Zerstreuungsfunktion, der Dampfung Rechnung 
getragen wird. In den Ausdriicken fiir die kinetische und potentielle Energie sowie die Zerstreu- 
ungsfunktion wurden die unbekannten Schwingungsformen des Stabes durch benachbarte 
Funktionen ersetzt. Die Methode kann demnach als eine Erweiterung des Rayleighschen Ver- 
fahrens auf Systeme, die unter Wirkung von Reibungskraften stehen, angesehen werden. 


(Eingegangen am 12. Mai 1951.) 
Anschrift des Verfassers: Dr.-Ing. habil. H, Eschler, Téging am Inn, Hauptstrafie 81. 


1 Vel. FuBnote 2 von S.1; ebendort S. 277, Gleichung (24). 
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Zusatz zu J. Dorr: Strenge Lésung der Integralgleichung 
fiir ein Fliigelgitter. 


Von K. Nickel. 
In seiner kiirzlich hier erschienenen Arbeit! betrachtet Herr Dérr die Integralgleichung” 
1 
if 
(2) = 5 | T(é) tg = (x—8) dé (1) 
= 


und zeigt, daB ihre allgemeine Lésung die folgende Gestalt besitzt: 


k Gin = Coj = ¢ ey be x 

ale | h 

(D(a ee : = dx. (2) 
al) Sin? + — 6ine = 5 


I | Sin? a — Sin? > é = (x <5) 


—1 


Dabei ist k eine beliebige Konstante. Bei der Anwendung von (1) auf die Strémung durch ein 

Fliigelgitter interessiert vor allem diejenige Lésung /'(€) von Ss die an der Stelle €=1 ver- 

schwindet (Abstrémbedingung). In der angefiihrten Arbeit wird dariiber nur gesagt, da man 

durch passende Wahl von k erreichen kénne, daB /’(1) = 0 wird. Dies sollim folgenden allgemein 
1 


durchgefiihrt werden, weiter werde noch eine Formel fiir den ,,Gesamtauftrieb‘‘ [T® d€ 
Si 
angegeben. 


Es gelten die Formeln 


Gin? | — Sin? «= Gin Z (1+ x) Gin Z(1— 2) (3) 
und 
Gin = —s) Goji -A—H 
1 Sit in —é SO i ——s 
; Bg = (x— 2). (4) 


sett: TayT pibeshd nos, Sas ic 
Gof | (*—&) Sin | (l—2) Sin 7, (l—*) 


Durch Einsetzen von (3) und (4) in (2) ergibt sich 
1 a 


/ cane 
| oe =e x) 


kh 
a ae a Coj - z &—2 of > (==) p20 (x) a 
Cities ~ (L—*) 
ad, —1 
ré ae St 
hein 7, (+4) Gin 1-8 (5) 
ee = : 
Sig (ee) [ein 5 Gi 
if : in ees 
+4|/= — w (x) tg 7 (x -9 pale ae 
Sin > (1 + ¢) Gin > 7 —x) 


==] 


Es laBt sich nun zeigen, daB im allgemeinen? das Integral im letzten Summanden von (5) 


fii 
& 1 endlich bleibt. Da dieser Summand damit bei = 1 verschwindet, erhalt man a ; 


us der 


Ede Darr, Ing.-Arch. 19 (1951) S. 66. 


2 Die hier mad im folgenden auftretenden uneigent 
aufzufassen. 


° Z. B. fiir alle stiickweise stetigen Funktionen w(x). 


lichen Integrale sind als Cauchysche Hauptwerte 
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Bedingung J'(1) = 0 fiir k die Bestimmungsgleichung : 


foe hs (6) 


‘ii os Sins (1 + x) 
Sin > (1 — x) 


Setzt man (6) in (5) ein, so kommt schlieBlich die gesuchte Gleichung: 


Sa = 
; Sin = (1—E) i‘ ‘ein (1 + x) 2 
P(g) = h : w(x) , stg 7 (x —&)—Bq h dx. (7) 
) 
= 


Sin= (ds) 


Fiir manche Zwecke ist es noch von Vorteil, einen geschlossenen Ausdruck fiir den ,,Gesamt- 


1 
auftrieb*’ A = {T® d& zu kennen. Bildet man A nach (7) und vertauscht die Integrations- 
—1 


reihenfolge, so 1aBt sich das innere Integral elementar berechnen. Man findet so 
1 1 A ee 
5 Gin (1 + «) 
w(x) eee ame te ax. (8) 

Sins (1 — x) 


Wendet man diese Rechnungen auf Fliigelgitter an, so bedeutet /'(&) die Zirkulationsverteilung 
iiber die Fliigeltiefe, w(x) die an der Stelle x eines Profils induzierte Abwindgeschwindigkeit. 
wahrend durch h die geometrische Gestalt des Gitters bestimmt wird (vgl. die Arbeit Dérrs). 
Wie auch in der Arbeit von Dérr bemerkt wird, ist ["(£) bei reellem w(x) nur dann reell, wenn h 
reell oder rein imaginar ist, entsprechend einem Staffelwinkel des betrachteten Fliigelgitters 
von 0 oder 90°. Analog zu den Rechnungen von Birnbaum! kann man dann in diesen beiden 
Fallen w(x) durch die Profilkontur des Fliigelgitters bestimmen. Z.B. entspricht w(x) = const 
einem Gitter aus ebenen Platten. Fiir w(x) = 1 findet man aus (7) und (8) das bekannte? 


Ergebnis ; 


> | /cnZa—4 
LEN el VE oe peer ge 
Coj = Gin > (1 + €) 


(Eingegangen am 8. Mai 1951.) 


Anschrift des Verfassers: Dr. Karl Nickel, Instituto Aerotecnico Cérdoba, Argentina. 


1 W. Birnbaum, Z. angew. Math. Mech. 3 (1923) S. 290ff. 
2 Vel. mit etwas anderen Bezeichnungen R. Grammel, Die hydrodynamischen Grundlagen des Fluges, 
S. 93—97. Braunschweig 1917. 
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Druckverteilung unter elastisch gelagerten Kreisplatten. 
Von H. Jung. 


1, Einleitung. Eine elastisch isotrope Kreisplatte liege auf einem elastisch isotropen Halb- 
raum auf. Wird die Kreisplatte durch eine rotationssymmetrische Belastung beansprucht, so 
ergibt sich eine Druckverteilung zwischen Platte und Halbraum. Die Druckverteilung hangt 
von der Art der Belastung, der Plattensteifigkeit und den Elastizitatskonstanten des Halbraumes 
ab. Zur Vereinfachung der Rechnung wurde bisher! angenommen, daf die Platte in ihrer ge- 
samten Ausdehnung auf dem Halbraum aufliegt. Neményi? hat darauf hingewiesen, dafi sich 
die Platte vom Halbraum abheben kann, wenn die Belastung, die Plattensteifigkeit und die 
Elastizitatskonstanten vorgegeben sind. 

Es soll nun gezeigt werden, wie dieses von Neményi gestellte Problem gelést werden kann. 
In Ziff. 2 wird das allgemeine Verfahren angegeben. An einem Beispiel soll dann in Ziff. 3 der 


Rechenaufwand gezeigt werden. 


2. Verfahren zur Ermitthung der Druckverteilung. Gegeben ist eine Kreisplatte vom Halb- 
messer a und der Dicke h. Die rotationssymmetrische Belastung der Platte werde mit p,(r) und 
die Druckverteilung zwischen Platte und Halbraum mit p,(r) bezeichnet. Es wird nun ange- 
nommen, daB die Platte im belasteten Zustand mit der Flache r < « auf dem Halbraum aufliege. 
Ferner macht man noch die Annahme h <a, so daB die technische Plattenbiegungslehre zur 
Anwendung kommen kann. 

Wird der elastische Halbraum z => 0 durch eine rotationssymmetrische Belastung beansprucht, 
so ist die Verschiebung der Oberflache z= 0 mit der Abkiirzung 


8 (A) = [pr (2) @ Jo(A@) do (1) 
gegeben durch? 
a hiee it 
(= Gm | &(A) Jo (Ar) de. (2) 
0 
Die Druckverteilung zwischen Platte und Halbraum ist 
Pi(r) = [é: (A) A Jo (Ar) dd. . (3) 


Die Durchbiegung der Kreisplatte mit der konstanten Biegesteifigkeit N ist durch die Diffe- 
rentialgleichung : 


AAw, = <-[p2(r) — px (r)] (4) 

mit den Randbedingungen 
r= ea, (T= 0 , (5) 
te 62 yea) (6) 


bestimmt, wobei 
A oie mat es ist und gq, di kraf' d i : me 

dr? Eas q, die Querkraft, m, das Biegemoment bedeuten. Die Lésung von (4) 
]aBt sich in der Form 


+ €@ Jo(a) + e144 Jy(Aa) + ea 2 Jg(a) + cy.4 2? Jy(Aa)|% Jo(Ar) dA 


[EXO @ J, (Ae) do — EXO @ Jy(A@) de 
0 0 : (7) 


sik Schleicher, Kreisplatten auf elastischer Unterlage. Berlin 1926. 
* P. Neményi, Z. angew. Math. Mech. 11, (1931) 5S. 450. 
° H. Jung, Ing.-Arch. 18, (1950) S. 178. 
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angeben. Das Integral (7) mu die Randbedingungen (5) und (6) erfiillen. Setzt man Creat), 


so ist die Randbedingung (5) befriedigt. Da w,(r) eine stiickweis stetige und stetig differentiier- 
bare Funktion ist, so muf 


a 


s2(2) = 72) | pole) 0 Ilo 0) do — gy(A) + ¢, aA J,(Aa) + 0, a2 Jy(Aa) + cg a2? J, (Aa) 
0 


beschrankt und fiir 0 </ < oo quadratisch integrierbar sein. Diese Forderung fiihrt auf die 
Bedingung 
lim g, (A) +0. (8) 
A+0 
Entwickelt man g, (A) in eine Potenzreihe nach J, so miissen simtliche Glieder mit 4” gleich 
Null sein. Da g,(A) eine quadratische Funktion in / ist, erhalt man aus (8) zwei Gleichungen, 
wobei die Glieder mit A~4 auf die Gleichgewichtsbedingung fiihren. Setzt man die Glieder mit 
A gleich Null, so erhalt man eine Gleichung zur Bestimmung von c,. Die Randbedingung (6) 
fiihrt auf eine Gleichung zur Bestimmung von ¢,. 
Macht man nun die Annahme, dafs die Platte reibungsfrei auf dem Halbraum aufliegt, so muB 


Weir) =—w (nr) fur er Sa, 
Pi(r) =9 raCi gay oreo" 
sein. Dies ergibt die Integralgleichung 


a 


| mo @ Jo(A@) do — g,(A) + ead Jy (Aa) 
0 (9) 
4 a2 Jg(Aa) + c, 022 tia Jo (Ar) da. 


nm — 1 fs i r 
"im | Bil4) Solr) a= x, | 


0 0 


Um diese Integralgleichung zu lésen, macht man fiir die unbekannte Funktion g,(A) den Ansatz 
1 


aS aUe era) Te A (10) 


Bestimmt man mit (10) die Druckverteilung zwischen Platte und Halbraum, so wird! mit (3) 


awe 
nae |> Jie) (r <a) 
— 2 — ye yh 
. 0 sa) 
Die Verschiebung w, der Oberflache des Halbraumes ergibt sich® mit (10) und (2) zu 
1 
= = ) l 3 
4, ae Fa. ns 1.53) ieee ip 
—T =n! 12 
(7) = . (12) 
rae let 3) a 
ee ee eee 
2 


Bei der Berechnung von w,(r) fiihrt die Bedingung (8) auf zwei Gleichungen zur Bestimmung 
von A, und c,. Wertet man das Integral (7) mit oD aus, so erhalt man mit 


[plo ) o y(t0) do= & B, K, A, In (@)) 


1 W. Magnus und F. Oberhettinger, Formeln u. Satze fiir die speziellen Funktionen der Mathematik 
und Physik, $.51. 2. Auflage Berlin 1948. 
2 Seite 50 von Fufnote 1. 
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die Durchbiegung der Platte in der Form 


Sm. caer 4 Semen +auey +2) (13) 


Setzt man (12) und (13) in (9) ein, so wird fir r Sa 


Gm >, An 1 


n+ — 
2 2 


(14) 


F 1 gil | == > Ka) aD L,,(r?) + ¢, U(r?) + ¢3 


Z 


n! < 


Verlangt man, da® (14) auf den Kreisen r= a/u mit w= 1,2,3... erfiillt ist, so erhalt man 
n-+Illineare Gleichungen zur Bestimmung der 4, und ¢3. Die Konstante c, wird aus der Rand- 
bedingung (6) ermittelt. . 

Gabe ann? ist nun noch der Druckkreishalbmesser «. Die durch (12) gegebenen Verschie- 
bungen bleiben auf dem Kreis r—o endlich. Es gilt, wie einfach nachgewiesen werden kann, 


Abb. 1. — Deformation des Halbraums fiir « = a. Abb. 2. Platte mit Einzellast. 
—+— Deformation des Halbraums fiir « < a. 


1 1 
n— n— — 
es Dial a ees 
ae Bese n1,1) a ; pono ask 
ae Se \ ee P| =) . 
2 Sy! a 2 
Bestimmt man dw,/dr auf dem Kreis r= a, so erhalt man 
at == 500 « 
dr |= 


Eine senkrechte Tangente kann sich aber dann und nur dann ausbilden, wenn a= a ist 
(Abb. 1). Die Druckverteilung p,(r) wird dann an der Stelle r= a unendlich. SchlieBt man 
diesen Fall aus, so mu, wie einfach gezeigt werden kann, 


A,(a) = 0 (15) 
gefordert werden. Da aus der Gleichgewichtsbedingung A, als Funktion der A, und a bestimmt 
wurde, ist durch (15) eine Bedingung zur Berechnung von « gegeben. 

Ergibt sich aus (15) « > a, so ist a= a zu setzen. Die Platte driickt sich dann in den Halb- 


raum ein und A, ist durch die Gleichgewichtsbedingung bestimmt. An einem Beispiel soll der 
Rechnungsgang noch gezeigt werden. 
3. Platte mit Einzellast im Koordinatenursprung, Durch Abb. 2 sei der Belastungsfall ge- 
geben. Bei Verwendung eines dreigliedrigen Naherungsansatzes wird 
1 3 5 
£1 (A) = Ag A * Say (A) + Ap A * Jgyg (0A) + Apa * Js,o (aA). (16) 
Mit (16) geht man in (2) ein und erhalt 


/ 


m— | 7 1 y/an ere a y/an Te 3 r4\' 
0) == V+ 4,4) —t4) 445/24 45) 


ea) 


208 2 Oe 8 a4 
Geht man mit (16) in (7) ein, dann ergibt sich die Durchbiegung der Platte zu 


He, AS; A) te, Ae ee ae (| rat locatnidas | 
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In (18) wurde zur Abkiirzung gesetzt 
ea ER: 
By= Sas ee | ag’ Jy(A 0) d 
Mit der Bedingung (8) erhalt man aus (18) die beiden Gleichungen 


a pe a ) Ez: é We 2% 4.= 0, | 


2 0 x? 2aat 4 2a a4 c | 
S | seal — 
yes Ag+ 22% 4 39 “17 i+/% B10 a aaa: 
Bestimmt man aus (19) 4, und c, und geht damit in (18) ein, so 14Bt sich das Integral aus- 
werten und ergibt 


(19) 


w,(r) = By Ky + A, L,+ 4, 1,+4+ 41+ ¢;, (20) 
wobei zur Abkiirzung gesetzt wurde: Pir T= 0 g 
5 r a? 1 1 
—— 2 2 | wad ise 2 
TN ele Peel ne 7b erie nigh oneal a 
(a? — r?)3 Ol Sp)4 (a? — r?)é ors 
r 3l-dlat  o8-3las + 9% 6las ! ; 
pelae— 1 ot 1 
Dope |) | 2) | 2)2 
1 22 ei 2 Gh 8 ne asa ee 
1 il 
| 2 2)3 _| 2 Pay alee 
BER fie a dUaaia em . 
(2B 2 WAY. a a? 1 
L= : | 2 2) | 2 2)2 
2 poe Se Pin! pf) aeemarie cn 
1 i 
ee 2 2)3 2 ENE eae 
Dionink ame ey aie =) | 
1—r? 
Re 
2 4 
dagegen fiir r >a: 
K ay 
. c= — Zine, 
i 2aa% 
1 a! 
24 2 oh 
L, Visa ™* 
1—r 
L, i 
Die Randbedingung m,= 0 fiir r— 1 fiihrt auf 
B 2 a? m — | [20203 mead pe 4 a4 eel 
Hae Me aiewr ae | ore rere Viz eee 


Fordert man nun, da® auf den Kreisen r= a, r= 4/2 und im Ursprung w,(r) = w,(r) ist, 
so erhalt man mit der Abkiirzung 


die drei Gleichungen 
By + fi() + Ai 81(%) -F A, h, (x) — ¢,= 0 
By + fa () + Ay 8 (0) + Aa hy (4) — ¢3= 9, 
By f(a) + Ay g3(«) + Ap hg (a) — ¢3 = 0 


(21) 
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Eliminiert man aus (21) ¢3, so wird in ausfiithrlicher Schreibweise 
— nfo te +east] a ES +) Eta SSallt 
ee Sosy 4 Soro ae a et ye2—i|=0, (22) 
aalieret fom? + Toman + Alas) +5) + tm) ees + 


7B OT 5704 1/20 ae ne m — | 
| 5 Be CY as ‘(== 0, 
alee wey 90-21) a m+1 
Streicht man in (22) die Glieder mit «? und héherer Ordnung, so erhalt man die Naherungs- 


gleichungen 


By/an 5 OT 109 
a Bt Aaah Betas: | 


5 By /an OT 3 LOL 

A, 7 2 | Ay 7554 «pyle 1 "398 

Aus (22a) bestimmt man sich erste Naherungswerte fiir A, und A,. Mit diesen A, und A, 

wird in die erste Gleichung aus (19) eingegangen und A,(a) bestimmt. Mit A,(«) = 0 erhalt 

man einen ersten Naherungswert a,. Mit diesem a, geht man in die Glieder héherer Ordnung 

von (22) ein und erhalt damit verbesserte Gleichungen (22a), aus denen dann auf dem beschrie- 

benen Weg ein neues a, ermittelt wird. Das Verfahren wird solange durchgefihrt, bis sich a, 
nicht mehr andert. Bei der Rechnung ist es nun zweckmafig, die Plattendicke 


(22a) 


h=x0 


zu setzen, wobei der erste x-Wert geschatzt werden mu. Bei dem vorgegebenen Beispiel ergab 


sich mit h= = aund a= 1 der Druckkreisradius zu 
1 
es 366 < 


(Eingegangen am 22, Mai 1951.) 
Anschrift des Verfassers: Dr. H. Jung, Stuttgart-N., Robert-Mayer-Str. 66. 
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Durchgang einer Potentialstérung durch einen Leitschaufelkranz. 
Von H. Séhngen. 


1, Aufgabenstellung. Bei Schwingungsuntersuchungen an laufenden Schaufelkrinzen axial 
durchstrémter Maschinen tritt éfters die Frage auf, wie sich eine Stérung, die z. B. von einer 
Stiitzrippenreihe erzeugt wird, auf einen dahinter bzw. davor befindlichen Laufschaufelkranz 
auswirkt, wenn zwischen beiden noch ein Leitschaufelkranz liegt. Um diese Frage beantworten 
zu kénnen, muf man zunachst einmal wissen, in welchem MaBe die Stérung durch den Leitkranz 
hindurchdringt. In der Praxis interessieren insbesondere zwei Falle: 

I. Stiitzrippenreihe vor dem Eintrittsleitrad (Abb. 1), d. h. die Stiitzrippenreihe liegt in einem 
Parallelstrom vor einem Beschleunigungsgitter, wobei Stiitzrippenreihe und Gitterachse zur an- 
kommenden Strémung senkrecht stehen. Unter den unten genannten Voraussetzungen wird 
gezeigt, daB in diesem Falle die Stérung nicht durch den Leitapparat hindurchgeht. 

II. Stiitzrippenreihe hinter dem Austrittsleitrad (Abb. 2), d.h. die Stiitzrippenreihe liegt 
hinter einem Verzégerungsgitter und es stehen Stiitzrippenreihe und Gitterachse senkrecht zum 
Abstrom. In diesem Fall wirkt sich das Leitgitter ungiinstig aus, es verdoppelt die Stérung und 
holt sie um die Tiefe des Leitgitters naher heran. 

Beide Falle werden unter den Voraussetzungen untersucht, da 

1) das strémende Medium inkompressibel und reibungsfrei ist, 

2) es sich um eine Potentialstérung (Verdrangungsstrémung) handelt, 

3) die Teilung a des Leitgitters klein im Vergleich zur Schaufeltiefe c ist, d. h., die Beschran- 
kung auf eine lineare Theorie in a/c. Dies ist gleichbedeutend mit der Annahme, daB in dem 
gesamten Schaufelraum die Strémungsrichtung den Schaufeln parallel ist, 

4) die Stiitzrippenzahl klein im Vergleich zur Blattzahl ist, 

5) die Héhe des Ringkanals, in dem der Kranz liegt, klein gegeniiber dem Durchmesser ist, 
man den Strémungsvorgang also in eine Ebene abwickeln kann, 

6) die Schaufeln des Gitters einander gleich und diinn sind, 

7) die Stiitzrippenreihe soweit von dem Gitter entfernt ist, daB ihre Anstrémung als homo- 
gen angesehen werden darf und mithin die Riickinduktion vernachlassigt werden kann, 

8) die Stiitzrippen gleiche und symmetrische Profile habent. 


2. Geschwindigkeitsfelder und Randbedingungen, Wir denken uns den Strémungsvorgang im 
Ringkanal in eine (x, y)-Ebene abgewickelt, und zwar mége dabei die Gittervorderkante mit 
der y-Achse zusammenfallen. Der Nullpunkt werde so gewahlt, daB die x-Achse Symmetrie- 
achse fiir eine Stiitzrippe ist. 

Hat der Schaufelkranz den Halbmesser R, und ist o die Anzahl der Stiitzrippen, so hat die 
Strémung in y die Periode 2 R/o. Es ist daher zweckmaBig, Fouriersche Reihen einzufihren. 
Wir benutzen diese hier in der komplexen Form und schreiben z. B. fiir die Str6mungsgeschwin- 


digkeiten 


i 

| 

i Me 

a 
ry 
3 
a 
eae 
® 
e 
—_— 
8 
= 
! 

| 

i Me 
3 
PA 
= 
d 


u(x, y) = 


wobei wir unter u stets eine Geschwindigkeit in der x-Richtung und unter v eine solche in der 


y-Richtung verstehen. 
Handelt es sich speziell um eine Potentialstrémung, so mu fiir n + 0 


o o 
(aye de Ba Blew BL, | as 


oO 
v,(x) =iA,e 8 —i Be R 


1 Diese Voraussetzung ist nicht wesentlich, Die beiden genannten Hauptergebnisse gelten auch dann, 
wenn die Stiitzrippen verschieden, unsymmetrisch bzw. angestellt sind, Die Beweisfiihrung ist dann ein 


wenig abzudndern. 


14 Séhngen: Durchgang einer Potentialstérung durch einen Leitschaufelkranz. Ingenieur-Archiv 


und fiir n= 0 
U(x) = Ag, —% (x) = By (2) 
sein, wobei die Ax und Bx Konstante sind. 
Betrachtet man insbesondere die Stérung, die im Falle I von den Stiitzrippen am Orte des 
Leitapparates erzeugt wird, wenn dieser nicht vorhanden ist, so kann diese, da die axiale 
Stérgeschwindigkeit in y gerade ist und mit wachsendem x abklingt, wegen (1) in der Form 


ul, 40, = Wa >) a, € eae (3) 
n=1 
geschrieben werden, wobei die a, wegen 8) reell sind. Die GréBe a, Wo gibt die Amplitude der 
no Harmonischen der Stérung auf dem Umfang an. Die a, sehen wir hier als gegeben an, sie 
sind durch die Stiitzrippenreihe und ihre Lage zum Gitter bestimmt. 
Ganz analog kann in dem Falle II die Stérung, welche die mit der Geschwindigkeit W, an- 
gestrémte Stiitzrippenreihe am Orte des Gitters erzeugt, falls dieses nicht vorhanden ist, in der 
Form 


fesse es W, > b, see (4) 
n=1 
mit reellen 6, geschrieben werden. Die 6, seien gegeben. 

In dem Leitgitter muS die Strémung an den Schaufeln diesen parallel gerichtet sein. Um 
dies zu erreichen, wird man die Schaufeln mit Wirbeln belegen. Da aber die Schaufelzahl als 
unendlich gro vorausgesetzt wurde, so bedeutet dies eine Belegung des gesamten Schaufel- 
raumes mit Wirbeln. Es ist zweckmafig, dabei diejenigen Wirbel, die auf einer Parallelen zur 
Gitterachse liegen, zu einer Schicht zusammenzufassen. 

Wir bestimmen zunachst das Feld einer derartigen Wirbelschicht, die mit der y-Achse zu- 
sammenfallen mége. Die Wirbeldichte langs der Schicht sei /'(y). Wir schreiben 


und berechnen das zu I', ound gehorige Feld, wobei wir noch n + 0 annehmen wollen. Wegen 
(1) ergibt sich fiir die von dieser Schicht induzierten Geschwindigkeiten, falls 


x <0 0< x | 
wee ae | 
Une elne ; Teenie ute 
[n| Z x in| = 
My = 1 Cue Sa t= — tT Dye ES 


Wegen der Kontinuitat beim Durchgang durch die Schicht mu8 C, = D,, weiter muB 


%»(—0) — 0, (4-0) = 237 Dy = ‘Ae 


|n 


sein. Damit ergibt sich fiir das von J', induzierte Feld 


<a ee 1 ies n ne % : ts 2 
TL isa Ts (5 ok (nm == 0). (5) 
Ist n= 0, so findet man, da eine derartige Schicht nur stromabwarts induziert, 

Ug — 1 V9 = aaa 


Insgesamt erzeugt also die bei x= 0 liegende Wirbelschicht das Feld 
BE Oe & See a) 
u eee Gaake : 


Beachtet man noch, daB J‘(y) reell und mithin I_,= J", ist — durch den Querstrich wird 


a n 
der Ubergang zur konjugiert komplexen Gré®e angedeutet — so kann man auch schreiben 


6G 
Sil 


Stes eit | ye SS ee ee 
in “Al | 1) Q al may Ine zinc i Shia ping R 


x+ty) 
. 


(6) 
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Mit derartigen Schichten hat man den gesamten Schaufelraum zu iiberdecken. Au®erdem 
muf} man noch, um die Randbedingungen befriedigen zu kénnen, eine isolierte Wirbelschicht 
langs der Vorderkante anbringen. Diese entspricht der Singularitét an der Vorderkante der 
Platte. Bezeichnet man die Wirbeldichte der Schicht lings der y-Achse wie oben mit J” (y) 
und die Flachendichte im Inneren des Gitters mit y(x, y), so erhalt man fiir das von der Wirbel- 
belegung des Schaufelkranzes induzierte Feld nach (6) die Darstellung 


ae Shon det ae i #\~ no (xtiy) 
le zPo(l +) x (! alge ag 


7 no (xtiy) i =p; 
at | Dee Fg | volte ee ae 
- (7) 
ixer x —E\ mo (x—é+iy) 
espinal @(1— ea) x dé + 
rea () 
Cy *—é Swe (e=€44y) 
Oo 
n=1 9 i 


Dabei ist die Tiefe des Leitkranzes gleich c gesetzt. 
Die Wirbelbelegung des Schaufelraumes ist so zu bestimmen, daf die resultierende Strémung 


U—iV=(u,+ u,+ u,,) —i (vt v.+ v,), (8) 
wobei u,— iv, das Anstrémfeld darstellt, also 
Ug —it,z= W, im Falle I 


und 

Ug — tVg= W,e-*’o im Falle II, 
zu den Schaufeln tangential verlauft. Ist die Skelettlinie der Schaufeln durch 

Vat (te eturee 0) ex is € 
gegeben und setzt man 
f(x) = tg B(x), 

bezeichnet also mit {(x) den értlichen Umlenkungswinkel, so muf8 im Schaufelraum, d. h. fiir 
Nee <Gics 

— Usin B(x)+ Vcos (x)= 0, 
also 

U—iV=(U+i V)e—?P (9) 


sein. 
Da in den beiden hier zu betrachtenden Fallen das Stérfeld in der Form 


Timi oye te (n+0) (10) 


geschrieben werden kann, so geht die Randbedingung (9) mit (8), (7) und (10) in die folgende 
Gleichung iiber: 


o bp 
<n & (x+iy) ‘ 


(x+ty) ; 7 ‘ 
U,z—tv,+ >» A,e ® Sid Bye at Ra Ei ome -++ i f y9(&) dé — 
n=—co n=1 0 
of n= (w—E+iy) ~ x 2 (x-§+iy) 
=) | (6) 2 ® de +i Ss fle * d= 
n=1x n=1 
, oo Page eS oo Soe (x—iy) 
ad ee Meer SPD ee ie ae 
n=—oco n=l 
‘ a ARE Te na (x—§—iy) cect #3 a an & (x— Py | 
—ify(écé+iy f ynl€)e dé—i S | mn(6) g| 
0 ee n=1 (0 
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vil Mivgt coke 
Betrachtet man hierin die Glieder mit dem Faktor e 8% 5 (n + 0), so miissen diese fiir sich 
genommen die Gleichung befriedigen, d.h. es muB in 0 <x <e 


ve Re if yim (é) OO Ra de pn aE) A PS Plage 
(11) 

gs x x% ST °F ops 

Ri fy Ge ay | 

0 


|n| 


—il,e 


Se i < 
gelten. Die Glieder mit dem Faktor-e — 'R” Jiefern die zu (11) konjugiert komplexe Be- 
ziehung, also keine neue Bedingung. Fiir n= 0 ergibt sich noch 


Hea Dey BSL on eis aaree hpe kin, = 410, i [ 70 (6) a8]. (12) 
0 


Die Gleichungen (11) und (12) stellen die Integralgleichungen zur Bestimmung der Wirbel- — 
verteilung im Schaufelraum dar. Die Gleichung (12) entspricht dabei dem ungestérten Fall, 
aus ihr folgt 


ug tg B(x) — va + Ty + f yo(é) dé = 0 
und daraus fiir x — 0 : 
Pi = 0, — iu," f’(0) ed 
bzw. durch Differentiation 
Yo (x) = — ua f’"(*). (13h 


Damit ist die Gleichung (12) gelést. Die Gleichung (11), d.h. den Stéranteil, behandeln wir 
fiir die beiden Falle I und II getrennt. 


3. Stiitzrippenreihe vor dem Eintrittsleitrad. In diesem Falle ist die Anstrémung durch 
U,= W, .und .v2=—0 
gegeben. Weiter ist nach (3) und (10) 


Wo ye A, == 0, ~wenn vit 0 | 


und 


A-— Was, Wenn prt <a 0. 


Geht man damit in die Integralgleichung (11) ein, so geht 


ies 
ae diese in die Gleichung 
‘Lee 


MN, e = é 2n 7 x—2iB(x) 
—n= —2n = x—2tB (x 
—fr(Gje "Ri dgae""® 
—— C—_ Ed 
Abb. 1. Stiitzrippenreihe vor ae no € 
Beschleunigungsgitter. W, Fa a 1 PS = if Vn (€) e R dé 
0 


nibie y= lA oon mil ) Kae Ke tilt, Dee Gleichung wird offenbar durch 
[,=—iWyan,,  yn(é)=0 (14) 
gelost. Damit ist die Wirbelverteilung fiir das Gitter vollkommen bestimmt und das resultie- 


rende Geschwindigkeitsfeld U—iV kann in der Umgebung des Gitters berechnet werden. 
Hierfiir erhalt man aus (8), (3) und (7) unter Beachtung von (13a, b) und (14) den Ausdruck 


0 1V=M—EM(I+ EH) FO $M f $f. + =F] ae + 
(15) 


1 x ~ ok ot ng tin) 
+2 M(l— 7) > = —eR Ces 


x| n=) 


Daraus ersieht man, daf hinter dem Gitter die ungestirte Abstrémgeschwindigkeit herrscht. 
Die Stérung dringt also in diesem Falle nicht durch das Gitter hindurch. 


XX. Band 1952. Séhngen: Durchgang einer Potentialstérung durch einen Leitschaufelkranz. a 


4, Stiitzrippenreihe hinter dem Austrittsleitrad. In diesem Falle ist das Anstrémfeld durch 
Uz— tv,= W, e *Po 
gegeben, und fiir das Stérfeld gilt 
B= oun tO, 
AL 0 fir n <0, 
so daB also die Integralgleichung (11) hier die Form 


Wb, —i f ya(B)e We weeds Vt fr ag (16) 


fursiti— leno und 0 <4 < clannimmet. 

Im Gegensatz zu Ziff. 3 wollen wir uns hier auf die Betrachtung einer speziellen Form der 
Skelettlinie beschranken und annehmen, dafs B(x) eine lineare Funktion von x ist. Weiter 
wollen wir voraussetzen, daS im ungestérten Fall der Eintritt in das Gitter stoBfrei erfolgt 
(I) = 0), also B(0) = By ist. Da weiter die Abstrémrichtung durch B(c) gegeben ist und hori- 


zontal sein soll, so ergibt sich 


B(x) = Bo: [1 — = ae — 
und mithin | } va ‘ 


f(y = g, meee By(l——). = ia: 


Um die Integralgleichung (16) zulésen, differenzieren 
wir sie nach x und erhalten durch eine leichte Um- 


formung By Vin 
Wy, 


= (» ue == 4 2) [r, +. j Vn (€)e tmp Hae] (17) Abb. 2. oor el upenrene hinter einem 


erzogerungsgitter, 


und daraus durch Differentiation und Division durch 1 + od Nd 


yala) + (252+ in) yale) tg B(x) = 


die Differentialgleichung 


also 
ese eh Sone are ed 
— |2— +in —| f(x) 
V(x Ce | ae n)/ : (18) 
wobei die Konstante C—y,(c) noch zu bestimmen ist. Dazu multipliziere man (16) mit 


oO 
2 = L x 
e (a cms ) und differenziere. Man erhalt dann 


nc dre x +21 B(x) , 27 f(x) 
2 |r ; | Wb, — if yo(6) oT ae] om SS rgey ee 
und daraus fiir x= c 
o B 
C= ule) =i (ne — i) Mi ba Lo (19) 
Damit ist y,(x) bestimmt. Weiter folgt aus (17) fiir x= 0 
2(nz iE) P= 2+ (0) #008 fi 
also 
, ifotn oe (24 rink) fO) _ ALA bn gibotm 5 (et -ih) i 
=i Wb, cos py ¢ aaa rer ae : (20) 
dabei ist (Abb. 2) f(0) = — h gesetzt. 
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Mit den Formeln (18), (19) und (20) ist die Wirbelverteilung im Schaufelraum bekannt und 
es kann das von den Wirbeln induzierte Feld berechnet werden. Wir beschranken uns hier auf 
die Bestimmung dieses Feldes vor dem Gitter und erhalten hierfiir nach (7) 


fo 2] 


nti ees res iS | née 


n=1 n= 


ne (xtiy) no (x—€+iy) 


dé. 


= 


Da aus (16) fir += 0 


e Ses 3 —2% By 


—ify(ée ® dé=—W,b,—il,¢ 


folgt, so kann die von dem Gitter induzierte Geschwindigkeit auch in der Form 


o . 
= (x+Tty) 


fo) n ( 
— W,: Sbe ® 
(ies 


ne (x+ty) 


—2ifp, 
Up — UVUyp — 


—i (1 +e by Ee 
n=1 
dargestellt werden. Setzt man darin den Wert (20) fiir [7, ein und beachtet (4), so folgt 


o 

" . oo n— (xtiy+e+ih) 

uy — ivy = — (u,—iv,) +2 W, Sd) b,e ® ; 
n=l 


Fiir das gesamte Geschwindigkeitsfeld vor dem Gitter ergibt sich also 


—iB, xtiyte+ih) 


Ee ER Sie 
I 

Das Stérfeld wird also in diesem Fall durch das Gitter nicht ausgeléscht, sondern verstarkt. 
Sieht man von dem Faktor 2 zundchst einmal ab, so verhalt sich die Stérung so, als ob die an 
der Hinterkante herrschende StérgréBe langs des Gitters bis zur Vorderkante gleitet und sich 
von dort aus nach dem Potentialgesetz ausbreitet, d.h., der Abstand des betrachteten Auf- 
punktes bis zur Stérquelle wird gewissermaBen um die Gittertiefe verkiirzt. AuBerdem tritt 
noch eine Verdoppelung dieses Feldes ein. Diese Aussage ist hier allerdings nur fiir den Fall be- 
wiesen worden, daB der értliche Umlenkungswinkel f(x) eine lineare Funktion von x ist. Ver- 
mutlich gilt sie aber auch noch fiir allgemeinere Verteilungsformen der Umlenkung. Weiter war 
noch vorausgesetzt, da die Schaufeln des Gitters alle einander gleich sind. Es laBt sich zeigen, 
da bei gegebener Stérung die Schaufeln so verformt werden kénnen, daB die Stérung durch 
das Gitter nicht hindurchdringt. 


5. SchluBbemerkung. Die wesentlichen Ergebnisse dieser Untersuchung, da eine Stérung 
vor dem Gitter durch das Gitter vollstandig ausgeléscht, eine Stérung hinter dem Gitter durch 
dieses nach vorn verschoben und verdoppelt wird, sind hauptsachlich durch die Voraussetzung 
bedingt, daB die Gitterteilung (a) klein gegeniiber der Tiefe (c) des Schaufelraumes ist, man sich 
also auf eine lineare Theorie in a/c beschranken kann. Dies ist gleichbedeutend damit, daB in 
dem gesamten Schaufelraum — insbesondere also auch am Austritt — die Strémungsrichtung 
vollstandig mit der Richtung der Gitterschaufeln iibereinstimmt. Um zu einer Abschatzung des 
Giiltigkeitsbereiches zu gelangen, kann man z. B. fiir den Fall eines homogen angestrémten 
Streckengitters die Ergebnisse, die die lineare Theorie fiir den Auftriebsbeiwert liefert, mit den- 
jenigen der exakten Theorie vergleichen!. Es zeigt sich, daB man bis zu Werten von a/t = 0,8 
(t= Schaufeltiefe) eine sehr gute Ubereinstimmung erhalt. Ist die Anzahl der Stiitzrippen nicht 
zu groB, die Anstrémung des Leitgitters also einigermaBen homogen, so wird man annehmen 
diirfen, dafs die lineare Theorie bis zu dem eben genannten Wert auch im vorliegenden Fall noch 
gut brauchbare Ergebnisse liefert. Ist die Teilung gréfSer, aber noch nicht zu groB gegeniiber der 
Fliigeltiefe, so hangt die Strémungsrichtung am Austritt nur wenig von der Richtung der Zu- 
strémung ab, so daB die gefundenen Ergebnisse wenigstens noch angenahert zutreffen werden, 


(Eingegangen am 24, Mai 1951.) 


Anschrift des Verfassers: Dozent Dr. Heinz Séhngen, Darmstadt, Technische Hochschule, Inst. fiir Mathematik. 


1 F, Weinig, Aerodynamik der Luftschraube, S. 139, Fig. 140, Berlin 1940. 


| 
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Vergleich zweier Verfahren zur Ermittlung der spezifischen 
Oberflaiche des Mahlgutes. 


Von G. Matz, 


1, Einleitung. Fiir alle Mahl- und Zerkleinerungsvorgange der Technik und auch fiir andere 
Kornverteilungen spielt die bezogene Oberflache, d.h. die Oberflachenentwicklung je Gewichts- 
einheit eine groBe Rolle, da sie das eigentliche Ma8 des Feinheitsgrades der Mahlung ist. Die rech- 
nerische Ermittlung der spezifischen Oberflache (m2/kg) des Mahlgutes ist kein an sich neues 
Problem, aber die Unstimmigkeiten in alteren Arbeiten und die Tatsache, daB es zwei Verfahren 
gibt, die z.T. recht verschiedene Oberflachenwerte liefern, machten eine exakte Durcharbeitung. 
des ganzen Gebietes notwendig?. 


2, Das Oberflichenintegral nach Rammler. Nimmt man mit Rosin, Rammler und Sperling an, 
da die Kornverteilung des Mahlgutes ein Exponentialgesetz befolgt, das die Sieb-Riickstande 
mit der Korngréfe bzw. Siebmaschenweite durch die Beziehung 


Re 100 en 4" (1) 
verknipft, so erhalt man nach Rammiler? fiir die spezifische Oberflache den Wert 
Imax 
On UI tnt | dn—2 gba" Jd), (2) 


Imin 
Hierbei ist R der Sieb-Riickstand in Gew.-%, d die KorngréBe bzw. Siebmaschenweite in ju, d’ 
die durchschnittliche (vom Material abhangige) KorngréBe, bei der der Riickstand auf et= 36,8% 
gesunken ist, n die GleichmaBigkeitszahl (je kleiner n, desto héher der Feinanteil des Gekérns; 
n= © bedeutet gleichkérniges, isodisperses Gut), h—=1/d’", O die spezifische Oberflache in m*/kg, 
s das spezifische Gewicht des Mahl]gutes in kg/m*, f ein vom Matcrial abhangiger Formfaktor 
(Heywood-Faktor). 

Wir wollen die spezifische Oberflache O fir n=0,6; 0,7; 0,8; 0,9; 1,0; 1,1; 1,2; 1,3; 1,5; 1,75; 
2,0; 2,5; 3,0 und oo berechnen. Der von Rammler? gefundene Ausdruck fiir das Integral in (2) 


= 


(0 106 b q’—} bd? n—l1 b? a n—l bp d’ n—l1 - ee 


s Ugo ly (ye Gn = 2! Ga hele d 


1 
IX ell 
Zur Auswertung entwickeln wir zunachst die Exponentialfunktion des Integranden in eine 
Reihe und erhalten 


befriedigt nicht, da er fiir n = 1, Lx emoug ews [allgemein i )wnendlich wird. 


k—0 
Durch Einsetzen in (2) ergibt sich 
ae A jar Sanh aia 
O=-——_fnb Dea Ie 
dmin ae 
Gnax 
oo i qh” 
Ne a) | dn—2 —- d(d), 
k=0 d 


min 


1 §, KieBkalt und G.Matz, Z. VDI, 93 (1951) 5S. 58. 
2 E. Rammler, Verfahrenstechnik ZVDI-Beiheft 1940 Nr. 5. 
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d 


max 


oo k 
= > 2 i dn(k+1)—2 d(d) 
k! 


min 


: 6- 10° 
mt Oe pea ND 
Jetzt miissen zwei Falle unterschieden werden: 


1) n= eq mite == 0 2a ere Danunee 


Ber ee (5% 1 Pree Dailtae (Gans | | 
aaa n(k+1)—1 He tacit K! aye 


Hierbei ist durch S”’ symbolisiert, daB das Glied k= K = a 1 in der Summe nicht auftritt. 


max 


is (3) 


muy 


d 
d 


1 F : 
2) n RE miteko—— 02D anmesh 
y 1 1 bey =a Imax 
oe C2 ae k! ,(&+1)—1 id Sy J Bree (4) 


Nur im Falle (2) kommt also kein logarithmisches Glied vor, das Rammler ganz allgemein weg- 
gelassen hat. Wir miissen uns indessen fiir den‘ haufigen Fall 1) der Formel (3) bedienen. 


3. Vereinfachung und Neuauswertung des Oberflaichenintegrals durch Substitution, Die Glei- 
chungen (3) und (4) sind aber noch zu schwerfallig, da sie den dimensionsbehafteten Faktor b 
enthalten. Wir lassen b= (d’)—” in (2) stehen, substituieren dort y = (d/d’)" und erhalten 

d 


Max 


d\n 
dem (7) d"— d(d), 


Ymin | 
Gleichung (5) kann abnlich wie (2) durch Reihenentwicklung des Integranden ausgewertet werden. | 
Auch hier miissen analog zu oben die beiden Falle 1) und 2) unterschieden werden: | 


1) n= EG mit K=0,1,2,.... Danniet 
Ymax Oh a _ (—1)" Ymax a & aa 1)! j - 12 rane | 
| ECE = a i Re eae i » fn 
k=0 k=0 k+1—— 
(— 1)< Ymax 1 i Ymax 4, (— 1\e | y* Ymax (— 10y 9 Ymax 
pS ln ly | ae 
K! Ymin Ymin a jon yes ls Ih Senin K! Ymin 
m™ 
2) n + ep mit K=0, To as RACES 
Ymax l Lys nas k 
a ae —y ai hs =|] k Ymax 
| of e dy= Ly J Ymin 1 Eales. . (7) 
Ymin wer Opeiae oy 


Nach (6) wurde die Oberflache fiir n = 1 berechnet, nach (7) fiir alle anderen n. 


4, Die Festlegung der Integralgrenzen, Von grundlegender Wichtigkeit fiir die Auswertung 
des Integrals erweist sich die Festlegung der Grenzen. Bei der oberen Grenze ist das recht einfach. 
Rammer schreibt: ,,Mit bester Annaherung kann man die obere GrenzkorngréBe setzen zu 


a 


Tee = a ; (8) 
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Daraus ergibt sich sofort 
Van = 1 L000 == 6,92 2 7, (9) 


Wir haben iiberall mit der oberen Grenze 7 gerechnet; der Fehler, der dadurch entsteht, da® man 
nicht d,.g,—= ©© setzt, ist duBerst gering. Der Wert nq, ~ 7 entspricht einem Grenzriickstand 
ees, 8 

Wahlt man (physikalisch nicht realisierbar) fiir die untere GrenzkorngréBe d,,;,,—= 0, so di- 
vergiert das Integral in Gleichung (2) fiir n< 1. Dies ist nicht etwa ein Widerspruch 
zwischen Rechnung und Erfahrung. Im Gegenteil, je feinverteilter das Korn 
ist (n< 1) und je kleiner die einzelnen Kérner werden (d+0), um so riesenhafter 
mu die Oberflache anwachsen. Da aber das Exponentialgesetz der Kornverteilung, das 
dem Integralausdruck zugrunde liegt, nicht fiir beliebig kleine d gilt, mu’ eine endliche untere 
Grenze gréBer als 0 festgelegt werden. Hier gibt es zwei Methoden: 

Nach Methode J setzt man die kleinste Korngréfe so ein, daB hierfiir R— 99,9°% wird. Nach 


Rammler errechnet sich dann die kleinste KorngréBe zu 


= 
‘0,0004 
d * — eed aaa: ° 10 
min b leg e ( ) 
Und hieraus folgt sofort in unserer einfachen dimensionslosen Schreibweise 


toes 0,0004 9,22 - 10-4. 
lg e 


Wir haben iiberall diese untere Grenze unseren Berechnungen zugrunde gelegt, weil sie sich der 
von uns benutzten dimensionslosen Schreibweise anpaBt und die untere Grenze yj, ZU einer reinen 
Zahl wird. 

Bei der Methode II wird d,,;, = 0,1 4 angenommen. Dadurch wird yj, = (0,1/d')” zu einer 
Funktion von d’, also wird auch das bestimmte Integral keine reine Zahl wie nach Methode J, son- 
dern eine Funktion von d’. Natiirlich unterscheiden sich die Ergebnisse beider Be- 
rechnungsmethoden und dirfen nicht miteinander verwechselt werden. 


5. Die endgiiltigen Oberflachen-Formeln, Da bei beiden Methoden y,,;, <1 ist, geniigt es, in 
den Reihenentwicklungen (6) oder (7) fiir die untere Grenze das gréBte Glied mit k= 0 heran- 
zuziehen. Dann folgt 


far n +1: GréBtes Glied =~—— 
er 
sure —— 1: GréBtes Glied = Iny. 
Also lauten die endgiltigen Oberflachenformeln 
1 1 
- 106 pe Caer tay 
gee f Vieni OY 
Sl 
9,22 10-4 
Fir n= 1 kommt 
See eee Ue 9,20 10-4 ul 
O= ; Bopte | Ra El n 9, P (11) 


fir n = 0,5; 0,7; 0,8; 0,9; 1,1; 1,2; 1,3; 1,5; 1,75; 2,0; 2,5; 3,0 u. co: 
6-106 f | ey (lens eee i as 
CS 7 re anche 1 Ps oD) 
0k ol | —— 


nm 


Ferner gilt mit O= 0,, fiir f= 1 und mit s= 10° kg/m? 
Ta 
a ERG 
Nach den Formeln (11) und (12) sind die Werte der beigefiigten Tabelle 1 berechnet. Bei der 
Summenberechnung ergeben sich folgende Schwierigkeiten: Da die Summen fiir groBe n (2,5; 


1 ET Pne F 
On . : 3 == 6000 = c . td in ; Oj, in m?/kg). 
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3,0) klein sind, miissen mindestens 19 Reihenglieder (k= 0, ... 18) berechnet werden, um eine 
Genauigkeit von + 0,05im Endergebnis zu erzielen. Eine wesentliche Vereinfachung erreicht man 
durch Einfiihrung des Eulerschen Integrals, das aber nur fiir n > 1 verwendet werden darf, wah- 
rend fiir n< 1 die Reihenentwicklung beibehalten werden muB. 

Die Tabelle 1 zeigt namlich (man vergleiche die Spalten 2 und 4), daB fiir n = 1,1; 1,2 und 1,3 


das Eulersche Integral 
one! 
Ey Dis 1 
jy tee=ne— ay my 
6 


nur wenig gréBer ist als die Summe der Spalte 2, d. h. das Integral 


Tabelle 1 (nach Methode I fiir Rmax = 99,9%). 
Spalte 1 Spalte 2 Spalte 3 Spalte 4 Spalte 5 Spalte 6 


1 
a(—5) 
es On,° d’- 10-3 
(9322. 05"), =6: Sp.3 fir n<l 
1 1 6+ Sp. 5 fir n>1 
fir n>1 

= 926,40 

0,7 — 3,63 43,15 = = 258,90 
0,8 — 4,89 18,07 — = 108,42 
0,9 — 9,73 9,71 = ne 58,26 

oo k _k co k _k 
(Cy a (Kay a 
1.0 m7+ >! a Tay mn7+ > pee ay 38,52 
if (Pesel a 
=—0,57 — In 9,22 - 10-4=6,43 | 

if + 10,49 4,63 10,64 4,78 28,68 
12 + 5,56 3,70 5,58 | 3,72 2232 
1,3 + 3,94 3,07 3,95> | 3,09 18,54 
1,5 + 2,76 2.46 2,68 | 2,39 14,34 
1,75 + 2,12 2,00 2,07 | 1,95 11,70 
2,0 + 1,80 1,74 177 17d 10,26 
25 + 1,54 1,52 1,49 | 1,46 8,76 
3,0 + 1,39 1,38 1,866 1,35 8,10 
20 — 1,00 1,00 1,00 6,00 


Fir n = 1,5; 1,75; 2,0; u. 3,0 ist die Summe sogar gri®er, das ist nur deshalb der Fall, weil die 
Entwicklung mit dem 19. Glied abgebrochen wurde. Wir sparen die ganze mihselige Auf- 


summierung der Reihe fiir n>1, wenn wir die Funktion JJ ome verwenden, deren 
nm 
Funktionswerte bequem aufgeschlagen werden kénnen. Diese Tatsache beweist 


wiederum nur, da es villig unerheblich ist, ob wir 7 oder co als obere Grenze annehmen. Dem- 
nach lauten die fiir die Berechnung zweckmafigsten Oberflachenformeln fir n Sil: 
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ayn ye 6519) 221 0% 
Bs 1 it | 
6108 ff | -= - Saget 
O= ; v| | y "e  dy— i 4 
d aa | 
m™ 
eke 
fe P 10° ne, be 1 (9,22 be Oma) n 13 
ies seods | wets 1 2 a 
it 
n 
| 1 9,92 10-4! a | 
O, = 6.- 10? . =~ oo, (9, ) fiir n>1 
: d | n 1 ? 
tee 
1m 
b) fiir din = 0,1 6 und Ymin = all 
eget _— (Onl ete asl . 
O, = 6-10 7 |Z = ae | ibang es es (14) 
| poe a) 
m™ 


Tabelle 2 (nach Methode IT fiir dmin = 0,1 pn). 


Spalte 7 Spalte 8 Spalte 9 
ee — aah They 1 =| ep ts for n>1 
i y > ea ‘ 1 Con 
ae Verges, if —-— (d@’) 
0.6 —4,01 + 3,77 (d’)%4 = 
0,7 — 3,63 + 4,65 (d’)%3 = 
0.8 —4,89 + 6,34 (d’)%? = 
0,9 — 9,73 + 11,34 (d’)% aa 
18 k k 
1,0 In7 4 Da —. : = In 0,1-+In d/=+1,73+In d’ SS 
Ii 
8,78 
iil — + 10,64 ~ (aynt 
3,78 
1192 — 5,08 (d’)o2 
2,17 
ees i 3,95 (d’)%8 
0,95 
1,5 ar ae 2,68 ~~ (d’/)o8 
0,41 
1,75 = al 2,07  (d’)%75 
0,20 
2,0 — + 1,77 — dv 
0,05 
2,5 a + 1,49 (ayes 
0,02 
3,0 a + 1,36 (dp 
oO — | + 1,00 
1 {Sp. 8 fir n<1) d’ in ps 
On, = 6+ 108+ 3; 


d’ \Sp. 9 fiir n>1) Ox, in m?/kg 
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Fiir n<1 miissen wir die Reihenentwicklung durchfithren, dann lassen sich die Gleichungen 
(11) und (12) nicht vereinfachen. Nimmt man fiir die untere Grenze des Integrals anstatt Ymin = 


/ 


rnin n 0,1) n 
9,22 - 10—* (entsprechend einem Riickstand R= NY fal a Rog (= ) = ( : 


d’} d 
so tritt dort an die Stelle des Gliedes 
1 
= 
= n n—l 
(Qp22 1008) a erin 9 222 10-4 seule eee 
if 1 (a= 
ego ieee 
Tm m1 


bzw. (In 0,1 — Ind’). Somit ergibt sich fiir die Oberflachen die endgiltige Wertetabelle 1 fiir dic 
Berechnung nach Methode J und die Tabelle 2 fiir die Berechnung nach Methode II (mit dyin = 
0.1 4). Die Multiplikation der Spalte 8 bzw. 9 mit 6 - 103/d’ ergibt die gesuchten Oberflachenwerte 
O;,, nach Methode II, die ersichtlich bedeutend umstandlicher als die Methode Jist. Spalte 8 zeigt, 
daB O,, fir n = 0,6 und d’ = 1 w negativ wiirde, was natiirlich physikalisch sinnlos ist. 

Dies hat seinen Grund in folgendem: Wenn bei der Auswertung des Integrals in Form einer 
Reihenentwicklung fiir die untere Grenze y,,;, nur das erste Glied beriicksichligt wurde, dann war 
Voraussetzung dafiir, daB ymin = (dmin/d’) < list. Wird nun d’ = 1 y, so ist diese Voraussetzung 
nicht mehr erfiillt, und die Vernachlassigung der héheren Reihenglieder wirkt sich dahin aus, daB 
man einen negativen Wert fiir das Integral erhalt. Fiir d’ = 10, 100, 1000 wv usw. ist jedoch die 
Voraussetzung gut erfillt und die nach der Tabelle gerechneten Werte stimmen aufs beste mit den 
graphisch von S. KieSkalt! umgerechneten Werten iibercin. Beide unterscheiden sich natirlich 
(und besonders fiir kleine n betrachtlich) von den nach der Methode I (entsprechend einem Riick- 
stand R von 99,9%) berechneten und gezeichneten Werten. In Tabelle 3 sind die Werte von Ox, 
fiir d’= 10, 50, 100, 500 und 1000 uw nach beiden Methoden zum Vergleich eingetragen. 


Tabelle 3. 
Sas 
ae 10 50 100 500 1000 
BS 
0.6 3264 1670 1181 498 B35 
92,640 18,528 9264 1853 926 
0.7 3390 1365 892 319 200 
| 25,890 5178 2589 518 259 
0.8 3096 1070 661 204 122; 
10,842 2168 | 1084 PA UT 108 
0,9 2730 850 4.96 137 Ue 
5826 1165 583 Tiley 58.3 
1.0 2418 677 380 95,4 51,8 
3852 770 385 17,0 38,5 
iil 2184 564 306 al Sitti 
2868 574 287 57,4 28,7 
1194 1914 462 245 54 Pa lesll 
2232 446 POPS} 44,6 22,3 
3 1716 394 204 43,4 D2 
1854 371 185 aunt 18,5 
Wes 1428 306 155 31,6 15.9 
1434 287 143 28,7 14,3 
175 1198 246 11733 24.9 12,4 
1170 234, 117 23.4 Junker 
2.0 1050 741 9 106 Zila 10.6 
1026 205 103 20,5 10,3 
2.5 894 179 89.4 17,9 8,9 
876 175 | 87.6 es 8,8 
3.0 816 163 81,6 16.3 8,1 
810 162 81,0 1052 8,1 
600 120 60.0 12,0 6.0 


F 2/4, /Oberer Wert nach Methode II entsprechend dmin = 0,1 
OG am tore \Unterer Wert nach Methode I entsprechend ‘ 


1 S. KieSkalt, Bauwirtschaft Heft 1 (1949). 


max = 99.9%. 
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6. Graphische Darstellung der spezifischen Oberfliche als Funktion von d’, Aus der Form der 
Gleichungen (12) und (13) ist unmittelbar ersichtlich, da® sich in einem logarithmischen Netz mit 
d’ als Abszisse und Ox, als Ordinate die Oberflaichen O; als Geradenschar mit dem Scharparameter n 
Haretellen. Alle diese Geraden fir verschiedene n sind einander parallel und haben eine Steigun 
von 135” gegen die positive d’-Achse, denn es ist z. B. nach (13), wenn wir den Faktor pas 


1 


m( s) [9.22- 10-4)» 
os 
nN 


O;,==; log 0, = log C,— log d’. 


cena a! 
ieee ante 


mit C, bezeichnen, 


Selbstverstandlich gilt diese Uberlegung nur fiir die untere Grenze Yin = 9.22 - 10-4. Fite die 
untere Grenze Y;,j, = 0,1/d’ erhalt man hingegen keine reinen Geraden mehr, wenn auch das Glied 
0,1" /(1 — 1/n) -1/d’™—! , sehr klein werden kann (namlich fiir hohe n). Je gréfer also n wird 
um so geradenahnlicher werden die Oberflichenkurven sein (vgl. Abb. 1). 


ni/igt spez.Obertlache tur Kugelkorn °=7 u.8=1000 kg rt, Oe, spez. Oberfische Ox, ' 
70 T T = 
8 + 
6 Cot 2 d 
cae a | {TTT 


| 06 | | 
p 07 B 2 
j 08 
70 = a“ 10° 
gt H 09. HH 8 
6 44 a 
10 | 9 
4g a 4 
| | ial 17 im 
42 | 
2b + ll 43 3 
15 
| 415 
70" 38 19" 
8 iE a = 30 alias ara |e 
6 [ i pee ; 
4 | ieee Ea Hf y 
2 2 
ae JE 7 1 5——_1— ae feline Arline | eae tel be 
2 4 68 PD IT EG ie tip = F 4 6 bw 
Abb. 1 


7, SchluBbemerkung. Die vorliegende Arbeit zeigt die Bedeutung einer hoéheren transzen- 
denten Funktion, der Gammafunktion, fiir siebtechnische Probleme. Die Methode der Auswertung 
eines Integrals durch Reihenentwicklung und die Verwendung der Gammafunktion fihren hier 
zur Liésung eines Problems, das fiir Fillstoffe auf dem Kunststoffgebiet, Baustoffe, Pigmente 
der Lack- und Farbenindustrie, Adsorbentien usw. von praktischem Interesse ist. Die numerische 
Berechnung der Oberflachenwerte ist deshalb notwendig, weil es bislang keine befriedigenden ex- 
perimentellen Methoden zu ihrer Bestimmung gibt. 


(Eingegangen am 6, Juni 1950.) 
Anschrift des Verfassers: Dr. G. Matz, Leverkusen-Wiesdorf, Gellertstr. 5. 
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Systematische Abwindmessungen an Pfeilfliigeln*. 
Von H. Trienes und E, Truckenbrodt, 


1, Einleitung. Die Anzahl der Arbeiten, die sich mit der Bestimmung des Abwindes in der 
Umgebung eines Tragfliigels befassen, ist sehr groB. Sowohl in experimenteller wie auch theo- 
retischer Hinsicht wurde das Problem des ungepfeilten Fliigels in zufriedenstellender Weise ge- 
lést. Eine Zusammenfassung aller diesbeziiglichen Arbeiten bis zum Jahre 1938 geben I. Fliigge- 
Lotz und D. Kiichemann}. In zwei weiteren Arbeiten amerikanischen Ursprungs? werden eben- 
falls MeBergebnisse mit der Theorie verglichen. Von den zahlreichen theoretischen Arbeiten, 
die alle auf der Anwendung der Prandtlschen Tragfliigeltheorie und des Biot-Savartschen 
Gesetzes beruhen, wollen wir nur die Arbeiten von H. Multhopp? und E. Truckenbrodt* erwabnen. 

In neuerer Zeit hat nun der gepfeilte Fliigel stark an Bedeutung gewonnen. Uber syste- 
matische Sechskomponentenmessungen und Druckverteilungsmessungen an Pfeilfliigeln hat 
kiirzlich W. Jacobs® berichtet. Diese Arbeit vermittelt einen recht guten Uberblick iiber die am 
Fligel auftretenden Krafte und Momente, sowie die Auftriebsverteilung in Spannweiten- und 
Tiefenrichtung bei symmetrischer und unsymmetrischer Anstrémung. 

Das Ziel der vorliegenden Arbeit soll sein, tiber Abwindmessungen hinter denselben von 
W. Jacobs benutzten Fliigeln zu berichten. Die Bedeutung dieser Messungen besteht darin, Unter- 
lagen dariiber zu liefern, wie sich der Beitrag des Héhenleitwerkes zur Stabilitat um die Quer- 
achse mit dem Pfeilwinkel des Flugzeuges andert. Durch die vom Fliigel abgehende Wirbelflache 
werden bekanntlich hinter dem auftrieberzeugenden Fliigel Abwartsgeschwindigkeiten hervor- 
gerufen, die den wirksamen Anstellwinkel des Héhenleitwerkes verkleinern. Bei jedem Fliigel 
hangen diese Abwindwinkel auBer vom Abstand des Héhenleitwerkes zum Tragfliigel noch 
sehr stark von der am Fliigel herrschenden Auftriebsverteilung ab; und zwar ist der Abwind- 
winkel am Ort des Leitwerkes um so gréBer, je gréBer das Verhaltnis des drtlichen Auftriebs- 
beiwertes in Fliigelmitte zum Gesamtauftriebsbeiwert des Fliigels ist. Ist der Tragfliigel nun ge- 
pfeilt, so andern sich einerseits die Abstandsverhaltnisse und zum anderen beeinfluBt die beim 
gepfeilten Fliigel geanderte Auftriebsverteilung auch stark die Abwindwinkel am Ort des Leit- 
werkes. Beim riickwarts gepfeilten Fliigel ist bei gleichem Gesamtauftrieb der Auftrieb in Flii- 
gelmitte kleiner als beim ungepfeilten Fliigel, was eine Verkleinerung des Abwindwinkels zur 
Folge hat. Das bedeutet aber, daB die wirksamen Anstellwinkel des Héhenleitwerkes vergréBert 
werden, und damit das Leitwerk zur Stabilitaét um die Querachse einen gréBeren Beitrag liefert 
als beim ungepfeilten Fliigel. Bei nach vorn gepfeilten Fliigeln wird die stabilisierende Wirkung 
des Leitwerkes entsprechend verkleinert. Ahnlich kann man sich auch iiberlegen, daf bei gleichem 
Abstand beim stark zugespitzten Fliigel das Héhenleitwerk erheblich weniger stark stabilisierend 
wirkt als bei einem Rechteckfliigel. Die Richtigkeit dieser Angaben wird im Verlauf der Arbeit 
durch unsere Messungen bestitigt. 


2, Durchfiihrung und Umfang der Messungen, Zur Bestimmung des Abwindes in der Um- 
gebung eines Fliigels kénnen grundsatzlich zwei Methoden verwendet werden. Die eine be- 
stimmt den Abwind punktweise mittels eines WinkelmeBgerates, z. B. Zwei- oder Vierfingergerat!.2. 
Dieses Verfahren hat den Vorteil, daB es auBer den brtlichen Abwindwinkeln auch die Staudruck- 
verteilung mitliefert. Es hat aber dafiir den Nachteil, daB es recht mithsam und langwierig ist. 


* Die Messungen stellen einen Auszug aus der von der Fakultat Maschinenwesen der Technischen 
Hochschule Braunschweig angenommenen Dissertation von H.Trienes dar (1946). (Hauptberichter: Prof. 
Dr. H. Schlichting, Mitberichter: Prof. Dr. H. Blenk). Sie wurden 1944/45 ausgefiihrt im Windkanal des 
friiheren Aerodynamischen Institutes der T, H. Braunschweig (Leitung Prof. Dr. H. Schlichting). Die 
Bearbeitung fiir die Veréffentlichung wurde von E. Truckenbrodt vorgenommen; yon ihm stammt auch 
der Abschnitt 4 tiber den Vergleich mit der Theorie. 

1 [, Fliigge-Lotz und D, Kiichemann, Jahrbuch 1938 der Deutschen Luftfahrtforsch., S. I 172. 

» NACA Report No. 648 und 651 (1939). 

° H. Multhopp, Luftfahrtforschung 15 (1938), S. 463. 

4 EF. Truckenbrodt, Ing.-Arch. 18 (1950) eSe2s08 

® W. Jacobs, Ing.-Arch. 18 (1950), S. 344 (I. Mitt.) und 19 (1951), S. 83 (IT. Mitt.). 
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Die zweite Methode, die wesentlich schneller zum Zicle fiihrt, bestimmt den Abwindwinkel mit 
Hilfe einer hinter oder vor dem Fliigel befindlichen ,,Fuhlflache“, die mit einem einfachen Ge- 
stange am Fliigel befestigt ist. Nach diesem Verfahren wurden die nachstehend mitgeteilten 
Messungen durchgefithrt. Der angewendeten Methode liegt folgende Uberlegung zugrunde: Es 


werden Dreikomponentenmessungen (Auftrieb, Widerstand und Kippmoment in Abhangigkeit 
vom Anstellwinkel) durchgefiihrt am 


Modell ohne Fiihlflache (Fliigel +- Ge- fufttlache 

stange) und am Modell mit Fiihlflache Flige/ aay grolbe tndscheibe 
(F liigel la Gestinge + Endscheiben Bezugspunkt=geometrischer Hh Meine Lndschelte 
-++ Fihlflache), vgl. Abb.1 und Abb. 2. Neutrajounkt io 


Wenn man die Messungen mit und ohne Hugelsehne 


; z GS 
Fithlflache vergleicht, dann gibt es | 
Me 
~~ Oo z | ——, 
—— — —| 
“ uinie 

geomerrischer 8 

< Neutrajpunkt oF 

on SI9q 


rintiiote =P 75 
i oe /nle 


a ey 


“| Gestinge om 
lige! verschraubt 


(oe in mm) 


by269 


3b! 


Qy= aE b=150 = 


Abb. 1, Anordnung der Fihlflache (Rechteckfligel). Abb. 2. Modellskizze des ungepfeilten Rechteckfliigels, 


einen ganz bestimmten Anstellwinkel des Fliigels x*, bei dem das Kippmoment mit und ohne 
Fihlflache gleich ist. Das bedeutet aber, daB® bei diesem Winkel die Fiihlflache unbelastet ist, 
d.h. der Anstellwinkel der Fiihlflache «= 0 ist!. Wie aus Abb. 3 hervorgeht, besteht zwischen 
dem Anstellwinkel des Fliigels a 


und dem der Fiihlflache ay, der Zu- Nullauttriebsrichtung 


sammenhang , Cog hunlflache 
ae _ Anstromungsriohtung. “b= 

On = O-> x — Oy. (1) Fliigelsehne Hep Hiilicke pte ek a 
Hierin bedeutet eg den Winkel i 5 ES i (aie 
zwischen der Fliigelsehne und der ve va Bee ie | 
Nullauftriebsrichtung der Fiihl- ~/ acess 
flache?, der auch als Schrankungs- a Ly el 
winkel bezeichnet wird. Die GroBe Flige/ Fuh Hache 


QO, ist der mittlere Abwindwinkel am 
Ort des Leitwerkes (der Fiihlflache). 
Er enthalt die Beeinflussung des 
Leitwerkes durch den Fliigel. Wir 
rechnen den Abwindwinkel positiv, 
wenn er den Anstellwinkel des Leit- 
werkes verkleinert. 
winkel zu 


Ow = 0 + Ex. 


Abb. 3. Erlaiuterungsskizze. 


‘Winkel am Fliigel und an der Fihlfliche. 


a = Anstellwinkel des Fliigels 
Or = Anstellwinkel der Fiihlflaiche (Héhenleit werk) 
&» = Abwindwinkel am Ort der Fihlflache (Hohenleitwerk) 


ey = Schrankungswinkel = Winkel zwischen Fliigelsehne und Null- 
auftriebsrichtung der Fiihlfliche (Héhenleitwerkssehne) 


Bei gegebenem ey, findet man aus &dy—0 den gesuchten Abwind- 


(2) 


1 Dieses ist richtig, wenn man annehmen darf, daB keine Riickwirkung von der Fiihlflache auf den 


tok) 


Hauptfliigel stattfindet. 


2 Hat die Fihlflache symmetrisches Profil, dann fallt die Nullauftriebsrichtung mit der Sehne zusam- 


men, 
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Die Auswahl der zu messenden Anstellwinkel wird so getroffen, da sich jeweils eindeutige 
Schnittpunkte fiir die Momentenkurven cy(a) mit und ohne Fihlflache ergeben. 

Die Abwindmessungen wurden an einer Reihe von unverwundenen gepfeilten Trapezfliigeln 
gleichen Profiles (NACA 23 012) und annahernd gleichen Seitenverhaltnisses (A = 5, ohne End- 
kappen) durchgefiihrt. Erginzend wurde noch ein unverwundener ungepfeilter Ellipsenfligel 
gleichen Profiles hinzugenommen. Die Abmessungen der Fliigel gehen aus Tabelle 1 sowie Abb. 4 
hervor. 


Tabelle 1. Geometrische Fliigeldaten. 


; sin Tief Tief Mittl 3) it it u. ; 
Zuspitzung Pfeilwinkel ooeion ee ices | chuesBndkeppes Fligelflache Seitenverhaltnis 
l, 
Z= — P ki lq lm b | b/ F 3 
My 


[m*] 


0.2 0° 45° 0,250 | 0,050 | 0,172 0,756 | 0,750 0.1128 5,06 
0.6 0° 0,188 | 0,113 | 0,153 0,764 | 0,750 0,1135 5,16 
1.0 a! 0,150 0.150 0,150 0,768 0,750 0.1146 5,18 


Abb. 4. Ubersicht der gemessenen Fliigel. 


Die Momentenbezugsachse geht jeweils durch den geometrischen Neutralpunkt des Fliigels. 
Unter dem geometrischen Neutralpunkt verstehen wir den Auftriebsschwerpunkt einer Fliigel- 
halfte bei konstanter ¢, ;,,-Verteilung und unter der Annahme einer drtlichen Neutralpunktlage 
auf 25% der értlichen Fligeltiefe. Fiir einen Trapezfliigel fallt dieser geometrische Neutral ae 
mit dem //4-Punkt im Flachenschwerpunkt einer Fliigelhalfte zusammen. : 

Der Momentenbeiwert ist folgendermafen definiert: 


ee 
‘ qF lm ° (3) 
b/2 
Hierin bedentet J,,— 4 [ 12(y) dy die e . e ; ‘a 
mR geometrische mittlere Fliigeltiefe. Fiir den Trapezfliigel ist 
—b/2 
L, gleichbedeutend mit der Tiefe im Flachen r kt ei P - : 
fiir J,, sind in Tabelle 1 angegeben. SUSE Pia tase cancers alent eee vacnee 
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verwendet. 


Die Spannweite der Fithlflache by betragt +/, der Spannweite des Hauptfliigels ohne Endkappen 


Systematische Abwindmessungen usw. 


Trienes u. Truckenbrodt: 
Als Fithlflache wird ein Rechteckfliigel mit symmetrischem Profil (NACA 0015 
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(Abb.2); es ist somit b= 0,250 m. Die Flache hat 15% der Flache des Rechteckfliigels. Die Fiih]- 
flache ist mit Endscheiben versehen. Diese Endscheiben tragen Bohrungen, die eine Einstellung 
der Fiihlfliche bei der kleinen Endscheibe in den Hochlagen 6,,= 0 und —0,2 erméglichen; werden 
die Endscheiben umgekehrt angebracht, so erméglichen sie eine Einstellung der Fithlflache in 
den Hochlagen Cy = 0 und --0,2. Die groBen Endscheiben wurden benutzt fiir die Hochlagen 
Cu = =F 0,4. Die Endscheiben sind auf ein Gestange aufgeschoben und in jeder beliebigen Riick- 
lage hinter dem Fliigel einstellbar. Das Gestinge selbst besteht aus Profilstangen und ist mit 
dem Fliigel fest verbunden. 
Der Umfang der durchgefiihrten Messungen geht aus Tabelle 2 hervor. 


Tabelle 2. Umfang der durchgefiihrten Messungen. Die Riicklage &y der Fiihlflache ist gemessen 


Raktila ohowe vom geometrischen Neutralpunkt des Fliigels 


Ricklage Hochlage Schrankungswinkel bis zum Einviertelpunkt der Fithlflache. 
; _ ty ‘on 2 sy r Die Messungen wurden im Windkanal des 
a HY i fritheren Aerodynamischen Institutes der T. H. 
| ———— Braunschweig (Strahldurchmesser 1,2 m)_ bei 
==(),5 a) ae : ; einer Windgeschwindigkeit von v= 40 m/sec 
jeweils 6 Winkel im daschestahteaD eT ales R ld 
0 Bereich gefiihrt. Das entspricht einer Reynoldschen 
150 0,2 —9°<en <6° Zahl fiir den Hauptfliigel von 
0,4 Re =a y= 4 10 
—1,5 0 Bei der Auswertung wurden Anstellwinkel und 
Zz Widerstand fiir den endlichen 
ee ee — 8" Seep Strahldurchmesser nach den iib- 
| PAP Aer by F5 by= lichen Methoden? korrigiert. Fiir 
| pes oy oe das Lingsmoment wurde keine 
i 6° res: Korrektur vorgenommen. 
=69 =-40 =0 , ; 
-30 3. Ergebnisse der Messungen, 
is H ap Fiir das Beispiel des ungepfeilten 
+ =-15| =0 — Rechteckfliigels ist in Abb. 5 das 
fe Ey 709 fo | -p0 Ergebnis einer vollstandigen Mes- 
“3 50 =-90| =- G2 sung dargestellt. Die ohne Mef- 
iyo Zo | 8° ee punkte eingetragene Kurve ist 
= 72 | jeweils die Momentenkurve ohne 
ae 50 | =-40|=-0¥ Wihlflache. Fiir jede Anordnung 
Le OF yy Fol po] eae Aas ergeben sich entsprechend den 
r 70 he =70 ie ah sechs verschiedenen Schrankungs- 
es a Z winkeln sechs Schnittpunkte zur 
oe | Bestimmung von a*. Aus diesen 
Ko 5 9°77 40 g ergibt sich nach (2) der mittlere 
for 70 2 . Abwindwinkel &,,. Wir erhalten 
pra aa? aI ea ae 7 He auf diese Weise Auftragungen 
6-5 Loo ght O,(0)?. Fiir das Beispiel des un- 
oe . gepfeilten Rechteckfliigels zeigt 


Abb. 6. Abwindwinkel in Abhangigkeit vom Anstellwinkel @,,(c) fiir den 


ungepfeilten Rechteckfliigel. 


Abb.6 das Ergebnis. Wir haben 
an dieser Stelle darauf verzichtet, 


fiir alle Fliigel die den Abb.5 und 6 


entsprechenden Auftragungen anzugeben und verweisen hierfiir auf die Originalarbeit von 


H. Trienes. 


Aus den Auftragungen 4,,(«) haben wir die Steigungen da,,/da bestimmt und diese in Abb. 7 
iiber der Riicklage —&y fiir die Hochlage €, = 0 und in Abb. 8 iiber der Hochlage Cy fiir &;—=— 1 
aufgetragen. Wir haben jeweils drei Auftragungen vorgenommen: 


a) Zuspitzung Z— 1,0; y veranderlich; b) y= 0; Z veranderlich; c) 


= 45°; Z verdinder- 


lich. Die Kurven in Abb. 8 zeigen eine systematische Unsymmetrie bei entsprechenden posi- 
tiven und negativen Werten Cy. Die Werte oberhalb der &-Achse (¢q<0) sind immer etwas 


eréBer als diejenigen unterhalb der &-Achse. 


1 Ergebnisse der Aerodynamischen Versuchsanstalt Gottingen, IT. Lieferung. 


2 Der Einfachheit halber haben wir fortan geschrieben a* = a. 
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Bereits in der Einleitung sind wir auf die flugtechnische Bedeutung der Abwindmessungen 
eingegangen. Es sei an dieser Stelle noch kurz angegeben, in welcher Form die von uns gemessenen 
Abwindwinkel in die Momentenbilanz des Flugzeuges eingehen. 

Ist uns der Abstand des Héhenleitwerkes (Einviertelpunkt) hinter dem Flugzeugschwerpunkt 
mit ry gegeben, dann erzeugt das Leitwerk den Momentenbeitrag 


My =— Con du Fury. (4) 
Der Index H gibt an, daB-es sich um die 
GréBen des Hohenleitwerkes handelt. 
Ist uns der Auftriebsanstieg des Leit- 
werkes de,7/day bekannt, dann kénnen 
wir, wenn wir noch den dimensionslosen 
Beiwert nach (3) einfiihren, fiir (4) auch 
schreiben 


au 


ties 3 da Fy Ty Ie, 7 


m dX py ets (6) G2 


In dieser Gleichung ist uns der Anstell- 
winkel des Leitwerkes a, durch (1) ge- 


geben. DenBeitrag des Héhenleitwerkes 9 9G 06 76 7 75 rE. “3 
zur Stabilitat um die Querachse erhalt ORS ; ; 
man durch Differentiation von (5) nach aie 
dem Auftriebsbeiwert des Gesamtflug- 5 
zeuges. Dieses ergibt die Verschiebung * p=0 
des Gesamtneutralpunktes durch das ; EES 
Leitwerk go 
deny ae (- = nu 10 
de, L. = 
(Pe pee Ole da 
cs Set a OH Tat H da ; (6) he 
q F 1, do, da de, > 
A | al B Fi/ipse 

Wir haben dabei angenommen, daf das : ¥ 

sig it H 
Staudruckverhaltnis qu!4 vom Anstell- og 3 78 7 E a + 
winkel unabhangig ist. Die entschei- - 
dende GréBe in (6), die den EKinfluB oad L 
des Fliigels auf das Leitwerk wieder- 06 


gibt, ist auBer im Staudruckverhaltnis 
qu/q im Wert day/dx enthalten, den 
wir unter Beachtung von (1) wie folgt o¥ 
aus unseren Messungen finden: 


do da, 


ie Lane (7) gg 


Wenn wir jetzt annehmen, da sich 
das Leitwerk jeweils vom geometrischen 
Neutralpunkt aus gemessen gleich weit "ng 06 
hinter demFliigel befindet, dann kénnen Abb. 7. Der Abwindanstieg da,,/dox in Abhangigkeit von der Riicklage 
wir aus Abb.7 und 8 in Ubereinstim- Stee? bei-dermtisablageel yoy = 6. 

mung mit dem in der Kinleitung 

Gesagten folgendes beziiglich des Stabilitatsbeitrages des Leitwerkes feststellen: 

1) Bei gleicher Zuspitzung (Z = 1,0 und Z = 0,2) riickt der Neutralpunkt bei zunehmender 
Riickwartspfeilung nach hinten. Fir die gemessenen Fliigel vom Seitenverhaltnis A = 5 ist die 
Neutralpunktverschiebung bei einem um 45° nach hinten gepfeilten Fliigel gegentiber einem 
ungepfeilten Fliigel um etwa 25% gréfer. Bei Vorwartspfeilung liegen die Verhaltnisse um- 
gekehrt. 

2) Bei gleicher Pfeilung (p = 0 und y = 45°) ergibt das Leitwerk, wenn es sich hinter dem 
am meisten zugespitzten Fliigel befindet, die kleinste Neutralpunktverschiebung. Sie betragt 


08 40 12 1¥ 16 


- 02 


OY 
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bei einem Seitenverhaltnis von A= 5 fir den Fligel Z = 0,2 etwa 25% weniger als beim Recht- 

eckfliigel (Z = 1,0). : 
Wie bereits gesagt, wurden samtliche Messungen gewonnen mit Fliigeln vom Seitenverhalt- 

nis A ~5. Wir kénnen den EinfluB des Seitenverhaltnisses naherungsweise bescitigen, wenn 


GH Z4=10 Cu 
—OY 
r) 
{+ p45 | Y= a 
30° Sa” 
15° ae —b2 
aiky 
La 
0 
G2 Ke G6 
Y 
7 lige 
a b 
GY 2 
Abb. 8. Der Abwindanstieg dx,,/da in Abhangigkeit von der Hochlage ¢y bei der Riicklage E47 = —1; A = 5. 


wir die gemessenen Ab- 
windwinkel durch den in- 
duzierten Anstellwinkel 
bei elliptischer Auftriebs- 
verteilung, d.i. ¢,/7A 
(c,= Gesamtauftriebsbei- 
wert) dividieren. Wir 
schreiben also 


Trienes (Fige/ mit Gestinge) oy aa A F 
Jacobs c,/% A da de,/dax ° (8) 


Hierin bedeutet de, /da 
den Auftriebsanstieg des 
jeweiligen Fliigels. In 
Abb.9 haben wir diese 
Werte nach den Messun- 
gen von H. Trienes (Fliigel 
+ Gestange)und nach den 
Messungen von W. Jacobs! 
miteinander verglichen. 
Wir haben das arithmeti- 
sche Mittel aus beiden 
Werten genommen und 
damit die Werte von 
x A/(de,/dx) nachTabelle 3 
bestimmt. Die hiermit 
nach (8) aus Abb. 7 und 8 
ermittelten Werte von 
O»/(C,/% A) sind in Abb. 10 
iiber dem Zuspitzungsver- 
haitnis Z aufgetragen; 
und zwar im oberen Teil 
fiir die Riicklage 
fn = 1,0 
und die Hochlagen 


Abb. 9. Auftriebsanstiege dc,/da der verschiedenen Fliigel 
a) in Abhangigkeit vom Pfeilwinkel y. b) in Abhangigkeit von der Zuspitzung Z. Cn ars 0, + 0,2 und 0,4; 


1 Siehe FuBnote 5 von S. 26. 


XX. Band 1952. Trienes u. Truckenbrodt: Systematische Abwindmessungen usw. 33 


und im unteren Teil fiir die Riicklage £y = —1,5 und die Hochlage Cy = 0. Mit Annaherung 
an den Rechteckfliigel (Z = 1,0) nehmen die Abwindwinkel ab. Zum Vergleich haben wir auch 
noch die MeBwerte des Ellipsenfliigels miteingetragen. In Abb.11 haben wir die Werte 
Ow/ (Cai A) itber dem Pfeilwinkel — fiir die Riicklage £;— —1 und die Hochlagen 6, = 0; = 0,2 
und - 0,4 dargestellt. Mit gréBer werdender Riickwartspfeilung nehmen die Abwindwinkel ab. 


4, Vergleich mit der Theorie. Tabelle 3. Werte von 2A/(dc,/do). 
Wir wollen jetzt die experimentell 
gefundenen Werte mit theoretischen wii 
Rechnungen vergleichen. Die von 
E. Truckenbrodt! angegebene Nahe- 2 
rungsformel zur Berechnung des fh 
Abwindes von ungepfeilten Fliigeln, Ellipse 
wollen wir auf den gepfeilten Fliigel 
ausdehnen. Dabei wollen wir uns auf den Fall beschranken, wo die Aufpunkte (&), ¢) sich inner- 
halb der Wirbelschicht (€ = 0) befinden. Zur Berechnung des Abwindwinkels einer beliebig ge- 
formten tragenden Wirbellinie greifen wir auf eine Formel von H. Multhopp? zuriick, die wir fiir 
unsere Zwecke wie folgt schreiben: 


=] 
: E 7 L(é, 7, 7 , a 
au (Es) = elim 2 ( rye f : 2) (a) a. (9) 
ers 


oar) 


Los 


0, 
0 


Hierin bedeutet y(7) = I'(n)/vb WES 
die dimensionslose Auftriebsver- 6, [A 
teilung langs Spannweite und 
L(é, 7,7’) eine Einflu8funktion, 
die sich zu 


L(E,9,4')=1 


past = e710) 

ee) G7") 
berechnet. Dabei sind §,7 die 
Koordinaten des Aufpunktes, und 
&’ = &'(n’) ist die Lage der tra- 

genden Linie (Abb. 12). 

Nehmen wir einen Pfeilfliigel 
an, bei dem die tragende Linie 
(Einviertelpunktlinie) auf jeder 
Fliigelhalfte cine gerade Linie 
bildet, dann kénnen wir fiir die 

EinfluBfunktion schreiben 


L(é,7,n')=1 (11) 7 (Tessung 
So i 11/1). te Theorie (S= 0) 


 VWEe—G, I) tegl? + — 7’ 


Die neuen Bezeichnungen €, und 
Na gehen aus Abb. 12 hervor. 
Unter 7, wollen wir den seitlichen pT etn ent ee 
Abstand des Auftriebsschwer- C7 A 
punktes einer Fliigelhalfte (aero- dae ce ak ae 
dynamischer Neutralpunkt) verstehen*®; er berechnet sich zu 


1 
ees 
Na = oie 
en 
1 Siehe FuBnote 4 von S.26. 7? Siehe FuBnote 3 von S. 26. 
3 Wir wollen im Folgenden immer annehmen, daf die Auftriebsverteilung zu 7 = 0 symmetrisch ist. 


3) 


Smee, é 
— in Abhangigkeit von der Zuspitzung Z; 


Se (12) 
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Entsprechend soll auch die Lage des Aufpunktes in €-Richtung (¢,) vom Neutralpunkt aus ge- | 


messen werden. 


Wir entwickeln jetzt Gleichung (11) in eine Taylorsche Reihe nach tg y und finden bei Vernach- | 


lassigung von Gliedern zweiter Ordnung 


, ca G1) Gay 
L¢(é, Y, 1 ) = 1 ; ae oS (13) 
il GieuP) ye Ga) 
Nehmen wir jetzt an, daB & > (j7—7’)?, dann kénnen wir noch weiter vereinfachen und finden 
, Pre heey ; 1p aia : 
LE. 9) =o tie ee Pe ee ee (14) 
rales Gas fi 
7 Wir setzen diesen Ausdruck in (9) ein und er- 
o, it A halten nach einer entsprechenden Rechnung, : 
48 . ; A ay, a 
j wie sie schon in der zitierten Arbeit von 
A= 50, 2=40, Cy= "10 E. Truckenbrodt durchgefiihrt worden ist, das 
fa Gt | folgende Ergebnis: 
- dsc oat ata 
+| 2 = aN Wy, We = 
= Ow (&.%) 0 zi} Oj () 4 eae TA P. (15) 
eS 
~ a 
oi vi Aufounkt 
12 OY, (aes 
48 E<0 ou 
4 —  Messung aerodynamischer 
eee Tove (c-0) DS 
Ce =k | 7 
re tragénde Linke 
| ' Na, CR 
=) =A eae 0° 10° Jy ff 49° 0°D Wi | \! 


Ow 


Abb. 11. Die Abwindwerte 


in Abhangigkeit vom Pfeil- 


Abb. 12. Erlauterungsskizze. 


a 
winkel g; Rechteckfliigel, Z = 1,0; A = 5. 


Hierin bedeutet 4;(7) den nach Prandil bekannten induzierten Anstellwinkel am Ort des Fliigels 
und c, den Gesamtauftriebsbeiwert des Fliigels. Unter Beachtung von (12) kénnen wir uns 
leicht davon iiberzeugen, da das letzte Glied in (14) keinen Beitrag zum Abwindwinkel 
liefert. Fiir den ungepfeilten Fliigel geht (15) in die bereits bekannte Formel iiber. Der maB- 
gebliche Abstand ist also vom Auftriebsschwerpunkt (Neutralpunkt) bis zum Aufpunkt zu 
nehmen. Nach unseren Messungen ist €y der Abstand vom geometrischen Neutralpunkt bis 
zum Kinviertelpunkt der Fiihlflache. Ist dic Lage des aerodynamischen Neutralpunktes 7, be- 
kannt, dann errechnen wir die Verschiebung gegenitber dem geometrischen Neutralpunkt 7, zu 


AS == (a — Ha tees (16) 


Bei den gemessenen Rechteckfliigeln ist 7, stets kleiner als 7,; daraus folgt also, daB der Auftriebs- 
schwerpunkt vor dem geometrischen Neutralpunkt liegt. 
Aus der erweiterten Traglinientheorie (Dreiviertelpunkt-Theorie) ist bekannt, daB zur Fest- 


legung des Anstellwinkels, im vorliegenden Fall des Anstellwinkels der Fihlflache, der Punkt auf 


75% der Fliigeltiefe maBgeblich ist (Pistolesisches Theorem). Die von uns gemessenen Abwind- 
winkel gehéren also zu den Riicklagen 


ly 
§ = Ey — = Eg — 0,091. (17) 


Fiigen wir noch den Anteil aus (16) hinzu, dann erhalten wir fiir den in Formel (15) bendtigten 
Abstand 8 


a= Su — | — (Ns — Na) 189 - (18) 
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Den Vergleichswert des Abwindwinkels mit der Messung (£,< 0) erhalten wir als Mittelwert 
tiber die Fiithlflachenspannweite. Nach (15) ergibt sich dann 


No 
er ee 
areal | oi(n) dn + 42 nA’ (19) 
wobei 7 = 2 yo/b = 1/3 das Spannweitenverhiltnis der Fiihlflache zum Hauptfliigel darstellt. 


Zur Berechnung von q; ist die genaue Kenntnis der Auftriebsverteilung erforderlich. Nehmen 
wir an, da uns diese langs Spannweite durch den Fourierschen Reihenansatz 


ve \ ,; | 
Y= 5 = 2 D>) Gen yi sin (2n+ 1)0 (20) 


n=0 


gegeben ist, wobei # durch die Beziehung 


2 
n=—~ = cosd (21) 


definiert ist, dann kénnen wir den induzierten Anstellwinkel in bekannter Weise berechnen. 


Bes (Up Palen Eke (22) 


sin? 


Fihren wir jetzt unter Beachtung von (21) die Integration nach (19) aus, dann erhalten wir, wenn 
wir noch beriicksichtigen, daB sich der Gesamtauftriebsbeiwert durch Integration der Zirku- 
lationsverteilung (20) iiber die Spannweite zu 


i 
tA {yen — ala, (23) 
Sil 
berechnet, 
Ow o i cos(2m-+1)0 %n+1 | ia 
c,/tA 4 2 cos By ce 4&2 8) 


Fiir unseren Fall mit 4) = 1/3 und entsprechend @ = 70,5° gilt 


Xw ee 5 en si as ° 
Sed = 22000 — 5,116 3 + 5,948 38 + (24) 


Auf eine theoretische Behandlung des Problems, wenn die Aufpunkte auBerhalb der Wirbel- 
schicht liegen, soll an dieser Stelle verzichtet werden. Die in (24a) auftretenden Fourier-Koeffi- 
zienten entnehmen wir fiir den ungepfeilten Fliigel einer unveréffentlichen Arbeit von H. Schlich- 
ting und W. Kahlert, diejenigen des gepfeilten Rechteckfliigels einer Arbeit von W. Mutterperl?. 
In beiden Fallen wurde die Auftriebsverteilung nach einer erweiterten Traglinientheorie (Drei- 
viertelpunktmethode) berechnet. Ein Vergleich der theoretischen Auftriebsverteilung nach 
Mutterperl mit den Messungen von Jacobs zeigt befriedigende Ubereinstimmung. Tabelle 4 gibt 


eine Zusammenstellung der benétigtenZahlen- Tabelle 4. Fourier-Koeffizienten der Auftriebsverteilung. 
werte. 


In Abb. 10 haben wir die nach (24a) er- Geig4 45/4, | Ns—Na 
mittelten theoretischen Werte fiir §, = — 1,0; —0,0651 0,0193 
—1,5 und —oo dargestellt. Es zeigt sich, 0,0125 0,0151 — 
daB die Ubereinstimmung bei dem Rechteck- 0,0581 0,0050 
fliigel nahezu vollkommen ist. Bei den zu- 0,0581 0,0050 0,061 
gespitzten Fliigeln (Z — 0,2) ergeben sich eek Dyas ae 
einige Abweichungen, die wohl im wesent- 071533 | —0,0010 0,036 


lichen darauf beruhen, daB die erweiterte 
Traglinientheorie Auftriebsverteilungen liefert, die fiir unsere Zwecke nicht genau genug sind. 
Wir folgern dieses aus der Tatsache, daB fiir &,;— —1,5 die gemessenen Abwindwinkel kleiner 
sind als die theoretischen fiir £;;== — oo, was theoretisch nicht méglich sein dirfte. 


1 W. Mutterperl, NACA Tech. Note 834 (1941). 
2 Interpolierte Werte. 
3% 
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In Abb. 11 zeigt die theoretische Kurve mit zunehmender Pfeilung gréfere Abweichungen 
gegeniiber der MeBkurve. Bei g = 45° betragt der Fehler etwa 20%. Der Grund hierfiir muf 
wohl auch darin gesehen werden, da die zugrundegelegte Auftriebsverteilung des Pfeilfliigels 
den Genauigkeitsanspriichen einer Abwindrechnung noch nicht véllig geniigt. 


5, Zusammenfassung. Es wird iitber Abwindmessungen hinter gepfeilten Fliigeln berichtet, 
die mittels der Fiihlflachenmethode ausgefiithrt wurden. Die Messungen vermitteln einen guten 
Uberblick iiber den Beitrag des Héhenleitwerkes zur Stabilitat um die Querachse bei gepfeilten und 
zugespitzten Fliigeln. Befindet sich das Leitwerk in gleichem Abstand vom geometrischen Neu- 
tralpunkt hinter dem Fliigel, dann ist der Beitrag des Leitwerkes zur Stabilitat um die Querachse 
um so gréfer, je starker der Fliigel bei gleicher Zuspitzung nach hinten gepfeiltist. Bei Vorwarts- 
pfeilung ergibt sich das umgekehrte Ergebnis. Je gréBer die Zuspitzung ist, um so kleiner wird 
bei sonst gleicher Pfeilung der stabilisierende Beitrag des Héhenleitwerkes. Die Messungen werden 
mit theoretischen Rechnungen verglichen. 


(Eingegangen am 13. Juni 1951.) 


Anschrift des Verfassers: Dr.-Ing. E.Truckenbrodt, Braunschweig, Technische Hochschule, 
Institut fiir Str6mungsmechanik, 
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Turbulente Stroémungen in konvergenten Kanalen. 
Von W. Szablewski. 


1, Einleitung. Fiir dié scheinbare durch den Vorgang der turbulenten Mischung erzeugte 
Schubspannung hat L. Prandtl! im Jahre 1925 den bekannten hypothetischen Ansatz 


t= 0 12| | ou 

dy | Oy 

mit dem Mischungsweg l aufgestellt. (Bei Betrachtung einer Grenzschichtstrémung bedeuten 
dabei in ttblicher Weise x die Bogenlange der Kérperkontur, y den Wandabstand, u, v die Ge- 
schwindigkeitskomponenten in x- und y-Richtung, @ die Dichte, js die Zahigkeit, y= lo die 
kinematische Zahigkeit, p(x) den statischen Druck, t die Schubspannung.) In zwei neueren 
Arbeiten hat der Verfasser? auf der Grundlage des Prandtlschen Ansatzes die turbulente Stré- 
mung im Rohr und langs der ebenen Platte in allen charakteristischen Einzelheiten in guter Nahe- 
rung berechnen kénnen. In der vorliegenden Untersuchung werden als weiteres wichtiges Stré- 
mungsmodell die turbulenten Strémungen in konvergenten Kanalen behandelt. Eine rationelle 
Theorie dieses Modells ist, soweit dem Verfasser bekannt, bisher noch nicht gegeben worden. 
Wir schildern zunachst die in den Arbeiten iiber das Rohr und die ebene Platte zur Anwendung 
gekommene Methode des Rechenganges. 

Die Differentialgleichungen der turbulenten Grenzschicht 


Go) Ow ldp , 1% Gn , Ce __ 
eee T sional o dx | 0 Oy’ anya 


werden von der Wand her in zwei Schritten integriert. Die so gewonnene Funktion fiir die Ge- 
schwindigkeitsverteilung dient im itbrigen Bereich der Grenzschicht als Approximation. Wand- 
schubspannung Tt, und Grenzschichtdicke 6, die in der Funktion in Form von Parametern auf- 
treten, werden dann mittels der Randbedingung u = u, (u, Potentialstrémung am Rande der 
Grenzschicht) und des Impulssatzes bestimmt. 

Die genannten zwei Integrationsschritte kénnen folgendermafen charakterisiert werden: Der 
erste Integrationsschritt ist, wie in Abschn. 2 noch naher ausgefiihrt werden wird, gekennzeichnet 
durch die Annahme konstanter Schubspannung t ~ T,) und in bezug auf den turbulenten Mi- 
schungsvorgang durch den modifizierten Prandilschen Ansatz 


L = 2 (y — 0p) (=~ angendhert) 
(6) Dicke der laminaren Unterschicht). 


Der erste Integrationsschritt erstreckt sich iiber die laminare Unterschicht, den Bereich 
gleicher GréBenordnung von laminarer und turbulenter Reibung bis in den Bereich der (iiber- 
wiegend) turbulenten Reibung hinein. Fiir diesen letzteren Bereich nimmt die berechnete Ge- 
schwindigkeitsverteilung die bekannte Gestalt des logarithmischen Verteilungsgesetzes an. 
Die so erhaltenc Lisung wird dann im zweiten Integrationsschritt weiter in den Bereich der 
turbulenten Reibung der Grenzschicht hinein fortgesctzt. Dieser zweite Schritt ist gekennzeichnet 
durch Vernachlassigung der laminaren Reibung, in bezug aut den Vorgang der turbulenten Mi- 


schung durch den Ansatz 


2 
bent (SY 
Bei diesem zweiten Schritt wird dem Einflu8® der Tragheits- und Druckkrafte, die beim ersten 
Schritt entsprechend der Annahme Tt ~ To vernachlassigt wurden, in erster Naherung Rechnung 
getragen. a - 

Als empirische Koeffizienten treten in der analytischen Entwicklung, dem derzeitigen Stand 
der Theorie gemaB, die Koeffizienten x, 1 der Mischungswegvertcilung auf, ferner die dimensions- 


los gemachte Dicke der laminaren Unterschicht v,,69/7, wobei vy, = tle die sogenannte Schub- 


1 [. Prandtl, Z. angew. Math. Mech. 5 (1925), 5. 137. ' : 
2 W. Szablewski, Z. angew. Math. Mech. 31 (1951), 5. 131 sewie 5. 309. 
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spannungsgeschwindigkeit ist. Diese Koeffizienten besitzen leider keinen universellen Charakter, 
sondern sind als von der Zahigkeit abhingige Gré®en anzusehen. Immerhin wird fiir sehr groBe 
Reynoldssche Zahlen nach den Messungen Nikuradses' der ZahigkeitseinfluB sehr gering. Des 
weiteren ist eine Abhangigkeit vom Druckgradienten der Grenzschicht denkbar. Sehr wesent- 
lich ist hier das Ergebnis, daB die Koeffizienten des ersten Integrationsschrittes und v,. do/? 
vom Druckgradienten der Grenzschicht unabhangig sind! Den experimentellen Nachweis fiir 
diesen Sachverhalt haben neuerdings H. Ludwieg und W.Tillmann? erbracht. Fiir groBe Re- 
Zahlen kann man demgem48 diese Koeffizienten als ,,universelle“’ Konstanten ansprechen. 

Bei der Untersuchung der hier behandelten turbulenten Strémungen in konvergenten Ka- 
nalen ergibt sich nun, daB bereits mittels des im ersten Integrationsschritt zu gewinnenden loga- 
rithmischen Geschwindigkeitsgesetzes alle wesentlichen Merkmale der Strémung in guter Ap- 
proximation berechnet werden kénnen. (Es ist anzunehmen, daB sich entsprechend einfach all- 
gemein turbulente Grenzschichten mit Druckabfall behandeln lassen.) Bei den verzégerten Stré- 
mungen in divergenten Kanalen wird es allerdings erforderlich scin, den oben erwahnten zweiten 
Integrationsschritt zu vollziehen. 

Als experimentelle Unterlagen dienen der Untersuchung die bekannten Messungen von J. Ni- 
kuradse®, Leider — in Hinsicht auf eine zu entwickelnde Theorie der turbulenten Grenzschichten 
fiir groBe Re-Zahlen — liegen diese Messungen bei verhaltnismafig kleinen Reynoldsschen Zahlen 
(Re = u, b/y ~ 70000, b Kanalbreite). 


2, Das logarithmische Geschwindigkeitsgesetz in der Reibungszone turbulenter Grenzschichten‘, 
Wahrend im allgemeinen in der Grenzschicht einer Strémung definitionsgema4B die Tragheits- 
krafte von derselben GréBenordnung sind wie die Reibungs- und Druckkrafte, stirbt nach der 
Wand hin zu die Bewegung wegen der Haftbedingung ab. Dagegen nimmt die Reibung hier 
infolge des sehr steilen Geschwindigkeitsanstieges betrachtliche Werte an. In einer der Wand 
anliegenden schmalen Zone werden wir daher die Tragheitskrafte gegeniiber den Reibungs- und 
Druckkraften vernachlassigen kénnen. Die Bewegungsgleichung einer Grenzschicht 


Ou Cte ee tisdp 1 or 
"ant Vay ede + e 
reduziert sich dann fiir diese Zone auf 
CDOT ae 
eee a oyna 0, (11) 
aus der 
d 
t=%+ i y (III) 


folgt. Hierbei ist 
: T = Tham =e Tturb 
mit (IV) 
Vien == fe Tigh == 0 aay 
lam b ay ’ turb — OQ dy 5 
Beschranken wir uns auf hinreichend kleine Wandabstande, so kénnen wir, T) + 0 voraus- 
gesetzt, noch in (III) den Einflu8 des Druckgradienten vernachlassigen und 


TT, (V) 


setzen, also konstante Schubspannung annehmen®. Bei hinreichend kleinen Wandabstanden 
ergibt demnach die erste Integration von (I) 


pe tor) =n. (VI) 


ote ( icneie Gesetzmafigkeiten der turbulenten Strémung in glatten Rohren. VDI-Forschungsheft 
2 H, Ludwieg und W.Tillmann, Ing.-Arch. 17 (1949), S. 288. 
* J. Nikuradse, Untersuchungen tiber die Strémungen des Wassers in konvergenten und divergenten Ka- 
nalen. VDI-Forschungsheft 289 (1929). 
4 Die Ausfiihrungen dieses Abschnitts decken sich teilweise mit den Darlegungen von J. Rotta. Inc.- 
Arch, 18 (1950), S. 277. is ee 
5 Diese Feststellung hat erstmalig L. Prandtl bei der Behandlung der turbulenten Rohrstré getroff 
Ergebn. d. AVA, Gottingen, IV. Lfg. (1932), S. 18. ara th eh 
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Fir den Mischungsweg J kann man nach Prandtl in Wandnihe aus dimensionsanalytischen 
Griinden 


ta= 4 
setzen. Traégt man dem Umstand Rechnung, daB auch turbulente Grenzschichten in unmittel- 


barer Nahe der Wand einen laminaren Strémungsbereich, die sogenannte laminare Unterschicht, 
besitzen, so ist dieser Ansatz zu modifizieren in 


l= x(y — 6). (VII?) 


Die weitere Integration ]aBt sich dann unschwer durchfiihren. Man erhilt eine Geschwindigkeits- 
verteilung, die innerhalb der gekennzeichneten Zone fiir ,,nicht zu kleine‘‘ Wandabstinde zu 
einem logarithmischen Geschwindigkeitsgesetz fithrt, wie es zuerst L. Prandtl? fiir das Rohr und 
die ebene Platte abgeleitet und fiir wandnahe Punkte experimentell bestatigt gefunden hat. Das 
logarithmische Geschwindigkeitsgesetz lautet 


# la > + Ve|—-1+m (4 x Res Vey) + x Res Ve, 3 


uy 

fiir 
1 y 
ee a 
T u,0 

—— Re; = + 

(« Qu} f = (VIII) 
oder auch 

Se Pg ele lg Agee 9c 8% 

Vx x v x v 
fiir 


ita Ux ¥ 
Tu << Reser. 


In einer neueren Arbeit haben H. Ludwieg und W.Tillmann® den experimentellen Nachweis des 
logarithmischen Geschwindigkeitsgesetzes wandnaher Punkte auch fiir allgemeine Grenzschichten 
mit beliebigem Druckgradienten erbracht, 

Wir wollen die betrachtete Zone der Grenzschicht die Reibungszone der Grenzschicht 
nennen. Sie ist durch annahernd konstante Schubspannung und in bezug auf den Mischungs- 
weg durch 1 ~ x (y — 6) gekennzeichnet. 

Analysiert man die Grenzschichtstrémung vom Standpunkt der GréBenordnung der lami- 
naren und turbulenten Reibung, so gelangt man zu einer weiteren, fiir die Theorie der turbu- 
lenten Grenzschichten bedeutungsvollen, Feststellung. Wir kénnen hier bekanntlich in der 
Grenzschicht drei Bereiche unterscheiden. Auf die laminare Unterschicht folgt ein Ubergangs- 
bereich mit gleicher Gré®enordnung der laminaren und turbulenten Reibung, der in den Bereich 
der (iiberwiegend) turbulenten Reibung miindet. Laminare Unterschicht und Ubergangsbereich 
bilden beide nur sehr schmale Zonen der Grenzschicht, der Bereich der turbulenten Reibung 
erfiillt fast die ganze Grenzschicht. Fiir die theoretische Behandlung von Bedeutung ist nun 
die Feststellung, da die in (VIII) genannte Einschrankung 


y 1 
ee eh Le 
fiir den Giiltigkeitsbereich des logarithmischen Geschwindigkeitsgesetzes sich als gleichbedeutend 
mit der Vernachlassigbarkeit der laminaren Reibung erweist! Integriert man (VI) unter der 
genannten Vereinfachung, so erhalt man wieder das logarithmische Geschwindigkeitsgesetz. Man 
kann den Sachverhalt auch formulieren, indem man feststellt, dafB die oben definierte Rei- 
bungszone der Grenzschicht sich bis in den Bereich der turbulenten Reibung hinein erstreckt. 
Uber die im logarithmischen Geschwindigkeitsgesetz (VIII) auftretenden empirischen Ko- 
effizienten haben wir schon in der Einleitung gesprochen. Sie hingen von der Reynoldsschen 
Zahl Re;— u,6/y des Strémungsvorganges ab, wenn auch diese Abhangigkeit fiir sehr groBe 


1 Vgl. hier und im folgenden W. Szablewski, a. a. O. 
2 L. Prandtl, a.a.O., s. FuBnote, S. 38. 
3 H. Ludwieg und W. Tillmann, a. a. O. 
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Re-Zahlen eine sehr geringe wird. Dagegen hangen sie nicht, wie H. Ludwieg und W.Tillmann' 
experimentell nachgewiesen haben, vom Druckgradienten der Grenzschichtstrémung bzw. der 


Widerstandsziffer 4 = — dp/dx - 6/(@/2)u% ab. 


3. Bewegungsgleichungen, Die konvergenten und divergenten Kanalstrémungen stellen ein 
wesentlich allgemeineres Modell der Strémungsmechanik dar als Rohr und ebene Platte. Wah- 
rend beim Rohr im voll ausgebildeten Strémungszustand die Trigheitskrafte in Fortfall kommen, 
kénnen wir bei der ebenen Platte konstanten Druck annehmen. Demgegeniiber treten bei den 
Strémungen in verengten und erweiterten Kandlen sowohl Druck- wie Tragheitskrafte auf. 

Wahlen wir das Koordinatensystem wie in Abb. 1, so lauten die Bewegungsgleichungen der 
Grenzschicht bei Beschrankung auf kleine Offnungswinkel 


Qu, du\ dp | | Ou , {| du| du) ; 
o(u t+ 5y)= dx | Maya + ay?" |ayl ay’ (1) 

Ou dv 

at sy=?- (2) 


Im Sinne der Grenzschichttheorie waren die Reibungskrafte durch die Ableitungen normal zur 

Kanalwandung, und nicht in Y-Richtung, zu bestimmen. Fiir kleine Offnungswinkel ist dieser 

Unterschied vernachlassigbar. Der Druck wird als 

Z Y konstant iiber dem Kanalquerschnitt angenommen. 

Diese Annahme ist von Nikuradse? experimentell als 

A ir /-LeeN gutreffend nachgewiesen worden. Theoretisch kénnte 

— = “ sie an Hand der zweiten Bewegungsgleichung fiir die 
Y-Richtung nachgepriift werden. 


oe 


Yor tgoun | b/2=tgoe-x # Wir transformieren auf das Strahlenbiischel 
(tg 0¢<0) (tga>0) “7 
konvergenter Kanal — divergenter Kanal n= 1+ oe 2, (3) 
(0 =Kanalbre/te) ice 2-3 
Abb. 1. wobei 7=0 dann mit der unteren Kanalwand zu- 
- sammenfallt. 
Mit dem Ahnlichkeitsansatz fiir den voll ausgebildeten Strémungszustand 
u= Uy (x) p (1) » v= 0, (x) y(n) (4) 
(u, (x) Geschwindigkeit in Kanalmitte) erhalten wir aus (1) 
b 1 du, . | dy 1, dp dp b/2 D Che d {/ 1 \2 /dy\2 
2 u, dx eo n) Pan oo Uy uy dj dxqu? ! u, b/2 dn? ! dn (ep) (7) | ; 2) 
Es folgt 
1 : 
ea es Vy ~ Wy (~ proportional). 


Mit dem Ansatz u,; = A/x, v, = B/x ergibt (5), wenn wir noch die Reynoldssche Zahl 


Re=— ae 
und die Widerstandsziffer 
_ BP, 
ae dx 
SCHATE 
einfiihren, One 
d i al A Dind= 1 Ll \2 /dq\2 
==} “m2 t -( al = WY g We 
BOSE E De AY ay tore die Peale) (Gr) | 2) 
Aus der Kontinuitat (2) folgt 
A 
p=— xtee-(L—n)¢. (7) 
Durch Elimination von yw aus (6) mittels (7) erhalt man schlieBlich die Gleichung 
R 2eo “2 ae 
v _— Be 
4 | Re dr? | dy (ea) i) = eR re (8) 


1 H. Ludwieg und W. Tillmann, a. a. O. 
2 J. Nikuradse, a. a. O., s. FuBnote 3 von Seite 38. 
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Integration iiber 7 ergibt 


ho Diy LM fap. — n 
Teale ae tle) e =— tga [ g dy, (9) 


dy j 


0 Out Quy 


2 dy To ree To 
ne ay (zr) =—. cosa ~ —. 
(bei Beschrankung auf Winkel cos « ~ 1). 


Die Randbedingungen lauten g = 1 und = 0 fir 7— 1. Mittels dieser Randbedingungen 
ergibt sich aus (9) fiir 7) = 1 sofort der Impulssatz 


Mae 
a= pt tga: (ody. (10) 


Grenzt man als Kontrollflache ein Kanalsegment der Lange Eins und der Hohe b/2 ab, so besagt 
der Impulssatz, da die Impulsanderung, welche die Strémung beim Durchtritt durch das Seg- 
ment in der Zeiteinheit erfahrt, gleich den auf die Kontrollflache ausgeiibten Druck- und Rei- 
bungskraften ist. 

Welcher Bereich der Kanalstrémung besitzt nun den Charakter einer Grenzschicht ? Ist nur 
ein Teilbereich als Grenzschicht anzusprechen, so haben wir jedenfalls in Kanalmitte Potential- 
strémung, so da 


bzw. 
igg=— —. (11) 


Mit dieser Beziehung lautet der Impulssatz (10) 
1 
= tea! fo dy —1 |. (12) 
0 


Da der physikalischen Bedeutung nach der Ausdruck in der Klammer negatives Vorzeichen hat, 
¢, positives Vorzeichen, so folgt aus (12), daB tga <0 sein muf. D.h. nur bei konvergenten 
‘Kandlen haben wir in einem der Kanalmitte anliegenden Bereich Potentialstrémung zu erwarten, 
wahrend bei divergenten Kanalen der gesamte Strémungsbereich den Charakter einer Grenz- 
schichtstrémung hat. Die Beziehung (11) ist nun in der Tat fiir konvergente Kandle nach den 
Messungen von Nikuradse? erfiillt: SE 


a tg a —Ala 
Im folgenden werden wir nur noch die konvergenten Stré- | | ——— 
mungen behandeln. Mit der dimensionslosen Dicke 75 der Grenz- ser ee ae 
schicht cae wir aus (9) fiir 7= 7,5 mittels der Randbedingungen at ae 07140 = 0139 
gp ~ 1 und i ~ 0 fir 4= 75, und wenn wir noch (11) beriick- 
1 
sichtigen, 
ns 
a= teal | gy? dn — | - (13) 
Fihren wir den Winkel a* durch 
6 = (tg a — tg a*) x 
ein, so kommt nach (3) 
ee 14 
(OE ase (14) 
Transformieren wir auf die Koordinate 
z a 15 
S fii ae tg a* ” Ce) 
eee 


1 J. Nikuradse, a. a.O., s. FuBnote 3 von Seite 38. 
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so schreibt sich (13) 
a= (ga —teat)| foe as — 1]. (16) 
0 


In unserem Formelsystem treten zwei unbekannte Parameter auf, die Schubspannungsziffer 
cy(= 2 ¢,) und der die Grenzschichtdicke 75 festlegende Winkel tg a*. Zu ihrer Bestimmung 
werden wir im weiteren Verlauf der Rechnung die Randbedingung y = 1 fiir 7 = y5 und den 
Impulssatz in der Form (16) heranziehen. (Wir ziehen bei unserer Naherungsrechnung den Im- 
pulssatz statt der zweiten Randbedingung dp/dyj = 0 fiir 7 = 75 heran, da es mehr verspricht, 
eine approximierende Funktion zu integrieren statt zu differentiieren.) 


4, Berechnung der konvergenten Kanalstrémung. a) Logarithmisches Geschwindig- 
keitsgesetz. Wir beschrainken uns gleich auf das logarithmische Geschwindigkeitsgesetz, das 
fiir den turbulenten Teil der Reibungszone gilt (siehe Abschn. 2). Dieses lautet nach Formel 
(VIII), wenn man beachtet, dai Re = 2 Rey), ist, 


a — 
p =n + ]— 14 mn Ox Re Ye) +e), (17) 
v 

oder mit 

log. Verfeilung mit #=06 ter Vig vxy _ Reley (18) 

Ux Vey mT | 

auch 
Sia a Best a d Oo 
; _ in ea 1+ In4x+% ale (19) 


Hier treten die empirischen Koeffizienten x und 
v* 6,/v auf. Den Koeffizienten % ermitteln wir aus 
der Auftragung (u,;—u)/v* der MeBpunkte iber 
In 7 (vgl. Abb. 2). Nach (17) wird 


uy —Uu 


=——Iny + Konst. (20) 


Man ersieht aus Abb. 2, dafS§ das logarithmische 
a OS 10 15 LO 46 Geschwindigkeitsgesetz fiir kleine 7 erfillt ist. 
Abb. 2. Experimentelle Geschwindigkeitsverteilungen Des weiteren ergibt sich, was fiir die Theorie von 

nach Nikuradse, VDI-Forschung 289 (1929). A ; ;: 
Bedeutung ist, dai die Neigung der Geraden un- 
abhingig von « bzw. dem Druckgradienten ist. Es ergibt sich 


x ~ 0,60. (21) 


Den Koeffizienten v* 6,/y ermitteln wir aus der Auftragung u/v* iiber v*y/y der MeBpunkte, 
Abb. 4. Den theoretischen Zusammenhang gibt (19) an. Man findet fiir die konvergenten Stré- 
mungen 


a | P40 | Re 
— 2° 9,60 69 000 
— 4° 9,60 67 000 
— 8° 9,30 66 000 
also v*59/y ~ 9,50 fiir Re ~ 67500. (22) 


Fiir die divergenten Strémungen ergeben sich, bei Re-Zahlen derselben GréBenordnung, 
gréBere Abweichungen. Nach den neueren Messungen von H. Ludwieg! und W. Tillmann ist 
auch v*/d)¥” als vom Druckgefalle unabhangig anzusehen. 

Wir verwenden nunmehr die — exakt fiir den turbulenten Teil der Reibungszone giiltigen — 
Verteilungsfunktionen (17) und (19) als Approximation fiir die Verteilung der Geschwindigkeit 
im ganzen Bereich der Grenzschicht. 


1 Vel. H. Ludwieg und W. Tillmann, a. a. O. 
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Unterwerfen wir (vgl. den Hinweis Abschn. 3) die Verteilungsfunktion (17) der Randbe- 
dingung 


Cea fiir ai bzw. oe 


so erhalten wir die mit einem gewissen Fehler behaftete Relation 


_ Vex icy = by 
T= “Ain m+ —|—1- In (2% Re yey) + x=). (23) 
Wir kénnen dann fiir die Verteilung (17) auch schreiben 
p=1+!eing. (24) 
b) Grenzschichtdicke. Setzen wir (24) in den Impulssatz (16) ein, so erhalten wir 
Se a eee 25 
Sade (25) 
% 


bzw. nach (14) 


(26) 
Nach (18) wird 
(27) 
| 
ke=67 500 
e ie dag Mikuradse Re 67500 
V. -Forschung 269 (1929) 70 © ae"? Experiment Mikuradse 
a COE o — -g9 l-Forschung 289 (1928) 
| G —/nheorie 
Ug ed 
0 0 0 rr) r) ° o eS ay he 
Q eZ. -4 -6 -8 -10 O OF _ &o 15 20 ie) G0 Oe 40 
Abb. 3. Abb. 4. 


Mit den in c) angegebenen theoretischen Werten von ¢, l= = “) ergeben obige Formeln fiir die 


Grenzschichtdicke 


ye) 
Hains ine 


2,626 


— 2" 0,460 | 2,883 
| 2,071 


Inwieweit die Rechnung zutrifft, mag man nach Abb. 4 und 5 beurteilen. 
c) Widerstandsgesetz. Aus (23) und (26) ergibt sich die Widerstandsformel 


is : eee: os Cy 1 1 The meee (28) 
bere? ( = i tg (—a)| 
o 
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Mie; —— ae Man erhalt 


e | of theor | Caap 
— 2° 0,00483 | 0,00456 
— 4° 0,00534 | 0,00522 
— 8° 0,00592 | 0,00634 


(Vgl. Abb. 3). Die experimentellen Werte wurden dabei von Nikuradse’ mittels des Impuls- 
satzes aus der gemessenen Staudruckverteilung berechnet. 

Die Ubereinstimmung ist recht gut. Um so mehr, wenn man beachtet, daB der MeBwert cy 
fiir « — — 8° relativ sicher zu hoch liegt. Man ersieht das aus Abb. 4: die MeSpunkte der ex- 
perimentellen Geschwindigkeitsverteilung fiir « — — 8° liegen hier im Bereich des logarithmischen 
Geschwindigkeitsgesetzes unterhalb der MeBreihen fiir a = — 2° und a= — 4°. Neuere Mes- 
sungen2 haben jedoch ergeben, daB die Koeffizienten des logarithmischen Geschwindigkeits- 
gesetzes unabhangig vom Druckgradienten sind. 

d) Geschwindigkeits- 
verteilung. In Abb.4 haben 
GesthwintgKeltsverrEHlung wir das logarithmische Ge- 
schwindigkeitsgesetz (19) mit 
den empirischen Koeffizienten 
(21) und (22) bis zur Grenze 
| | V4, 6/v, die wir in (27) berech- 


net haben, dargestellt. Fiir das 


Re ~67 500 Restintervall bis zum End- 
° &= 2) Lyperiment Nikuradse punkt 4 
0 pe VAI- Forschung 289 (1929) — 0,,(b/2)/y = Reye,/2 
Theorie | haben wir dann Potentialstr6- 
mung und konstanten Verlauf 
| | | | | y UO 1/ye, (Cres. nach c)) 
if 7" tp anzunehmen. Auf die ab- 
: ee ee eS ee ao : 
0 OI G2 OF fy 05 06 07 eb 09 = weichende Lage der MeBwerte 
ieee fiir « = — 8° haben wir bereits 


oben hingewiesen. 
Die Geschwindigkeitsverteilung in natiirlichen Koordinaten enthalt Abb. 5. Die theoretische 
Verteilungsfunktion lautet nach (24) unter Beriicksichtigung von (26) 


Vex 
[ee 
p= ling ! F ! ina(2 * tg ( a) (29) 


x x Vey 
(Citheor nach c)). Man sieht, daB die theoretischen Verteilungen verbesserungsfahig sind. Eine 
bessere Approximation kénnte dadurch erreicht werden, dafs man die Tragheits- und Druck- 
krafte in erster Naherung beriicksichtigt®. Dabei ware es dann auch erforderlich, die Mischungs- 
wegverteilung in zweiter Naherung entsprechend dem Ansatz 


Sanit pe 


zu erfassen, was mit 4 die Hinfithrung eines dritten empirischen Koeffizienten mit sich bringen 
wiirde. Da jedoch bereits mittels des logarithmischen Geschwindigkeitsgesetzes die wesentlichen 
Ziige der hier behandelten Strémung in guter Naherung herauskommen, verzichten wir hier 
darauf. 

e) Verdrangungsdicke. Formel 


1 ts 


6 

sa | Ao an~ [ Lao) dy (30) 
0 0 

J. Nikuradse, a. a.O., s. FuBnote 3 von S. 38. 


1 
2 H. Ludwieg und W. Tiilmann, a. a. O. 
3 Vel. W. Szablewski, a. a. O. 
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lautet bei Transformation auf & nach (15) 


Oy tg a* = 
a -() aa (1—) dé. (31) 
0 
Mit (24) erhalt man bei Beachtung von (25) 
Oy Cy 
= ae 2 
b/2 Ve (32) 
OAS [em a) tg (— a) 
Py 
Mit Ci4.0, nach c) ergibt das a 
5 aaa Ye 
Oy Oy 
babes (Fa)exp oe Verdrangungsidicke 
° Ke 67500 
ae ms pe Hee an ° Experiment Mikuradse 
3° 0.012 ’ VOI -Forschung 289 (1929) 
= ; ON 4030 
(Vgl. Abb.6). Auch diese charakteristische Strémungs- ees 
eigenschaft kommt theoretisch sehr gut heraus. 0020 ; 
AbschlieBend bemerken wir, daB man nach dem 
hier entwickelten Formelsystem selbstverstandlich ggg 
auch konvergente Kanalstrémungen bei groBen Re- 
Zahlen berechnen kann. Bei der turbulenten Stré- ou : : : = &, 
mung langs der ebenen Platte haben wir fiir sehr Div S24 nian sll Ga Ee 
groBe Re-Zahlen x ~ 0,39, v, do/y © 5,521. Abb. 6. 


5. Zusammenfassung. Die Geschwindigkeitsverteilung in dem Teil der konvergenten Kanal- 
strémung, der den Charakter einer turbulenten Grenzschicht hat, wird durch das logarithmische 
Geschwindigkeitsgesetz approximiert. Die empirischen Koeffizienten dieses Gesetzes werden der 
Messung entnommen und erweisen sich als vom Druckgradienten unabhangig. Mittels einer 
Randbedingung und des Impulssatzes werden dann die Wandschubspannung sowie Grenzschicht- 
und Verdrangungsdicke berechnet. Gute Ubereinstimmung mit den Messungen. 


(Eingegangen am 14. Juni 1951.) 


Anschrift des Verfassers: Dr. W. Szablewski, Berlin W 8, JagerstraBe 22/23. 
Deutsche Akademie der Wissenschaften, Institut fiir Mathematik, Abt. ,,Angewandte Mathematik‘. 


1 Vel. W. Szablewski, a. a. O. 
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Wiarmespannungen in einer Scheibe infolge einer wandernden 
Warmequelle. 


Von E. Melan. 


Im Folgenden sollen die Warmespannungen in einer diinnen Scheibe, die sich allseits ins 
Unendliche erstreckt, ermittelt werden, wie sie durch eine lings einer Geraden mit gleichfér- 
miger Geschwindigkeit wandernden, punktférmigen Warmequelle entstehen. 

Das Material sei sowohl hinsichtlich der thermischen wie der elastischen Eigenschaften iso- 
trop; alle das Verhalten des Materiales festlegenden GréBen seien konstant, also weder von der 
Temperatur noch von der Spannung abhangig. LaBt man diese vereinfachenden Annahmen 
gelten, so bietet die Spannungsverteilung in einem hinreichend grofen Blech, die bei der Her- 
stellung einer SchweiBnaht entsteht, ein praktisches Beispiel fiir diese Aufgabe. 

Voraussetzung fiir die Bestimmung der Warmespannung ist in allen Fallen die Kenntnis der 
Temperaturverteilung. Bei einem Kérper von unendlicher Ausdehnung und gleichférmiger Ge- 
schwindigkeit der Warmequelle langs einer Geraden tritt ein quasistationdrer Temperaturzustand 
auf, der wieder eine quasistationare Spannungsverteilung zur Folge hat. Fiir einen mit der glei- 
chen Geschwindigkeit v wie die Warmequelle sich bewegenden Beobachter wird sich weder der 
Spannungszustand, noch die Temperaturverteilung andern. Bezieht man sich daher auf ein mit 
der Warmequelle wanderndes Koordinatensystem, wie dies im Folgenden geschieht, so kommt 
in den Gleichungen die Zeit t nicht vor. 

D. Rosenthal und H. Cameron haben sich mit der Temperaturverteilung bei beweglichen 
Warmequellen befaBt und bei einem quasistationéren Zustand die Differentialgleichung fiir die 
Temperatur T 

EOD BOT Sten oe ee eO de 
A ii "ay? | az? a 0x 
abgeleitet. Darin bedeuten x, y und z die Koordinaten eines Punktes in einem mit der gleichen 
Geschwindigkeit v wie die Warmequelle wandernden Bezugssystem, a = k/oa, wobei k die innere 
Warmeleitfahigkeit, 9 die Dichte und o die spezifische Warme vorstellt. Nimmt man die Tem- 
peratur langs der geringen Dicke der Scheibe konstant an, so verschwindet 0T/0z, und T ge- 
nigt mit v/2a== A der partiellen Differentialgleichung 
BP ES aT 
BN 5 T ay? BY a 
Das von Rosenthal angegebene partikulare Integral dieser Gleichung 
T= Ne-*i M(t) GFSyey®, i=y—1), 
worin Hj die Hankelsche Funktion erster Art und der Ordnung Null vorstellt, geniigt bereits 
den Randbedingungen, daB namlich 0T/dr fiir x > co und dT/dy fiir y > co verschwindet. Die 
Integrationskonstante N ist aus der Bedingung zu bestimmen, daB die Warmequelle in x = 0 
und y = 0 die Intensitat q besitzt, daB also fiir r— 0 


oT 
wird. Nun ist lim ir 2 Hi (iAr) = — =. und damit erhalt man 
r>0 


T= q/4k e—4*i Hi (iAr). 

In einer fritheren Arbeit? hat der Verfasser gezeigt, wie bei einer diinnen Scheibe mittels des 
thermischen Verschiebungspotentials @ die Spannungen infolge eincr Temperaturverteilung T 
ermittelt werden kénnen. © geniigt der partiellen Differentialgleichung 

#6) #0 m+ 
Ox? dy? am —] 


ole 


1 D. Rosenthal, Trans Amer. Soc. mech. Engr. 68 (1946), S. 849, sowie Rosenthal und C 
(1947), S. 961, ferner Max Jakob, Heat Transfer oa a. York. ae arr aaa 


* E. Melan, Osterr. Ing.-Archiv Bd. IV (1950), S. 153 und ebenda Bd. VI (1952), S. 1. 
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und fiir die Spannungen ergeben sich die Ausdiiicke 


a0 0 a0 
Ses — esos fe. Y = 
Ong? are en ee 2655 » On, =— 2GA0, Vig 2 © ae 


Erfiillen diese Werte nicht die vorgeschriebenen Randbedingungen, so mu8 durch Uberlage- 
rung eines weiteren Spannungszustandes mit T = 0 fiir die Befriedigung der Randwerte Sorge 
getragen werden. 


In dem vorliegenden Falle ist es méglich, die Auflésung der partiellen Differentialgleichung 
fiir © zu umgehen. Denn aus der Gleichung 


NO Nr und aus AT 2 784, also a ee) 
m— 1 Ox Ox 2A 
folgt 
0 Use! CLA Ree it net AT 
7549 Re Sy cary een SP bay Ty 
also 
0Ow ml Tes Reo Lo ie a Se 
pete a ne Goh peed ee 


Hieraus erhalt man wegen der Beziehung 
d 
4 (2) =— HI (e) 


(H} bedeutet die Hankelsche Funktion erster Art von der Ordnung 1) 


COS Wma eg Tl ee pr ys - OH (iAr) x 
cee 1 aes Ai HG (id1) + ED 
See dg a rer 7 _ &% my 
aap ae HG (ian : Hi (iAr)|. 
Fiir 
ao 20 
pele ines 
ergibt sich mit den vorstehenden Werten von 
eo icsee! _m+laq ans TTL (s 
Ox? und NO al teed uy (Ar) 


By ed ae HB EAD) TT Ge |b (id) — FH Gian 


dy? m—1 4k m—1 8k 


- m—I1 8k 
Endlich ergibt sich 


— mT) 446s] 5 Ht (éAr) apse (iar)]. 


pp — DE He Lay (an). 
Die Ausdriicke fiir die Spannungen lauten demnach 
ee Ce iH (iar) + (iar)], 
Cy, = 2c ee ts [iHg iar — 7 HY An), 
Osy = — 2G REA oie Ht (idr). 


Damit ist die Randbedingung da die Scheibe im Unendlichen spannungsfrei ist, bereits 
erfillt; denn die Hankelschen Funktionen verschwinden dort wie e—/". Die Spannungen 0,, 
hingegen kénnen durch Uberlagerung der Spannungen —o,, zum Verschwinden gebracht werden. 
Bei einer diinnen Scheibe werden hierdurch die Spannungen 0,,, Oyy und O,y nicht geandert. 
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Um einen Uberblick 


spannungen aus 


O12 


> 


itber die Spannungsverteilung zu erhalten, berechnen wir die Haupt- 


1 


1 
= 5} (Oxx Cyy) +) 


4 


(Gx — Oyy)” + Oxy 


— ~2G¢M@t1 %4 «ix [j Hi GAr) + At (éAr)] 


m—1 8k 


<\ 


Abb. 1. Hauptspannungstrajektorien. 


tg@P 

0,0000 
0,3287 
0,7500 
1,3333 
2,2913 


4,8990 


—A4,8990 
—2,2913 
—1,3333 
—0,7500 
—0,3287 


0,0000 


e= 0,1 

= 4,275 
0, = —7,071 
07 == 4,297 
On —7,107 
07 — 4,362 
0, = 7214 
oT 4,451 
Os = —7,363 
Oy = 4,540 
Os = —7,509 
= 4,631 
Os = 7,660 
O7 == 4,725 
On == —7,815 
0; = 4,820 
02, = 15913: 
= 4,918 
Og = —8,134 
0, = 5,012 
0, = —8,298 
0, = 5,119 
02 = —8,466 
= 5,196 
02 = —8,594 
0, = Does 
02 = —8,637 


und die Richtung von o, gegen die x-Richtung 


aus 
Ca One Orie COs RG 
Sees oxy y  rsing 3° 
Mit tg v= dy/dx=y’ erhalt man fiir die Haupt- 
spannungslinien die Differentialgleichungen 


y’ — peed bzw. yi = 


vy (po 


mit den Lésungen 


Sie yao Se 
a und x= He ih} 
also zwei Scharen confocaler Parabeln wie in 
Abb. 1 dargestellt. 

In der untenstehenden Tafel sind die Werte 


der Hauptspannungen und die Richtung von 0; 
gegen die x-Achse angegeben. Dabei wurde 


os ee aq(m+1)__4 
Ai Nia—e Ay = > Ar = Q und 8k(m—1) — 
gesetzt. 
0,5 1,0 1,5 2,0 
02826 0,0424 0,0090 0,0022 
—0,9965 —0,2396 —0,0698 —0,0219 
0,2898 0,0446 0,0112 0,0025 
—1,0220 —0,2518 —0,0862 —0,0242 
0,3124 0,0518 0,0122 0,0033 
—1,1010 —0,2926 —0,0942 —0,0326 
0,3452 0,0632 0,0165 0,0050 
— 207) —0,3574 Iai —0,0487 
0,3815 0,0772 0,0222 0,0074 
—1,3451 —0,4365 Al MAG —0,0726 
0,4216 0,0943 0,0300 0,0111 
—1,4866 —0,5331 —0,2316 —0,1083 
0,4660 0,1152 0,0405 0,0165 
—1,6430 —0,6521 —0,3127 —0,1616 
0,5150 0,1407 0,0547 0,0247 
—1,8158 —0,7954 —0,4221 —0,2410 
0.5692 0,1719 0,0738 0,0368 
—2,0068 07S —0,5698 —0,3595 
—0,6291 0,2099 0,0996 0,0549 
ee ee ae) —1,1866 —0,7691 —0,5363 
0,6952 02564 0,1345 0,0819 
—2,4510 —1,4492 —1,0382 —0,8002 
0,7493 0,2979 0,1483 0,1105 
—2,6419 —1,6838 —1,1452 —1,0801 
0,7698 0,3131 0,1815 0,122] 
—2,7088 Sle —1,4016 11937 


(Eingegangen am 9, Juli 1951.) 
Anschrift des Verfassers: Prof. Dr.-Ing. E. Melan, Wien IV, Technische Hochschule. 
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Die Stabilitatskriterien der Elastomechanik. 
Von H. Ziegler. 


1, Einleitung. In der Kinetik wird die Gleichgewichtslage eines Systems als stabil bezeichnet 
wenn sich dieses unter dem Einflu8 hinreichend kleiner, im iibrigen aber beliebiger Stérungen 
nur wenig von ihr entfernt (und damit eine kleine Schwingung um sie ausfihrt). Auf dieser 
Definition beruht das Verfahren der kleinen Bewegungen, das nur die Linearisierbarkeit der Be- 
wegungsdifferentialgleichungen voraussetzt und im Falle konstanter Beiwerte auf dem Wege 
itber den Exponentialansatz, die Wurzeln der charakteristischen Gleichung sowie die Routh- 
Hurwitz-Kriterien* die Stabilitat eines Systems zu beurteilen gestattet. 

In der Elastomechanik besteht das Stabilitatsproblem darin, die kritische Belastung zu be- 
stimmen, unter der eine urspriinglich stabile Gleichgewichtslage instabil wird. Bei seiner Lisung 
kann man sich unmittelbar auf die Definition der Stabilitat berufen und erhalt damit das voraus- 
setzungslos giiltige 

kinetische Stabilitatskriterium: Die kritische ist die kleinste Belastung, unter der 


eine geeignete Stérung zu einer nicht in der nachsten Umgebung der Gleichgewichtslage ver- 
laufenden Bewegung fiihrt. 


In der Praxis wird diesem kinetischen meist das bereits von L. Euler verwendete, in der Hand- 
habung wesentlich einfachere, aber nicht ohne weiteres mit ihm gleichwertige 
statische Stabilitatskriterium: Die kritische ist die kleinste Belastung, unter der 
neben der urspriinglichen (trivialen) erstmals eine weitere (nichttriviale) Gleichgewichts- 
lage existiert 
und in neuerer Zeit auch das 
energetische Stabilitatskriterium: Die kritische ist die kleinste Belastung, unter 
der die gesamte potentielle Energie des Kérpers nicht mehr positiv definit ist 
vorgezogen. 

Die Erfolge, die mit den zuletzt genannten Kriterien erzielt worden sind, haben die Frage 
nach ihrer Legitimitat in den Hintergrund gedrangt und zu einer rein statischen Konzeption des 
seinem Wesen nach kinetischen Stabilitatsproblems gefiihrt. Damit, d. h. mit der relativen Sel- 
tenheit nichtkonservativer Probleme, bei denen diese Konzeption zu Irrtiimern fiihrt, erklart 
sich, daB — um nur ein groBes Beispiel zu nennen — S. Timoshenko in seinem beriihmten Stan- 
dardwerk? ausschlieBlich — und ohne kinetische Begriindung — die beiden letzten Kriterien 
verwendet. 

A. Pfliiger hat in seinem Buche iiber Stabilitatsprobleme der Elastostatik® — freilich ohne 
Beweis und ohne die Konsequenzen aus dieser Bemerkung zu zichen* — darauf hingewiesen, daB 
bei nichtkonservativen Systemen nur das erste Kriterium in Frage kommt. Seither hat der Ver- 
fasser® nachgewiesen, daB beispielsweise das Knick- und das Drehzahlproblem der fliegenden, 
durch ein axiales Moment auf Torsion beanspruchten Welle nichtkonservativ ist und da hier 
das bisher verwendete statische Stabilitatskriterium versagt. 

Angesichts dieser Sachlage drangt sich die Frage auf, ob es nicht moglich sei, das statische 
oder das energetische Kriterium durch Umformulierung auch fiir nichtkonservative Systeme zu 
legitimieren. Zweck dieser Arbeit ist die Beantwortung dieser Frage, und zwar im negativen 
Sinne. 


2, Das kinetische Kriterium, Ein typisches (eindimensionales, nichtgyroskopisches® und kon- 
servatives) Stabilitatsproblem ist die Eulersche Knickaufgabe fiir den (beliebig, nach Abb. 1 


1 E. J. Routh, Advanced Rigid Dynamics, 6. Aufl., S. 228; London 1930; A. Hurwitz, Math. Ann, 46 
(1895), 5S. 273. * 

2 S. Timoshenko, Theory of Elastic Stability, New York u. London 1936. ; 

8 4, Pfliiger, Stabilitatsprobleme der Elastostatik, 5.67. Berlin, Gottingen, Heidelberg 1950. 

4 Trotz seiner Bemerkung auf $.67 behandelt er auf S. 217 ein nichtkonservatives (auch bei L. Collatz, 
Eigenwertaufgaben mit technischen Anwendungen, S. 41, Leipzig 1949 zitiertes) Problem mit dem 
statischen Kriterium. 

5 H. Ziegler, Z. angew. Math. Phys. 2 (1951), S. 265. \ 

6 Vel. H. Ziegler, Zum Begriff des konservativen Systems, Elem. der Math., erscheint demnachst. 
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etwa beidseitig gelenkig gelagerten) Stab mit der Masse (u(x) je Langeneinheit und der Biege- 


steifigkeit a(x). Seine Bewegungsenergie ist 
t 


Ta [ usta, 
0 


und die potentielle Energie setzt sich gemaf 


l 1 

1 D 

v-u—PM,, v=4{ aya, Wi= + | ytde (2.1) 
0 0 


aus der Forminderungsenergie und dem Potential der axialen Druckkraft zusammen. Ist ¥,( x) 
die zur Kigenkreisfrequenz o, gehérige k-te normierte Eigenfunktion des unbelasteten, quer- 
schwingenden Stabes, so kann man die Schwingungen unter der Druckkraft Pin der Form 


y(t) == 3 4 (0) pe(2) 


darstellen und die q, als Lage-, fiir P== 0 als Normalkoordinaten deuten. In ihnen driicken sich 
die drei Bestandteile der Gesamtenergie in der Gestalt 


U 
lw. LS LSS ‘Mh 
bere U=— > ota. Wy= > D>, 9% | ivi ds 


1 I cree 
oe (ii — ok ai) +P Sam | vivian. 
k=1 t,k=1 0 
und die Lagrangeschen Bezichungen 
LAST PE MART pe 4 /aL\ at 
ruckter, beidsel 1g gelen <1 ge ager er abd. AS a fai 
EP (Fac) 3a = ° al bein aes) ek 24) 
fithren auf die Bewegungsdifferentialgleichungen 
ie 1 
dit ohge— P >’ a | gigide=o (kea=) 2 rele (2.3) 
ae | 
Im Falle des prismatischen, homogenen und beidscitig gelenkig gelagerten Stabes ist 
oy 1B tae 08 Pe ae aes 
ee ian MC aE? w(x) =|, sin j 5 


und da hicr neben den Eigenfunktionen g, auch ihre Ableitungen gy im Intervall 0. ..1 ortho- 
gonal sind, zerfallt das simultane System (2.3) in die unabhangigen Differentialgleichungen 


i. Ls ai (kone 

qi -| (“i Fi uP P) a= 0 (eS 2) ee) ie (2.4) 
Die Eigenkreisfrequenzen unter der Belastung P sind mithin durch 

, LR apt ANE aed 
og =F i —— (aah 2) 
gegeben, und die Forderung, da sie samtlich reell seien, fithrt auf die Eulersche Knicklast 
me 0 
PS Tee (2.5) 


Das kinetische Stabilitatskriterium gilt, da es unmittelbar auf der Definition der Stabilitat 
beruht, voraussetzungslos. In der Tat ist leicht einzusehen, da® sich das hier verwendete Verfahren 
ohne weiteres auf mehrdimensionale Aufgaben und dadurch, da® man die rechten Seiten der 
Lagrangeschen Gleichungen (2.2) durch verallgemeinerte Krifte Q, erganzt, auch auf gyrosko- 
pische und nichtkonservative Probleme iibertragen la8t. Da aber das simultane System (2.3) nur 
ausnahmsweise zerfallt, ist die Methode im allgemeinen umstandlich. 


XX. Band 1952. Ziegler: Die Stabilitatskriterien der Elastomechanik, 51 


3. Das energetische Kriterium, Werden die Lagekoordinaten eines konservativen elastischen 
Systems von der trivialen Gleichgewichtslage (d. h. derjenigen des unbelasteten Systems) aus 
gemessen und die potentielle Energie in dieser Lage Null gesetzt, so gelten die Beziehungen 


V(0,...)= De 
( 2 ) 0 ? 0q;, 
Nach dem Satz von der Erhaltung der Energie nimmt die Bewegungsenergie T mit zunchmender 
Entfernung von der Gleichgewichtslage mit Sicherheit oder bestenfalls unter Umstianden ab, je 
nachdem V(q,, gp, - . -) positiv definit ist oder nicht. Bei konservativen Systemen kann also das 
kinetische durch das energetische Kriterium ersetzt werden. 
Oft 1aBt sich die potentielle Energie in der Umgebung der Gleichgewichtslage durch die ab- 
gebrochene Potenzreihe 


(0, ...)=0 (k= 1,2,...). (3.1) 


Sareea lw ey 

V(qu» Ge» ---)= V(0,...)+ ye aU eke mrap ay rr (3.2) 
darstellen, und da hier zufolge (3.1) neben dem ersten Term auch der zweite, namlich die erste 
Variation von V verschwindet, kann das energetische Kriterium auch dahin formuliert werden, 
daf} der dritte Term, der — vom Faktor 1/2 abgesehen — als zweite Variation von V bezeichnet 

wird, positiv definit sein mu, wenn das System stabil sein soll. 
Im Beispiel von Abschn. 2 ist die potentielle Energie V= U -— PW, nur fiir geniigend kleine 
Werte von P positiv definit. Beim prismatischen, homogenen und beidseitig gelenkig gelagerten 
Stab folgt aus (2.4), daB die q, auch fiir den belasteten Stab Normalkoordinaten sind, und die 


Bedingung, daB 
1 > ié 
V= oy & OK qi: 


positiv definit, die o,2 mithin samtlich gréBer als Null seien, fiihrt wieder auf die Eulersche 
Knicklast (2.5). 

Der Giiltigkeitsbereich des energetischen Kriteriums ist offensichtlich insofern enger als der- 
jenige des kinetischen, als es die Existenz einer stationaren, eindeutigen potentiellen Energie 
voraussetzt. Es bleibt indessen denkbar, daf es sich fiir nichtkonservative Systeme durch eine 
andere Uberlegung statischer Natur, z. B. eine Arbeitsbetrachtung ersetzen 14Bt. 

4, Das statische Kriterium. Lat die potenticlle Energie eines konservativen Systems die 
abgebrochene Entwicklung (3.2) zu, so kann sie, indem man die Darstellung (2.1) durch Einfiih- 
rung eines fiir die Belastung charakteristischen Parameters A verallgemeinert und (3.1) sowie 
(3.2) beizieht, in einer der Formen 


f°) 


; fl 
Va EN a ay a A qi Us ay) 
dargestellt werden; zudem ist sie mit der Energie U der inneren Krafte fiir geniigend kleine Werte 
von A positiv definit. Wenn man jetzt statt T und V die quadratischen Formen U und W, auf 
Hauptachsen transformiert, erhalt man an Stelle der q, neue Normalkoordinaten p,, in denen sich 


U und W, durch . 
ee a es a ea 
Ua paes Mo= > 2 uph (4.2) 


ausdriicken. Es ist daher 


cee : OV dt 

amie canne ag, 7 (LA) Pe (k=1,2,...). 

=] 
Solange/ geniigend klein ist, sind die Klammerausdriicke in diesen Beziehungen samtlich positiy, 
und da die Relationen 
a 1 — Ae’ py 20 (Ree 1 not) (4.3) 
04), 

die Gleichgewichtsbedingungen des Systems darstellen, existiert keine nichttriviale Gleich- 
gewichtslage. Wird andererseits mit zunehmendem / erstmals cine der Klammern Null, so ver- 
liert V seinen positiv definiten Charakter ; gleichzeitig 1aBt aber (4.3) eine nichttriviale Lésung zu. 

A* 
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Beim prismatischen, homogenen und beidseitig gelenkig gelagerten Stab von Abschn. 2 
fallen die Normalkoordinaten p, mit den q, zusammen. Die nichttrivialen Gleichgewichtslagen 
sind durch das Verschwinden der Klammern in (2.4) gekennzeichnet, und durch Nullsetzen der 
kleinsten erhalt man wieder die Eulersche Knicklast (2.5). 

Diesen Uberlegungen zufolge gilt das statische Stabilitatskriterium fiir konservative Systeme, 
deren potentielle Energie sich gema8 (4.1) baw. (4.2) darstellen 1a48t. Es bleibt aber denkbar, daB 
sich sein Giltigkeitsbereich — unter Umstianden sogar auf nichtkonservative Systeme — er- 


weitern laBt. 


5, Nichtkonservative Probleme. Schwingungsfahige Systeme sind schon mit Riicksicht auf 
die stets vorhandene, vom Bewegungszustand abhangige Dampfung nicht konservativ. Imallgemei- 
nen kommt man bei Vernachlassigung der Dampfung auf ein konservatives 
y System zuriick, das sich mit jedem der drei Kriterien behandeln Ja4Bt. Daneben 
kommt es indessen — und dieser Fall soll hier allein untersucht werden — 
vor, daB® ein elastisches System bei der Belastung seinen konservativen 
Charakter deshalb verliert, weil sich die vom Bewegungszustand unabhangigen 
Lasten nicht von cinem stationaren, eindeutigen Potential ableiten lassen. 
Hierher gehért der von A. Pfliiger’ angefiihrte Stab, der eine stetige, in die 
Tangente der elastischen Linie fallende, bei der Deformation in dieser ver- 
bleibende Belastung tragt, ferner die vom Verfasser® behandelte, am freien 
Ende durch einen axialen Momentvektor belastete fliegende Welle. Der 
Abb. 2. Nichtkonser- Grund fiir dieses anormale Verhalten liegt darin, da weder die mit einem 
T apeetial cearanes Kérper starr verbundene (mithin an seiner Drehung tcilnehmende) Kraft mit 
Ben peumeel I gieh ss konstantem Betrag noch der konstante Momentvektor konservativ ist*. 
ee Als Musterbeispiel fiir diese Problemklasse, die zwar nicht sehr umfang- 
reich, aber praktisch doch von Bedeutung ist, kann der einseitig eingespannte, homogene und 
prismatische, am freien Ende durch eine tangentiale Einzelkraft P belastete Stab (Abb. 2) be- 
trachtet werden. Die bei der Verschiebung in 
die (triviale) Gleichgewichtslage von P geleistete 
Arbeit hangt von der Art ab, wie diese Verschie- 
bung vorgenommen wird; die Last P ist also nicht 
konservativ. Das energetische Kriterium versagt 
in der in Abschn.1 gegebenen Form, und da man 
miihelos zeigen kann, daB keine nichttriviale Gleich- 


1 


gewichtslage existiert, mite mit dem statischen 
Kriterium auf eine unbegrenzte Stabilitat des 
Stabes geschlossen werden. 

Die in Abschn.3 und 4 offen gebliebene Frage, 
ob sich das statische oder das energetische Kri- 
terium so modifizieren lasse, daB es auch auf nicht- 
konservative Systeme anwendbar wird, 148t sich 

‘~_ schon an einfachen Beispiclen beantworten, so etwa 
ae an einer Gelenkstabkette, die nach K. Marguerre* 
Abb. 3. Doppelpendel als zweigliedriges Modell des 


RSet Go kes eae ALE o als vereinfachtes Modell des elastischen Stabes 
(Abb. 2) dienen kann. 


6, Ein Modell, Abb. 3 zeigt ein Doppelpendel mit den (als klein vorausgesetzten) Drehwin- 
keln 1, @g, bestehend aus zwei starren Staben der Lange /, deren Massen m,, m, in den Abstanden 
a,, a, von den Gelenken konzentriert gedacht sind. In m,, m, greifen zwei vertikale Krafte G,, G, 
an, die als Gewichte interpretiert werden kénnen, am freien Ende des zweiten Stabes die axiale 
Kraft P und in den Gelenken die Riickstellmomente cq,, ¢(p~, — 7) sowie die Dampfungsmomente 
bg, BG. — 9) 


Die gesamte Bewegungsenergie ist bis auf GréBen zweiter Ordnung genau 


Wa 


Oe F,)#0(Go-P) 


1 ; ae ; 
fe => [(m, af + m, I?) pi + 2 m, Lay Py Pz + mM, a2 GY] ’ 
2 


1 A. Pfliiger, a.a.O., S. 41. 

2 H. Ziegler, Z. angew. Math. Phys. 2 (1951), S. 279. 

3 Vel. H. Ziegler, Zum Begriff des konservativen Systems, a. a. O. 

* K. Marguerre, Neuere Festigkeitsprobleme des Ingenieurs, S. 211. Berlin, Géttingen, Heidelberg 1950. 
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die potentielle Energie der Krafte G,, G, sowie der Riickstellmomente 


1 
yee 7 [(G,a, + G1 2 c) pi — 2¢9, %.+ (G, a, + c)@3] » 


und die verallgemeinerten nichtkonservativen Krafte (die von P und den Dampfungsmomenten 
herrithren) berechnen sich zu 


Q, = Pl(p, — $2) — b (2g, — Gr), Q2 = bp, — 2). 
Die Lagrangeschen Relationen 

d /o0L OL 

e (55.) ape (Le TS V fe shay 2) 


fiihren auf die Bewegungsdifferentialgleichungen 
(m, aj + my 2?) Gy + m, Lay Gp, + b(2 Gy — 2) + [G, a, 4 (G,— P)l+ 2c] p,+ (Pl—)y,=0, 
My La, 1 - m, a; P, — b (py — H2) — €9 + (6, 4,-+ ¢)p,= 0, 


‘ 


und der Exponentialansatz 


GP, = A, e™* (k = 1, 2) (6.1) 
auf die charakteristische Gleichung 
PoM + P14? + Po? + psd + Py= 9, (6.2) 


deren Beiwerte durch 


Po = MM, a} a, 
Pi = [m, a} + m, (I? + 2 la, + 2 43)) b, | 
P2 = (m, a} -- m, I?) (G,a,-+ ¢) + m, a2 [G,a,-+ (G,— P)1-+ 2c] — m, la, (P1— 2c)-+ b?, (6.3) 
Ps = [G,4,+ 6, (14 24,) + 2¢] 6, | 


Ps = [G,4,-+ (G, — P)l-+ 2c] (G,a,+ c) + (Pl—c)c 
gegeben sind. 


Das Doppelpendel ist solange stabil, als simtliche Wurzeln A,, ...,A, der charakteristischen 
Gleichung nichtpositive Realteile aufweisen. Dann und nur dann bleibt namlich die allgemeinste 
Bewegung, die sich aus vier Lésungen der Form (6.1) zusammensetzt, beschrankt. 

Sind — wie dies im folgenden mehrfach zutreffen wird — je zwei Wurzeln entgegengesetzt 
gleich (etwa A, = —/,,/, = —A,), dann ist der Schwinger solange stabil, als die Werte 1} = 43, 
22 =} negativ oder Null sind. 

Verschwindet eine der Wurzeln, so existiert nach (6.1) eine Lésung mit konstanten, nicht 
gleichzeitig verschwindenden Winkeln ¢,, ¢, d.h. eine nichttriviale Gleichgewichtslage. Um- 
gekehrt verlangt die Existenz einer nichttrivialen Gleichgewichtslage das Verschwinden einer 
dieser Wurzeln. 

7. Folgerungen, a) Setzt man b= 0, m, = 2m, m,=m, a4, =4,=—1,6,= G, = 0, so 1aBt 
sich das Pendel der Abb.3 als zweigliedriges Modell des (um gedrehten) Stabes von Abb.2 auf- 
fassen, wobei die Masse auf die,Knoten verteilt und das Eigengewicht samt der Dampfung ver- 
nachlassigt ist. Die charakteristische Gleichung geht dann mit 


Do 2 m7 tt, Po mi (ic— 2 Pl), Pee (7.1) 
in die biquadratische Gleichung 


Po + P27? + Py = 9 (7.2) 
mit der Diskriminante 
a= pe ope pean (41 c2-— 28 PL + 4D?) (7.3) 
iiber. 

Fir P= 0 sind die Wurzeln 72 = 22, 22 =/72 der charakteristischen Gleichung, da neben der 
Diskriminante (7.3) alle Beiwerte (7.1) positiv sind, kleiner als Null (Abb. 4). Die allgemeinste 
Bewegung ist mithin beschrankt und setzt sich in bekannter Weise aus den beiden Normal- 
schwingungen zusammen. Lat man P anwachsen, so nimmt A ab und wechselt fiir P= 
(7/2 — y2) c/l= P,, das Vorzeichen. Die Wurzeln 72, 72 nahern sich daher mit zunehmendem P, 
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allen fiir P — P, zusammen, und da sie fiir gréBere Werte von P konjugiert komplex werden, ist 


P, = (7/2 — y2) 5 = 2,086 — (7.4) 


die kritische Belastung. 


Da nach (7.1) p, >0 ist, sind die Wurzeln von (7.2) fiir beliebige Werte von P von Null ver- 
schieden; es gibt mithin — wie man iibrigens auch durch Betrachtung des Kraftespiels in Abb. 3 


bestatigt — keine nichttriviale Gleichgewichtslage. Damit folgt schon aus diesem einfachen 
Beispiel 
Satz 1: Bei nichtkonservativen Systemen kann das statische Stabilitatskriterium ver- 
sagen. 


Man koénnte aber noch vermuten, daf es nur bei Abwesenheit einer nichttrivialen Gleichgewichts- 
lage versage. 
b) Modifiziert man das Beispiel durch Beriicksichtigung des Eigengewichtes, indem man statt 
G,= G,= 0 jetzt 6, = 26, 6, = © setzt, so erhalt man statt (7.1) 
Po ean ls P2= mP(7ce— 2 Pl-+ 661), | 


P= + 5e6l— PGP +3CP f (7.5) 


und an Stelle von (7.3) 
A = m?I4 [41 c2?— 4c(7 P— 11G6)1+ 4 P? 22 — 16 PGP 12671). (7.6) 


Fix P= 0 sind die Wurzeln J? = /2,/2 = 72 von (7.2) wieder negativ, das Pendel also erwar- 
tungsgema$ stabil. LaBt man P anwachsen, so geht die Stabilitat verloren, sobald (Abb. 4) ent- 
weder A* = 0 oder /? = 7? wird, denn dann verlaft mindestens /? die negative reelle Achse. Die 
kritische Belastung bestimmt sich daher aus der scharferen unter den beiden Forderungen 


Ps= 9, A= 0, Cit) 


von denen die erste gleichzeitig die Bedingung fiir die Existenz einer nichttrivialen Gleich- 
gewichtslage ist. Fithrt man (7.7) vermittelst (7.5) und (7.6) aus, so erhalt man auf Grund der 
statischen Stabilitatsbedingung allein die kritische Belastung 


: Lee meet es 
| [A] ee 2 G? 
— =. --+3 
LO a ro ae Ao 
cp aS ha der kinetischen zufolge aber zusatzlich 
\) 
Abb. 4, Verschiebung der Wurzeln der ah G / me z 
charakteristischen Gleichung mit —- | | | 2 
zunehmender Belastung. P, FF 26 V 22 3 1 G | G?. 


Da P;, > O ist, gibt es eine Belastung, fiir die eine nichttriviale Gleichgewichtslage existiert. 
Da indessen P;, << P, ist und P, fiir G— 0 iiber alle Grenzen wiachst, ist aber P, und nicht die vom 
statischen Kriterium gelieferte, unter Umstanden sehr viel gréfere Belastung P, die kritische. 
Ks gilt mithin : 


Satz 2: Das statische Stabilitatskriterium kann auch bei Vorhandensein einer nicht- 
trivialen Gleichgewichtslage versagen (und unter Umstanden eine viel zu hohe kritische Be- 
lastung liefern). 


Mit diesem Satze ist die Unbrauchbarkeit des statischen Kriteriums nach Abschn. 1 fiir nicht- 
konservative Systeme nachgewiesen; die Frage nach der Méglichkeit, P, auf anderem, rein sta- 


tischem Wege — z. B. durch eine Betrachtung der Arbeiten — zu bestimmen, bleibt indessen noch 
offen. 


c) Vernachlassigt man mit b — 0, CG, = G, = 0 die Dampfung und das Eigengewicht, so er- 
halt man mit m, = m, = m, a,= a,= 1/2 ein zweigliedriges Modell des (um z gedrehten) Stabes 
von Abb. 2, das dem unter (a) besprochenen gleichwertig ist und sich nur darin von ihm unter- 
scheidet, dafs jetzt die Masse nicht in den Knoten, sondern in den Mittelpunkten der Stibe kon- 
zentriert gedacht ist. Die charakteristische Gleichung ist auch hier (7.2), wobei aber (7.1) und 
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(7.3) durch 


1 
Po= gmbh, Po= mi? (Ile — 3 Po), ee 


A = = ml! (39 c? — 22 ¢ PI 3 P22) 


zu ersetzen sind. 


Da wieder p,> 0 ist, erhalt man die kritische Belastung durch Nullsetzen von J, und zwar 
ergibt sich mit 


Cc 
B= 3 


ein wesentlich héherer Wert. als er mit (7.4) im Falle (a) erhalten wurde. Es gilt also 


Satz 3: Bei nichtkonservativen Systemen hangt die kritische Belastung bei sonst gleichen 
Verhaltnissen (unter Umstanden sehr stark) von der Masseverteilung ab. 


Bei der Vereinfachung solcher Systeme ist mithin Vorsicht geboten. 
Da die Masse in statische Untersuchungen nicht eingeht, folet aus Satz 3 


Satz 4: Nichtkonservative Stabilitatsprobleme lassen sich auf statischem Wege nicht 
lésen. 


Damit ist gezeigt, da® sich auch das energetische Kriterium nach Abschn. 1 fiir nichtkonser- 
vative Systeme nicht legitimisieren la8t. Mit der Feststellung, da8 nur das kinetische Kriterium 
anwendbar bleibt, riickt aber die Frage nach dem Einfluf der Dampfung in den Vordergrund. 


d) Mit G, = G,= 0, m,=2m, m,=m, a, = a,=/ erhalt man ein zweigliedriges Modell 
des Stabes von Abb. 2, bei dem auch die innere Dampfung beriicksichtigt ist. Die charakteri- 
stische Gleichung hat die Form (6.2), wobei die Koeffizienten durch 


Poem, Daeg ouea 
Po mE (ice—2 Pl) bh, p= 2eb., (7.8) 
Poe 


gegeben sind. Da die Gleichung (6.2) im Gegensatz zu (7.2) nicht 
mehr quadratisch ist, empfiehlt sich die Verwendung der Routh- 
Hurwitz-Kriterien!, und diesen zufolge ist das System jedenfalls 
stabil, solange die Ausdriicke 

__ Po Ps 


=, PiPs 
Por Pir Pa a Sa 


Ps 


das gleiche Vorzeichen besitzen. Setzt man hier die Werte (7.8) ein, 
so kommt man auf die Stabilitatsbedingungen 


41 ¢ 1 b 
EE ee Lael ee 
ks SO 3 e 28 1 ae D sale ° \ 
Abb. 5. Stabilitatsbereich des 
il b2 Modells bei Beriicksichtigung der 
2. P < 45 Ee — = 6 4. ¢? = 0 > Dampfung (b Dampfung, 
14 1 DD pl P Belastung), 


von denen die vierte auf alle Falle erfiillt ist. Aus der ersten schlieBt man, daB das System fir 
negative Dampfungen (Abb. 5) instabil ist; aus der dritten, die scharfer als die zweite ist, ergibt 
sich fiir positive Dampfung die kritische Belastung 


tee 


ee r 2 mB’ 


und diese geht fiir b +0 nicht gegen (7.4), sondern gegen 
Pe= 1,464 —. (7.9) 


1 A. a. O. zitiert. 
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Mathematisch laBt sich die Diskrepanz zwischen P, und Pf damit verstehen, dafs die Routh- 
Hurwitz-Kriterien die Bedingungen fiir Abklingen der allgemeinsten Liésung darstellen, die 
stationaéren Schwingungen unter der Belastung Pf <P<P,, im Sinne dieser Kriterien also be- 
reits alsinstabil gelten kénnen. Dem entspricht die mechanische Feststellung, daf stets mit einer 
— wenn auch noch so geringen — Dampfung zu rechnen ist, und da diese im Intervall Pf < P< P, 
abilisierend wirkt, mu vom physikalischen Standpunkte aus auch fiir das praktisch ungedampfte 
System P¥ als kritische Belastung angesehen werden. Man hat also 


Satz 5: Bei nichtkonservativen Systemen kann die Dampfung labilisierend wirken, 


sowie schlieBlich 


Satz 6: Bei nichtkonservativen Systemen kann die geringste Dampfung die kritische Be- 
lastung erheblich modifizieren. 


Systeme der hier betrachteten Art lassen sich also nur mit dem kinetischen Stabilitatskri- 
terium behandeln. Dieses ist in der Anwendung wesentlich umstandlicher als die beiden anderen 
und wird mit der als notwendig erkannten Beriicksichtigung der Dampfung noch schwerfalliger. 
Die sich hieraus ergebenden Komplikationen diirften sich indessen kaum umgehen lassen. Sie 
liegen wohl in der Natur des Problems, und zwar in der Tatsache begriindet, daB die Energie- 
aufnahme des nichtkonservativen Systems durch geringe Variationen im Bewegungsablauf emp- 
findlich beeinflu8Bt werden kann. 


(Eingegangen am 28. Juli 1951.) 
Anschrift des Verfassers: Professor Dr. H. Ziegler, Riischlikon bei Ziirich, Weierweg 6. 
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Bewegung des langgestreckten Rotationskérpers 
in einer zur Langsachse geneigten Richtung. 


Von W. Miiller. 


1, Die grundlegenden Apestiee fiir die Querbewegung, In einer friiheren Mitteilung! ist die 
Darstellung der Langsbewegung eines Rotationskérpers mittels einer stetigen Quellsenkenver- 
teilung f(t) langs der auf der Achse gelegenen Grundstrecke (— 1 <t <-+ 1) gegeben worden. 
Die Querbewegung in Richtung senkrecht zur Achse 148t sich bekanntlich aus einer Verteilung 
von Dipolen gewinnen, die nach der Bewegungsrichtung, etwa der y-Achse orientiert sind. Mit 
der Verteilungsfunktion F(t), die langs derselben Grundstrecke definiert sein soll, erhalt man 
das Potential 

1 


rm 
oS — Fre SN (1) 
—~] 


wo r? = (x — t)?-+ 7? zu setzen ist, py den Polarwinkel gegen die y-Achse in der Parallelkreis- 
ebene und # = arc sin (y/r) den Winkel zwischen der Lingsachse und der Verbindungslinie des 
laufenden Dipolpunktes mit dem Aufpunkt P bezeichnet (Abb. 1). 

Die Lésung fithrt durch Umformung 
auf tesserale oder zugeordnete Kugel- 
funktionen P}(u) und Qj(¢) mit dem 
Index 1, wobei die Argumente jz, € mit 
den friither definierten Ellipsoidkoordi- 
naten identisch sind. Man erhalt dann 
die Darstellung 


D =F cosy SA, Pru Q(0) 2) 
oder mit Hilfe der bekannten Beziehun- 
gen P,(w) ad jl— 1? P; (M4); Q, (6) 

= — y—103() 
P= — SB, PW) O00) (3) 


Abb. 1 Grundachsen und Koordinaten, 


Man wird auch hier zu der Bedeutung und zur Bestimmung der Koeffizienten /, gelangen, wenn 
man den Grenziibergang zu Punkten durchfiihrt, die in unmittelbarer Nahe der mit Dipolen 
besetzten Grundstrecke liegen. Es ergibt sich dabei 


= 5 : B co = 
lim Q,(6) = 1, foe tore oa Ty, lim B= —5, TT by Pt). 


Der letzte Ausdruck stimmt aber formal iiberein mit dem ebenen Potential 


BF(t)cosp____—s mcosep 
277 27 


od — 


einer senkrecht zur (y,z)-Ebene gelegenen und nach der y-Richtung orientierten unendlich 
langen Dipolstrecke mit dem Moment m. Wir haben daher 


Fi > 6,2, (0; (4) 
oder die GréBen f, stimmen iiberein mit den Koeffizienten der nach den iberstrichenen 


Kugelfunktionen P, entwickelten (reduzierten) Verteilungsfunktion. Zu ihrer Bestimmung be- 


1 Ing.-Archiv 19 (1951) S. 282. 
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nutzen wir die Orthogonalitatsbeziehungen fiir P;, 


S| +1 
2n(n +1) 
Pi(t) Pi) dt=0 (n+), | [Ps(e)P de ED. (5) 
peG | =] 
vee ni avers les ycunles Dy P, 
Dann ergibt sich aus (3) durch Multiplikation mit P; = aa fase und Integration 
—t car 
41 i #3 
/ = Ph (t) = Ey wae 1 ; 1 2n(n + 1) 
| F(t)ePA) dt [#2 Da0- PPh (0) dt Pose Ph(t) dt= By 
=i =i =i 
oder mit geindertem Index 
an 
_ eer i i Pe 
Baas) [Fo P;(t) dt. ( ) 


=i 


Man kann daher das Potential (2) auch in der geschlossenen Integralform schreiben 


fl 
cop > morn (u) Q2(C) | F(t) Pi(t) dt. (7) 


Aus gewissen formalen Griinden werden wir im Folgenden haufig statt der Koeffizienten £, die 
eréBeren Koeffizienten /; = v(v-+ 1) 6, verwenden. 

Im allgemeinen beschranken wir uns ferner auf solche Funktionen F(t), die in den End- 
punkten t= + 1 der Grundstrecke Nullstellen haben. Dann kénnen wir setzen 


P(t) = (1— #) g(¥), 
wo g(t) eine ganze rationale Funktion bedeuten soll, die keine Nullstellen zwischen —+- | besitzt, 
so daf F(t) in dem Definitionsbereich iiberall positive Werte annimmt. Mit Einfithrung von g 
1aBt sich (6) auch in der Form 


eee oe) 
20(v 1). g(t) P, (t) dt (8) 


—l 


bp 


schreiben. Entsprechend ]a8t sich auch die Formel (7) auf die Funktion g(t) reduzieren. Man 
hat nur statt F(t) die Funktion g(t) und an Stelle von P/ die iiberstrichene Funktion P, zu 
schreiben. 

Wenn wir g(t) etwa in der Form (t-+ 1)™ ansetzen, so erhalten wir bei Benutzung der Funk- 
tionen 


— = — 3 Se i ae nwa 
YR We! pavers 1 RE 55 Pe nO aah) late Dames eel es (9) 
folgende Tabelle 1 fiir die Koeffizienten 6 in den Fallen m= 0 bis m= 5. 
Tabelle 1. Die Koeffizienten B. Damit ergeben sich z. B. fiir die zusammen- 
=e aL ALA Lm [oa sesetate Form 
g(t) = valt + 1) + ye(t+ 1)? + ya(t+ 18 
0 1 + Ya(t +> 1)? 
1 1 = folgende Koeffizienten: 
6 8 16 
es) cme ie ey i ee es 
5 3 15 i 
8 8 2 9 To , 2 ow , 40 
3 ; : : 7 — = Bs Sl eg) 2c er) Sele og) ae 
eon ees: 18 8 8 pagans Pipnatal te ree 
Tele) 08 Gael e lasts eae ERC tee ce sey 
5 24 200 16 48 8 8 pr 2 8 
TES | Oar nes os enemas ena ee 


Wir werden im folgenden auf diese Darstellung noch zuriickkommen. 
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2, Berechnung der Geschwindigkeiten, Nach den in der vorigen Mitteilung angegebenen 
Formeln erhalt man bei Einfithrung der Koeffizienten / fiir die Komponenten der Absolut- 
geschwindigkeit der Fliissigkeit folgende Ausdriicke : 


=< B cos p ; B cos p 
On > »G, (Cu = > é ( 
’ 22 (C2? — pn?) B; G, (Cu), Vy 2 7 9 (62 — w?) Bb; H, (Cu). (11) 
Dabei wird, wie die Rechnung ergibt, 


Gm Y= CPE) — HPA =LEP Gn Cal | 
i Aan ak (12a) 
ord wo, Keres Q, Se PP, Q,44] ° | 
0 1 dQ; mM 
H, = 565 H (6 PH) SE —p Sao 
=— B fee) P, Obs STG) PEON aeRISE + I (C2 — pu?) P, Q, (12b) 
Nes Bop — — = = | 
aie PPO Ow) PO, ta 1) Peel. | 
Mit den Abkiirzungen 
pr PaO) PnQn Tn Pn (u) Qn(£) 
*. n(n + 1) n(n+ 1) 7? ne” "—~ (n—I) n(n+ 1)(n+ 2) 


erhalt man ferner die Rekursionsformeln = 
G, = Gr» — (2m — 1) n? (0? — w*) T,_,, 
eet 78) mea ie Fag on Tym (ne Dus | (13) 
ed ee) (Ce ie ol) 9 Se i SE 


Im Besonderen wird 


G,=—pw(l—p), G,=(1 “)| Cn @ 1) 0; | 


: (14) 
(c? 1) Q); ee eh ye 


1 — 

5 (212) 0. | 
3. Beispiel des Ellipsoids, Wie bekannt, 1aBt sich die Langsbewegung des langgestreckten 

Ellipsoids, dem ein konstanter Wert von € = Cy zugeordnet ist, aus der Quellverteilung f(t) = t 

ableiten. Es 1a48t sich nun zeigen, daB die Querbewegung durch eine Dipolverteilung 


F(i)=1—#®, g=1 ul 
erzeugt werden kann (Abb. 2). 


Wir erhalten hier nur einen 
Koeffizienten 


1— il 
H, =" |t—1)— 5 


18 


+1 
—— | (hae) di 1 
= 
unserer Entwicklung und da- 
her das Potential 


= — Seen = G, (0). (15) 


Die Absolutgeschwindigkeiten 


werden 


Abb. 2. Die Verteilungsfunktionen f(t) und F(t) fiir die R, (Ellipsoid). 


Be Aa eB 2c) —(C?—p*) 1] (16) 
3) Pere ae Geel 2 0(¢? — wu?) (C2? — 1) 
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Wenn man zum Ausdruck bringt, daB im Schnitt mit der Oberflache (¢ = Co) mit der Meri- 
dianebene y = 0 das Verhiiltnis der Relativgeschwindigkeiten (v, — v,)/v, der aus der Betrach- 
tung der Lingsbewegung hervorgehenden Tangentenneigung 


a MWe est 17) 
ct ea Go Vale ( 
gleich sein mu, so erhalt man durch eine einfache Rechnung 
= ag (18) 
eee Alo) 
Co 


Da die Lingsbewegung des Ellipsoids bekannt ist, so erhalt man durch Uberlagerung mit der 
Querbewegung auch die schiefe Bewegung, d.h. die Bewegung in einer Richtung y gegen die 
Langsachse, fiir welche die Bezichung tg y = v,/v, gilt. 


4, Der allgemeine Fall der zusammengesetzten Bewegung und das Hauptproblem, Wir haben 
in der ersten Mitteilung gezeigt, wie man aus einer stetigen Quellsenkenverteilung langs der 
Grundstrecke bei Zutreffen der SchlieBungsbedingung die Kontur des zugehérigen Rotations- 
kérpers und alle Daten der Langsbewegung dadurch ermitteln kann, daf} man die Koeffizienten 
a, aus der Entwicklung der Verteilungsfunktion f(t) nach einfachen Kugelfunktionen P, ab- 
leitet. Insbesondere sind mit Hilfe dieser Koeffizienten unmittelbar die Potentialfunktion @, 
die Stromfunktion Y und die Gleichung der Meridiankontur zu gewinnen. Die Querbewegung, 
die sich aus einer axialen Verteilung von Dipolen ableiten ]a8t, bringt nun aber eine nicht 
mehr rotationssymmetrische Strémung hervor. Da insbesondere keine Stromfunktion existiert 
wie im Falle der Langsbewegung, so ist es zunadchst nicht méglich, bei gegebener Verteilungs- 
funktion F(t) die Kontur in expliziter Form darzustellen. Man muf also nach einer indirekten 
Methode suchen, um die Meridiankontur oder die Quellsenkenverteilung aufzufinden. Der ein- 
zige Fall, bei dem diese Beziehung zwischen Lings- und Querbewegung in geschlossener Form 
angegeben werden kann, ist der Fall des Ellipsoids, wo die normierten Funktionen die einfache 
Gestalt f(t) = t, F(t) = 1— # haben. In diesem Falle ist aber die Funktion f(t) der negativen 
Ableitung der Funktion F(t) nach t proportional: 


S(t) ee eb (7): (19) 


Unter dem Hauptproblem der Hydrodynamik des Rotationskérpers werde ich die Fragen 
zusammenfassen nach der Auffindung der erzeugenden Dipolverteilung eines durch die Quell- 
funktion f(t) bestimmten Rotationskérpers oder nach der Ermittlung des Zusammenhangs beider 
Verteilungsfunktionen f(t) und F(t), wenn die Meridiankontur des Kérpers entweder graphisch 
oder rechnerisch, also etwa durch die Koordinaten der Randpunkte gegeben ist. Nach Beant- 
wortung dieser Fragen kann auch der ganze Verlauf der raumlichen Strémung, bzw. die ent- 
sprechende Druckverteilung an der Oberflache als bekannt gelten, die entsteht, wenn der Dreh- 
kérper unter irgendeinem Winkel gegen die Langsachse translatorisch in der reibungslosen 
Flissigkeit sich bewegt. Obwohl die mathematisch strenge Lisung des Hauptproblems bis jetzt 
nicht aufgefunden worden ist, so sind doch mehrere, fiir die Praxis ausreichende Naherungs- 
lésungen angebbar, die der vollstandigen Liésung beliebig nahe kommen. 


5. Erste Naiherungslésung. Fiir den Fall, daB der Meridianschnitt des Koérpers nicht allzu 
stark von der Ellipse abweicht, kann man in erster Annaherung annehmen, daf die Funktion f(t) 
der negativen Ableitung der Funktion F(t) proportional ist. Im allgemeinen ist damit in der 
Tat die charakteristische Beziehung zwischen beiden Funktionen richtig wiedergegeben, wie z. B. 
auch das von Karman durchgerechnete Beispiel eines Luftschiffkérpers zeigt. Die Nullstellen 
von f(t) fallen dann angenahert mit den Maximalstellen von F(t) zusammen, und diese Werte 
entsprechen ungefahr den Abszissen fiir die Extremwerte des Parallelkreishalbmessers des Dreh- 
kérpers. Bei einer kegelartigen Spitze am hinteren Ende hat z. B. f(t) eine einfache Wurzel 
und F(t) eine Doppelwurzel in t= — 1. Die Funktionen f und F verhalten sich also angendhert 
wie die Querkraftverteilung und die Momentenlinie eines stetig belasteten linearen Balkens. 

Die Berechtigung, die Gleichung (19) als erste Annaherung der Problemlésung zu betrachten, 
ergibt sich aus folgender Uberlegung. Zunachst kann man zeigen, daf die Quellsenkenbelegung 
f(t) gleichwertig ist mit einer gewissen Verteilung von axial gestellten Dipolen entsprechend der 
Funktion F, deren negative Ableitung mit f(t) ttbereinstimmt. Setzt man f(t) = — F’ (t), so 
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ergibt sich aus dem Potential der Quellbelegung durch partielle Integration 


Os 1 Vee dt 1 Ga 1 | Us (20) 


47 Tan iod len _ an r2 


2 et 2 . . . 
or — (x — t)?-+ 7? und cos } = (x — 2)/r zu setzen ist. Wenn, wie vorausgesetzt wird, F(t) 
in den Endpunkten 1 und 2 der Grundstrecke verschwindet, so wird 


Dee 1 bees ce (21) 


Am r’ 


Das ist aber der Form nach das Potential einer Verteilung von in die x-Richtung gestellten Di- 
polen mit dem auf die Langeneinheit bezogenen Moment F(t). Der einzelne Dipol vom Moment m 


m 
270, 
Bei Drehung der Dipolachse in die y-Richtung erhalt man bei seitlicher Paralleleteomoescinane 


gibt nun, mit dem Parallelstrom vereinigt, als Grenze eine Kugel vom Halbmesser R = 


m 


3 
digkeit v, in Richtung der negativen y-Achse wieder eine Kugel vom Halbmesser Vx . Man 
Op 


braucht also nur das Moment entsprechend abzuandern, um die gleiche Kugel zu gewinnen wie 
im ersten Fall. Bei einer Verdrehung aller Dipolelemente um 90° mit entsprechender Anderung 
der Momente wird im allgemeinen die Grenze nur in erster Annaherung erhalten, da in Wirk- 
lichkeit infolge gegenseitiger Beeinflussung eine Anderung der Geschwindigkeitsvektoren und 
damit eine Verschiebung der Grenze eintritt. DaB aber diese Verschiebung nur gering ist, ergibt 
sich aus der Eigenschaft des Ellipsoids, bei dem bereits die als erste Naherung bezeichnete An- 
nahme von F(t) mit der genauen Lisung zusammenfallt. 

Wir wollen zunachst, d. h. vor der Verfolgung weiterer Naherungslésungen, noch einige Fol- 
gerungen aus der Beziehung f= — c F(t) ableiten. Wenn man die Reihen fiir f und F einfihrt, 
so ergibt sich aus 


f= > P= — eF')=— ch, tae DBP, 


daB die a, und /, in einfachster Weise wie folgt voneinander abhangen: 


evy(y+ 1) $,=ch,—a,. (22) 


Dann wird aber 
A a Py 
P— = >'f, P,(u) 0, (6). (23) 
ferner das Potential der Querstrémung fiir y = 0 
B — 
= — = — 2 
®,= — 5 DF PQ. (24) 
Aus der Gleichung fiir die Kontur R 


oe > bP, (0G) 770 


folgt aber 
DP Oi (25) 
daher ist in der Meridianebene y — 0 E 
5 ea 


Die Potentialwerte ®, fiir die Punkte der Meridiankurve sind also den Ordinaten von R pro- 
portional. Dieser Satz ist besonders fiir eine erste Abschatzung des Tragheitskoeffizienten k, 
von Wichtigkeit, der sich unmittelbar aus dem Quotienten 


_ _ f ®,n dx a7 
fa Vo fn? dx Gy) 


berechnen 148t1. In erster Annaherung wiirde sich dann ergeben 


ky = Bee (28) 


1 W. Miiller, Ing.-Archiv, 18 (1950), S. 341. 
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Das Verhiiltnis der Gré®en A und B hangt von den Abmessungen des Rotationskérpers ab. 
Beim Ellipsoid, wo diese Verhaltnisse genau zutreffen, hat man c= 1/2 ferner 


A 62 — ] B Gr II 
oy = a P) = = o____, (29) 
2a" QilEo) ee ee easy 
So 
daher 
i= ee: (30) 
Z— Q1(60) 
0 


Zur Aufstellung weiterer Naherungen im allgemeinen Fall kann man eine erweiterte Verteilungs- 
funktion mit willkiirlichen Konstanten zugrunde legen, die dann nach einem noch zu besprechen- 
wie-1) den Verfahren zu ermitteln sind. 
Wenn wir etwa den aus der qua- 
dratischen Quellsenkenverteilung 


y F(t) = (¢-F 1) (¢— 1/3) 

gewonnenen R, ausgehen, so wird 
3 die erste Annaherung der zugehéri- 
gen Dipolverteilung durch die Funk- 


2 tion 


Fi =— 3, f[ f(t) dt | 

31 
=yet— Hey) © 
Af dargestellt. In zweiter Anndherung 


kann man dann eine Funktion hohe- 
ren Grades einfiihren, die in t= 1] 


tov 


-02 
einfach und in t= — 1 doppelt ver- 
ty schwindet, im iibrigen aber Kon- 
j | stante enthalt, die den besonderen 
Bedingungen des Problems anzu- 
sad passen sind. Wir setzen etwa 
-08 | Py (1) Ty (e -1) | 
Te 2Ae lp aoaleaie3 Vt va 
rH + ya (t+ 14]. 
a 6, Anpassung der Relativstré- 


Abb. 3. tg t und €y-Kurve tur die Punkte der R,. mung an die gegebene Kontur, Die 
vervollstandigte Lésung kann nun 
dadurch erreicht werden, daB die Koeffizienten /, unserer Entwicklung der besonderen Form- 
gebung des vorliegenden Rotationskérpers R, angepaBt werden. Die Ebene y= 0 ist die 
Symmetrieebene fiir die Seitenstrémung, die von der Bewegung des Kérpers in der y-Richtung 
herrihrt, und zugleich die einzige Ebene, welche vollstandige Stromlinien, sowohl der Langs- 
wie auch der Querstrémung enthalt. Die notwendige und hinreichende Bedingung fiir die Iden- 
titat der Grenzen in beiden Teilstrémungen l48t sich dahin aussprechen, daB die Verhaltnisse 
der relativen Geschwindigkeitskomponenten in beiden Fallen in allen Punkten der Meridian- 
kontur iibereinstimmen. Fihren wir die Geschwindigkeitswerte nach (11) ein und bezeichnen 
wir den als bekannt vorauszusetzenden Neigungswinkel der Tangente mit t, so laBt sich die 
»langentialbedingung“ in der Form schreiben 


DB, H, (Cu) + tg tS B,C, (C uw) = 92 0, (C2 — 2). (33) 
Wenn man diese Bedingungsgleichung fiir m Profilpunkte in Ansatz bringt, so lassen sich die 
Koeffizienten b, = B fj,/2 1 v, fiir »= 1 bis m aus einem System von m linearen Gleichungen 
eindeutig bestimmen, und damit ist nicht nur die stetige Verteilung der Dipole, sondern auch 
das Potential der Querstrémung und wenn man von einem durch eine Quellverteilung erzeugten 
Drehkérper ausgegangen ist, auch das zusammengesetzte Stromfeld mitsamt den Druckverhalt- 
nissen als angenahert bekannt anzusehen. 
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Tabelle 2. Die zu R, gehérenden Funktionen Gn und Hn. 
ie ae 
g LU x” n tgt G, G, CG, G, A, H, H, A, 
a i a el el 
1,004 | —0,863 | —0,8662 | 0,0453 | 0,3083 0,2203 | 0,1355 0,0851 0,0108 | 0,7399 | —0,6337|  0,4860 | —0,3715 
1,008 | —0,703 | —0,7086 | 0,0901 | — 0,2610 0,3550 | —0,1283 | —0,0687 0,1632 | 1,448 | --1,047 0,5741 | —0,2146 
1,010 | —0,617 | —0,6233 | 0,1114 | — 0,2400 0,3821 | —0,0573 | —0,1670 0,2041 | 1,776 | —1,138 0,4420 | 00,0746 
1,015 | —0,389 | —0,3948 | 0,1600 |  0,1854 0,3302 | 0,1913 | —0,2975 0,0667 | 2,325 | —0,9613 | —0,1278!  0,4656 
1,020 | —0,142 | —0,1450 | 0,1985 |  0,1246 0,1393 | 0,3900 | —0,1575 | —0,2049 | 2,5865 | —0,3873 | —0,5650 | —0,2790 
1,023 | 0,0143 | 0,0147 | 0,2156 | 0,0871 | —0,0143 | —_0,4180 0,0164 0,2397 | 2,5101 0,0378 | —0,6375 | —0,0276 
1,028 | 0,288 0,2961 | 0,2281 | 0,0015 | —0,2643|  0,2713 0,2470 | —0,0420 | 2,0145 |  0,6137 ; —0,2308 | —0,3443 
1,030 | 0,402 0,4144 | 0,2255. | —0,0406 | —0,3380]  0,1560 0,2590 | 0,07148 | 1,7030 | 0,7247 | —0,0544 | —0,3048 
1,036 0,760 | 0,7874 | 0,1758 | —0,2635 | —0,3210 | —0,1587 | —0,0185 0,0616 | 0,5210 | 0,4250|  0,2800} —0,1310 
1,039 | 0,945 0,9840 | 0,0897 | —0,7714 | —0,0975 | —0,0825 | —0,0635 0,0461 | 0,0265 | 0,0363} 0,0493 | —0,0397 


Wir legen den Rotationskérper 
R, mit der oben angegebenen 
Quellkurve und dem Streckungs- 
verhaltnis av, :A= 6,67 zu- 
grunde, fiir den in der Abb.3 so- 
wohl die Abhangigkeit der Ellip- 
soid-Koordinaten € (die einer 
Geraden sehr nahe kommt), wie 
auch die Kurve tgt als Funktion 
von & dargestellt ist. Ferner 
haben wir in den Abb. 4 und 5 
den Verlauf der Funktionen G, 
nnd. Hetur y— I bis 74 in 
Abhangigkeit von € veranschau- 
licht und schlieBlich in Tabelle 2 
fiir einige Profilpunkte die Werte 


dieser Funktionen zusammen- 
gestellt. 
Um das Wesentliche des Ver- 
fahrens hervortreten zu lassen, 17th‘ «70 «7% ~—«0R~S«C«T OP SBSS«TOBS«SOSESS«S ON 
beschranken wir uns auf vier Abb. 4. Die der R, zugeordneten Funktionen G. 
geeignet gewahlte Punkte der 2ehtt 
Meridiankurve, die zu den Werten ~ 
c= 1,004; 1,008; 1,023; 1,036 é 
gehéren. Mit Benutzung der Ta- ee 
belle ergeben sich dann aus der 4 
Tangentialbedingung  folgende g 
vier Gleichungen zur Bestimmung 
der Groen 0b’ de - 
7, 
0,8078 b, — 0,6756 b, 12 
+ 0,5122 b; — 0,3682 b; 
= (MILD, 28 
1,541 6; — 1,0805 6, ty 
+ 0,5561 b, — 0,172 b, g4 
= 0,0469 , 
| (34) g 
2,509 b, + 0,0742 b, 0 
— 0,6361 b; — 0,0485 by My : 
1052250) -O4 
0,6056 bj + 0,4668 65 = Me , 
+ 0,285 b; — 0,1472 b, ¢ 
== 0,087 
“427g 7008 ‘tule ‘10% ‘102 +102" = 008 ~—«t2S TOON 
ie der R, zugeordneten Funktionen H. 
Daraus berechnen sich die Werte BP ery ie ope guingmecotanet cu Buskt 
/ oe E> / == 3 
bi = 0,08814, 6, = 0,07830, 6,—=—0,004198, 6, = 0,01163, 35 
b, — 0,04407, b= 0,01605, b,—=—0,0003498, 6, = 0,000582. f (62) 
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Durch Multiplikation mit v, erhalt man dann auch die GroBen B£, und damit die vollstandige 
Darstellung des Potentials. 


Um nun die Reduktion auf die Form 
Fy (t) = (1— ®) [y, (+ D+ ye (t+ 1)? +73 @+ YF 4+ m4 (t+ D4] (36) 


durchzufithren, haben wir nach Tabelle 1 zu setzen 


6 8 16 
bien tere 8 Slaneay sles 22 


WR ec 
bay | 3 V2 7 3 | a1 4’ 
2 


(37) 


Dann wird 
V1, = 0,04971, y.=0,0083, y,;——0,0208, y,s=—0,00772. 
Zum Vergleich mit der ersten Naherung K (1 — #?) (t+ 1) miissen wir beide Funktionen auf 


die gleiche Flache oder Gesamtergiebigkeit reduzieren. Man erhalt dann, wie man durch eine 
einfache Rechnung feststellt, K = b, = 0,0441. 


Abb. 6- Die Verteilungsfunktionen f(i) und F(t) der R, . 


Wir haben das Ergebnis der Rechnung in einer Tabelle und einer Zeichnung (Abb. 3) dar- 
gestellt, aus denen hervorgeht, dai die Abweichung der beiden Funktionen Fy und Fy vonein- 
ander héchstens 4 bis 5% des gréBten Wertes betragt, so daB auch damit die Berechtigung, 
F, als erste Naherung zu betrachten, bestatigt wird. Es ist ohne weiteres verstandlich, daB 
die Lésung durch Hinzunahme weiterer Punkte der Kontur beliebig verbessert werden kann, 
was bei entsprechender Vollstandigkeit der Tabellen fiir die Kugelfunktionen keine wesentliche 
Schwierigkeiten bereiten wiirde. 


t | —1 | —o8 | 0,6 | 0,4 | 0,2 | 0 | 02 | 04 | 06 | 08 | 1,0 


0,00319 
0,00310 


0,0113 
0,0128 


0,0339 | 0,0441 | 0,0508 
0,0360 | 0,0448 | 0,0495 


0,0515 
0,0490 


0,0450 
0,0430 


0,0286 
0,0293 


0 


Zusammenfassend ergibt sich also folgendes Ergebnis. Wenn f(t) und das Streckungsver- 
haltnis A: 7 v,, ferner die Geschwindigkeit v, gegeben sind, so kann man unmittelbar A und 
a, = Aq,/2 v, bestimmen. Ferner ergibt sich aus v, und den Koeffizienten b, die GréBe 


| 


| 
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Bb,=27v,b,. Man hat z. B. mit v,:0,—=tg yp 


Damit sind dann alle GréSen bekannt, welche fiir die vollstindige Berechnung der aus der Langs- 
und Querbewegung entstehenden zusammengesetzten Strémung erforderlich sind. 


7, Die Differentialgleichung der Kontur bei der Querstrémung. Bei gegebener Dipolverteilung 
ist, wie ausgeftihrt, das System der Potentialflachen bekannt. Fiir die Relativstrémung hat man 


; = cos p >’, Pius) Q4(C) — vy = konst. (38) 


Die Stromlinien sind réumlich doppelt gekriimmte Kurven, fiir die zunachst keine Darstellungs- 
modglichkeit vorhanden ist. Nur in der Meridianebene y = 0 sind die Stromlinien in ihrem ganzen 
Verlauf ebene Kurven, und wenn in dieser Ebene die orthogonalen Trajektorien fiir das System 
der Potentialschnitte auffindbar ist, so ist damit auch der Meridianschnitt des Rotationskérpers 
als singulare Stromlinie als bekannt anzusehen. 

Wir kénnen daher jedenfalls die Differentialgleichung fiir dieses Kurvensystem angeben. 
Zu diesem Zweck miissen wir die Bedingung fiir die Orthogonalitaét zweier sich im Punkte ¢ us 
schneidender Kurven aufstellen. Die Tangentenrichtung in einem Punkt der Meridiankurve R ist 


dye le ee (GoD) 
dx y(E Fue) (39) 


wenn = ¢’ gesetzt wird. Daraus ergibt sich dann fiir zwei Richtungen C; und €3 als Bedingung 
LL 


des Senkrecht-Stehens 
CREST Ge iG a1o=— 


eine Beziehung, die leicht am System der Potential- und Stromlinien im achsensymmetrischen 
Fall der Langsstrémung nachgepriift werden kann. Wenn man nun die Gleichung (38) der Po- 
B 


Ga" 


Le (40) 


tentialschnitte in der xy-Ebene nach yw differentiiert, so ergibt sich mit eee =D 
i aap 
Dat wl 
i oH eae 41 
oe rhe dQ eg Pere ai ( ) 
b, P! i Le | Tanke 
pede CANT 
Dann gilt fiir die Orthogonalkurven, wenn wir statt C3 wieder C’ schreiben, 
dQ? 
2 pe oie 
eGR aay oe raat oy ie 
(iS eee eT pil — ( ) 
L—p) >) br * + uy 
(Uy peace 
oder bei Benutzung der Beziehungen P} = P, je? 5,02 — 0, yc2—1 und 
dP! 
ey Pay (Vee, eS. 
Me (43) 
Qu 


fea — 
VP—1 =O SGa= 0 +NAZ—6@ 
fir die auf die Kugelfunktionen P, und Q, reduzierte Form der Differentialgleichung 
d Fe dM nag te 
Db re YZ (Pe) QI — Ds Pr) GC) gE) = gy CH) (44) 


Fir den Fall des aus dem einfachen Quellsenkenkeil ableitbaren Ellipsoids ist die Losung 
leicht angebbar. Die auftretende Reihe beschrankt sich namlich auf ein einziges Glied. Wir 


haben also 


22 )— FEM +H OI=0. (45) 
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Setzt man b OQ; +1=C, also 26 4Q0,—Cuc—90, 
so wird auBer u = 0 

26,Q,— 6,60,—¢C¢=0 
und daher 


7 
= Ae ee at 46 
b= 50.) — bo Or (Co) io 


Benutzt man schlieBlich die Beziehungen 


La 20, , O. 
Oe Aaa = ~ a, ’ (47) 


so folgt wieder wie oben 
pee ey (48) 
ae Q: (60) 
Auf die Konstruktion der Lésung der Differentialgleichung im allgemeinen Fall hoffe ich an 
anderer Stelle zuriickzukommen. 


(Eingegangen am 15. September 1951.) 
Anschrift des Verfassers: Prof. Dr. W. Miiller, Miinchen 9, Taubenstr. 41. 


D 57 275 4022 0,075 (K.B./Z. 040). 


‘ 


Verstérker und Empfiinger 


Bearbeitet von Dr. eetha: Dr.-Ing. eh, M. J. 0. Strutt, ord. Prof. fiir theoretische Elektro- 
technik an der eidgenéssischen technischen Hochschule Ziirich. (Lehrbuch der drahtlosen 
Nachrichtentechnik, herausgegeben von Nicolai v. Korshenewsky, Stockholm, und Wilhelm 
T. Runge, Berlin Bd. IV.) Zweite, verbesserte Auflage. Mit 425 Abbildungen. XV, 
422 Seiten. 1951. - : | is) Ganzleinen DM 46.50 


nN 


Aus den Besprechungen:....Das vorliegende, in zweiter Auflage erscheinende Lehrbuch behandelt di 
grundlegenden physikalischen Eigenschaften sowie die Berechnung von Vorstarkeimn und Empfangern. Cogontibey 
der ersten Auflage wurden verschiedene kleinere und gréBere Textanderungen vorgenommen.. Dem Charakter eines 
Lehrbuches entsprechend, legte der Verfasser jedoch den gréBten Wert auf die eingehende Behandlung der Grund- 
lagen ... Der vorliegende Band kann Studenten zur Einfiihrung und dem praktisch tatigen Ingenieur als Nach- 
schlagewerk bestens empfohlen werden. » Technische Mitteilungen‘ 


Einfihrung in die theoretische Elektrotechnik. Von Dr.-Ing. Eb. Kat 


Kiipfmiller, Hon.-Professor an der Technischen Hochschule Stuttgart. Vorstandsmitglied 
der Standard Elektr.-Ges. Vierte, verbesserte und erweiterie Auflage. Mit 474 Ab- 
bildungen. V, 441 Seiten. 1952. ’ Ganzleinen DM 27.60 


... An Kipfmiillers Theoretischer Elektrotechnik schatzen wir Ingenieure, daB sie aus der Ingenieurtatigkeit heraus. 


fiir die technische Anwendung geschrieben ist, um mit Hilfe der Theorie technische Aufgaben in Entwurf, Konstruk- 
tion und Herstellung mit einem Minimum an Aufwand zu lésen. Auch seine Auffassung Schwachstrom- und Stark- 
stromtechnik einheitlich in gemeinsamer Sprache darzustellen, ist ein erfreulicher Fortschritt. |. 

- Fi »Hlektrotechnik und Maschinenbau‘ 


* 


Stromrichteranlagen der Starkstromtechnik. cintinrung in theorie 
und Praxis. Von Dr.-Ing. Helmut Anschitz, Mit 169 Abbildungen. VII, 213 Seiten. 1951. 
. Ganzleinen DM 18.— 


Aus den Besprechungen: Die Technik der Stromrichter zur’ Umformung elektrischer Energie von einer 
Stromart in die andere hat sich im Laufe der letzten Jahre besonders auf'dem Gebiete der Starkstromtechnik erheblich 
.erweitert. Sie beginnt u. a. auch im Rahmen der Energietibertragung durch hochgespannten Gleichstrom vermehrte 
Bedeutung zu erhalten. Andererseits harren besonders auf diesem Gebiet noch viele Aufgaben der Lésung, Der Ver- 


' fasser legt die sich daraus ergebenden Richtlinien in recht iibersichtlicher und theoretisch gut fundierter Weise dar. 


Er gibt einleitend eine Ubersicht iiber die verschiedenen Stromrichterarten und deren Aufgaben. In einem 2, Ab- 
schnitt gelangen die Gleichrichter, ihre Theorie und Berechnung bis in alle Details zur Darstellung. Ein weiteres 
Kapitel ist den Wechselrichtern zur Umformung des Gleichstromes in Wechsel- oder Drehstrom und die Umrichter 
zor Umformung von Wecehsel- oder Drehstrom in solchen yon anderer Frequenz gewidmet. Alsdann folgen die 
theoretischen Grundlagen und Erlauterungen zum Priifen und Messen yon Stromrichtern und Stromrichteranlagen. 
Weitere Abschnitte orientieren uber die Schalt- und Schutzeinrichtungen suwie tiber die Projektierung von Strom- 
richteranlagen fir die Starkstromtechnik. Einige spezieil die Praktiker interessierende Angaben und Hinweise zur 
Anwendung des Stromrichters beschlieBen diesen wohlformulierten Beitrag; daraus seien erwaihnt deren Gebrauch 
fir Licht- und Kraftnetze, Bahnanlagen, fiir elektrochemische Werke, Sendeanlagen, Ladestationen fiir Akkumulae 
toren usw. Das Werk schlieBt ab mit einem Literatur-, einem Formel- und einem Sachverzeichnis. 


Zum Verstandnis dieses Buches sind gute Vorkenntnisse der gesamten Elektrotechnik notwendig; es hilft dem Stu- 
dicrenden und dem in der Praxis stehenden Ingenieur zum raschen und griindlichen Kinarbeiten in das Spezialgebiet 


’ der Stromrichtertechnik und -theorie, Zu diesem Zwecke eiguet es sich in hervorragender Weise und kann zur Au- 


schaffung empfohlen werden. } = » Technische Rundschau* 


Stahlmaste fir Starkstrom-Freileitungen. Bercchnung und Beispiele von 


Wilhelm Taenzer. Zweite, neubearbeitete Auflage. Mit 264 Abbildungen. LV, 98 Seiten. 
1952. | DM 21.— 


... Das’ vorliegende Werk, aus der Praxis und fiir sie geschrieben, bildet einen der bisher wertvollsten Beitriige zur 
Férderung wirtschaftlicher Stahlbauten und kann zur lohnenden Anschaffung als ebenso wertvolles Nachschlagewerk 
wie als unentbehrliches Riistzeug zur Vertiefung in das Sondergebiet des Stahlmastenbaues und zur Schulung des 
technischen Nachwuchses angelegentlich empfohlen werden... Stahl und Hisen* 


Pa a eS Ee eT ay AMT enn Nee 
SPRINGER-VERLAG / BERLIN. GOTTINGEN- HEIDELBERG 


Maschineneleniente. Entwerfen, Berechnen und Gestalten im Maschi- 
nenbau. Ein Lehr- und Arbeitsbuch, Von Dr.-Ing. G. Niemann, Professor an der Tech- 


nischen Hochschule Braunschweig. 3 


Erster Band: Grundlagen, Verbindungen, Lager, Wellen und Zubehor. 
Mit 795 Abbildungen. VIII, 308 Seiten. 1950. Ganzleinen DM 28.50 
’ / J 4 
Aus den Besprechungen: Das Feblen eines modernen Buches tiber Maschinenelemente wurde seit langem 
als empfindlicher Mangel empfunden. Das Erscheinen des 1. Bandes Niemann darf deshalb mit Freuden vermerkt 
werden um so mehr, als bereits der Name des Verfassers fiir die Giite des Werkes garantiert. Die Gliederung des 
Buches lehnt sich an die bekannte~und tibliche Einteilung an. Nach einer Einfihrung in die Methodik des Kon- 
struierens werden Gestaltungsregeln, Festigkeitserechnung und Werkstoffe behandelt, wobei dem so wichtigen Leichtbau 
ein besonderes Kapitel gewidmet ist. Von den eigentlichen Maschinenelementen werden im 1. Band besprochen: 
Verbindungselemente, Lager, Wellen und Achsen sowie die dazugehérigen Wellenverbindungen. Ein ausfihrliches 
Verzeichnis zu jedem Abschnitt verweist auf das gerade fiir dieses Gebiet manchmal sehr verstreute Schrifttum. 
Es darf als ein Vorteil des Buches angesehen*werden, da8 es nicht nur ein Lehrbuch ist sondern — von einem er- 
fahrenen Praktiker geschrieben — auch stark auf die Bediirfnisse der Praxis abgestellt ist. Von irgendwelchen 
Ableitungen wurde nur sparsam Gebrauch gemacht und nur da, wo sie fiir das Verstindnis unumganglich nétig 
sind; dagegen werden die zahlreichen Erfahrungsangaben und Zahlenunterlagen dankbar aufgenommen werden. 
DaB8 bei dem umfangreichen Stoff manches nur kurz behandelt werden konnte, ist verstindlich. Der unter diesem 
Zwang oft 2um Telegrammstil gekiirzte Text wird vom vielbeschaftigten Fachmann sicher begriiBt werden, wenn er 
auch dem Lernenden einige Schwierigkeiten bereiten kann. Diesem werden aber die zahlreichen, durchgerechneten 
Beispiele wieder das Eindringen in den Stoff erleichtern. .. 


Das Buch kann als modernes Standardwerk der Maschinenelemente sowohl dem Studierenden als auch dem in der 
Praxis stehenden Ingenieur auf beste empfohlen werden. » Zeitschrift des Vereins Deutscher Ingenieure.* 


« 


Berechnung der Maschinenelemente. von Dipt.-ing. M. ten Bosch f, Pro- 


der ,, Vorlesungen iiber Maschinenelemente“*. Mit 926 Textabbildungen. X, 534 Seiten. 
1951. Ganzleinen DM 45.— 


ae 


AusdenBesprechungen: Dieses moderne Lehrbuch der Maschinenelemente vermag auch dem Bauingenieur, 
insbesondere auf dem Gebiete des Stahlbaues, vielerlei Aufklarungen und Unterrichtung zu bieten, hat doch roh 


Angewandte Festigkeitslehre, Verbindungen und Wellen bilden wertvolle Erganzungen zu den bekannten Fach- 
biichern der Bauingenieure, besonders wertvoll nicht zuletzt deshalb, weil viele Fragen von etwas anderem Gesichts- 
punkt aus betrachtet werden, als wir Bauingenieure dies gewohnt sind. Dieses Buch verdient auch yon Bauingenieuren: 
beachtet zu werden. Der Preis ist miSig zu nennen in Hinsicht auf den groBen Umfang und die vorziigliche Aus- 
stattung des Buches. ‘ 4 »»Der_ Bauingenieur.* 


Py 


j 
; 


e 


Maschinenelemente. Leitfaden zur Berechnung und Konstruktion fiir Maschinen- 
bauschulen und fiir die Praxis mittlerer Techniker. Von Dipl.-Ing. W. Tochtermann, Pro- 
fessor an der Staatl. Ing.-Schule EGfingen. Seehste, véllig neubearbeitete Auflage. 

_Mit 641 Abbildungen, XII, 515 Seiten. 1951. Ganzleinen DM 34.50 


Au’s den Besprechungen: Das véllig neubearbeitete und ausgezeichnet zusammengestellte, Werk ist als 
Lehrbuch fiir technische Lehranstalten und als Handbuch fiir den in der Praxis stehenden Techniker bestimmt, 


Der Verfasser stellt das Wesentliche der Ingenieurausbildung in den Vordergrund, nimlich die Maschinenelemente 
und die in ihnen ruhenden Probleme. Ihre Berechnung und zeichnerische Darstellung werden hier unter: Weglassung 
veralteter Konstruktionen wirklich erschépfend behandelt, wobei allerdings vorausgesetzt wird, daB der Studierende 
die Grundlagen der Mathematik, Mechanik und Festigkeitslehre beherrscht. Maschinenclenente bilden das Funda- 
ment jeder Konstruktion, der zukiinftige Techniker muB sie richtig berechnen und gestalten kénnen, nicht nur in 
bezug auf Sicherheit, sondern auch nach Gesichtspunkten rationeller Herstellung. Aus diesem Grunde legte der 
Autor bei der Neubearbeitung des Buches besonderen Wert auf eine Vielzahl aus der Praxis gegriffener Berechnungs- 
beispiele und auf eine Fille aufschluBreicher Abbildungen, die das Gesagte klar und deutlich veranschaulichen, , , 


j i »»Maschine und Werkzeug,“* 
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Diesem Heft liegen 3 Beilagen des Springer-Verlages, Berlin - Géttingen . Heidelberg, bei. 


es 


Springer-Verlag / Berlin Géttingen. Heidelberg. Printed in Germany, D 57 275 4022 0,8 K. B./Z. 040 


fessor an der Hidgenéssischen Technischen Hochschule Ziirich. Dritte,erganzte Auflage 


gerechnet etwa die Hialfte des Buches auch fir ihn unmittelbar Bedeutung. Besonders die drei groBen Kapitel: 
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